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ABSTRACT 


Complementary  binary  sequences  were  invented  by  Golay  in  an  in¬ 
vestigation  of  infra-red  multi- slit  spectrometry.  This  dissertation 
formalizes  the  basic  results  obtained  by  Golay  and  develops  new  concepts 
and  techniques  for  examining  the  characteristics  of  these  special  binary 
codes.  This  work  has  developed  new  understanding  of  the  structure  and 
methods  for  the  decomposition  of  complementary  sequences. 

Complementary  sequences  have  the  property  of  an  infinite  correla¬ 
tion  peak  to  ambiguity  ratio  when  detected  with  a  matched  filter.  These 
binary  sequences  should  find  much  application  as  pseudo  -  random  noise 
modulation  signals  for  both  radar  and  communications  systems. 

A  discussion  of  the  need  for  such  sequences  is  included  in  the  intro¬ 
duction  and  is  followed  by  a  state  of  the  art  description,,  An  operations 
group  on  the  sequences  is  formulated  and  the  proofs  of  several  theorems 
concerning  the  operations  group  are  given  in  a  rigorous  manner.  One 
reason  for  developing  the  operations  group  is  the  application  to  elimina¬ 
tion  of  redundance  in  the  computer  search  for  new  codes. 

Several  invariant  properties  of  complementary  codes  are  proved 
through  the  use  of  the  Hamming  distance  concept.  Many  more  invariant 
properties  of  the  sequences  are  demonstrated  through  the  introduction  of 
a  Hamming  vector  The  concept  of  a  Hamming  vector  is  extremely  use¬ 
ful  as  a  complementary  code  decomposition  tool.  A  large  number  of 
theorems  are  proved  to  enhance  its  use  in  this  field. 

Several  computer  searches  for  complementary  sequences  are 
described  and  the  actual  computer  programs  for  the  CDC  1604  are  in¬ 
cluded  in  the  Appendix. 
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CHAPTER  I 
INTRODUCTION 


Search  radars  of  today  are  required  to  obtain  smaller  and  smaller 
targets  at  ever  increasing  ranges.  This  has  led  to  many  different  tech¬ 
niques  of  summing  radio  frequency  energy  to  obtain  these  long  ranges. 
Complementary  sequences  are  well  adapted  for  coding  the  RF  energy 
pulses  to  give  a  summing  of  the  returning  RF  energy  in  the  echo.  The 
present  work  is  more  concerned  with  the  characteristics  of  complemen¬ 
tary  codes  such  as  their  operation  groups,  Hamming  distances  and 
Hamming  vectors,  than  it  is  with  their  application  as  modulating  signals 
for  RF  carriers.  However,  the  introduction  will  develop  the  radar 
problem  as  a  vehicle  to  emphasize  the  use  of  these  codes  in  future  appli¬ 
cations  in  radar.  Their  application  will  probably  be  equally  important 
in  the  communications  field.  The  initial  work  in  complementary  sequen¬ 
ces  was  directed  toward  an  infra-red  application,  but  it  is  felt  that  the 

eventual  application  of  these  sequences  will  spread  to  many  different 

1  Z 

fields.  * 

The  usual  approach  to  the  problem  of  increasing  radar  range  is  to 
look  for  ways  of  varying  some  of  the  parameters  in  the  radar  range 
equation  which  will  lead  to  a  greater  maximum  range. 

The  radar  range  equation  is 
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where  P  =  peak  transmitted  power  in  watts, 

P  =  minimum  peak  detectable  signal  in  watts, 
mm 

S  =  scattering  cross  section  of  target  in  units 
consistent  with  range, 

G  =  gain  of  transmitting  antenna. 


A  =  effective  area  of  receiving  antenna  in 
consistent  units. 
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D  =  a  composite  loss  factor  for  transmission  lines, 
atmospheric  losses,  etc* 

The  parameters  under  the  control  of  the  design  engineer  are  A^,G, 

P  and  P^.  A  and  G  are  fixed  by  various  considerations  such  as 
min  t  e 

antenna  beamwidth,  side  lobe  suppression,  physical  size  of  the  antenna, 

etc.  ,  and  are  not  truly  available  as  adjustable  parameters  beyond  small 

variations.  Therefore  the  two  parameters  for  the  designers  to  optimize 

are  P  .  and  P  - 
mm  t 

P  .  is  basically  limited  by  three  factors:  Johnson  noise  in  the 
mm  7 

input  circuits,  shot  effect  and  other  noises  in  the  first  tube,  and  cosmic 
noise  picked  up  by  the  antenna.  P  ^  approaches  KBT^  in  a  perfect 
receiver,  where  K  is  Boltzmann’s  constant,  B  is  the  bandwidth  and  T 

a 

is  the  absolute  space  temperature.  P^  n  is  currently  being  reduced 
toward  the  above  value  through  the  utilization  of  parametric  amplifiers, 
masers  and  other  low  noise  devices,.  Another  method  of  reducing  P 

mm 

is  through  pulse  integration,  since  noise  voltage  averaged  over  a  period 
of  time  has  a  value  approaching  zero  while  the  signal  values  are  additive. 
However,  in  most  cases  this  integration  of  pulses  can  also  be  considered 
as  an  increase  in  transmitted  energy.  It  will  be  considered  from  that 
viewpoint  in  this  paper*  There  are  two  basic  types  of  integration,  co¬ 
herent  and  noncoherent.  Coherent  integration  gives  a  gain  of  approxi¬ 
mately  N,  where  N  is  the  number  of  pulses  integrated  and  noncoherent 

integration  gives  a  gain  of  approximately  N*  Both  of  these  values  are 

3 

compared  to  a  single  pulse  .  The  values  are  respectively  the  upper 
bound  for  coherent  integration  and  the  lower  bound  for  noncoherent 
integration. 

There  are  basically  three  different  integration  time  bases  possible 
in  a  radar  system.  One  is  the  antenna  scan  to  antenna  scan  where  the 
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operator  or  a  computer,  decides  if  the  echoes  from  previous  scans  were 
signals  or  noise.  The  second  method  of  integration  is  the  interpulse 
method  where  the  cathode  ray  tube,  a  delay  line,  or  a  computer  sums  up 
the  hits  within  one  antenna  scan  of  the  target.  The  third  method  is  intra¬ 
pulse  integration,  or  pulse  compression,  where  different  portions  of 
one  radar  transmitted  pulse  are  summed  up  Complementary  codes  are 
used  in  one  scheme  of  intrapulse  integration. 

An  examination  of  how  the  transmitted  power  or  energy  since 
Ej_  =  Pj_  A  T,  can  be  varied  is  now  in  order.  From  this  relation,  it  is 
seen  that  there  are  two  basic  methods  for  increasing  the  transmitted 
energy:  either  increase  peak  pulse  power,  or  increase  the  pulse  width. 

A  common  method  of  increasing  radiated  energy  has  been  to  increase 
the  peak  power.  This  method  ha  s  proven  both  expensive  and  wasteful  as 
many  radar  transmitter  tubes  are  peak  power  limited  while  only  dissipa¬ 
ting  a  very  small  percentage  of  their  average  power  capability,  and 
radar  modulators  are  taxed  with  a  very  similar  problem.  The  alterna¬ 
tive  did  not  look  very  promising  at  first,  since  to  increase  pulse  duration 
normally  reduces  the  high  frequency  content  of  a  pulse,  which  in  turn 
reduces  range  resolution  capabilities.  However,  intrapulse  modulation 
retains  the  high  frequency  content  of  a  narrow  pulse,  while  increasing 
the  energy  content  to  that  of  a  pulse  of  long  duration.  Intrapulse  modula¬ 
tion  allows  the  summing,  required  in  intrapulse  integration,  for 
detection  in  a  radar  receiver. 

There  are  two  methods  of  intrapulse  modulation,  the  digital  method 

and  the  analog  method.  Although  the  emphasis  in  this  paper  is  on 

digital  modulation  schemes  -  the  analog  system  will  now  be  briefly 

mentioned  since  it  is  a  part  of  the  overall  picture.  Figure  1.  1  contains 

4 

the  graphs  used  in  the  explanation  of  the  chirp  radar  scheme  In 
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Frequency  modulation  of  a  chirp  radar* 
FIGURE  1.  1 
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Signal  addition  through  the  use  of  a  delay  line. 
FIGURE  1.2 
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this  analog  method  of  intrapulse  modulation  the  frequency  of  the  outgoing 

pulse  is  swept  linearly  upward  with  time  as  is  shown  in  Figure  1.  1-b. 
This  gives  a  large  high  frequency  content  to  the  pulse  while  at  the  same 
time  maintaining  a  long  pulse  duration.  Upon  reception  the  signal  is 
passed  through  a  reversed  time  delay,  delaying  the  low  frequencies  more 
than  the  highs,  as  is  shown  in  Figure  1.  1-d*  This  allows  all  the  fre¬ 
quency  components  of  energy  to  arrive  at  the  output  simultaneously, 
giving  a  large  and  narrow  output  pulse.  The  type  of  filter  used  in  the 
delaying  process  is  known  as  a  matched  filter  and  is  discussed  in 
Appendix  I,  The  procedure  here  is  seen  to  be  quite  simple:  first  we 
generate  energy  components  with  various  frequencies  at  different  times, 

<w 

and  then  we  delay  all  components  so  that  they  arrive  at  the  output 
s  imultaneous  ly. 

As  a  preliminary  to  the  digital  type  intrapulse  modulation  consider 
a  transmitted  waveform  such  as  in  Figure  1.  2.  This  waveform  con¬ 
sists  of  three  pulses  of  unit  length  with  spacing  between  them  of  one  and 
three  time  units  respectively.  This  waveform  is  then  fed  into  a  delay 
line  in  the  receiver  which  has  the  delays  shown  in  Figure  1.2-b.  The 
output  is  either  a  one  or  a  zero  except  at  the  exact  match  of  signal  to 
filter;  the  output  is  then  three  as  is  shown  in  Figure  1.  2-c.  The  three 
pulses  have  therefore  been  summed  as  though  they  were  one,  at  the 
time  between  6  and  7  units.  These  pulses  still  contain  all  the  high 
frequency  components  of  the  short  pulse  of  unit  length  but  have  effec¬ 
tively  three  times  the  energy  transmitted  in  a  single  pulse  at  the  time 
of  exact  match.  There  are  two  features  to  be  noted  about  this  output, 
first  it  has  been  delayed  in  time  and  second  there  are  secondary  pulses 
at  various  undesired  times  which  might  be  classified  as  coding  noise. 
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FIGURE  1.3 


Fig.  1.3-  Matched  filter  detection  of  a  random  code. 

Fig.  1.4  -  Matched  filter  detection  of  an  optimal  noise  code. 
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I  will  consider  any  undesired  signal  which  might  be  mistaken  for  a  true 
signal  as  noise.  These  two  features  will  be  discussed  in  some  detail 
in  the  true  digital  intrapulse  modulation  schemes  to  be  discussed  later. 

The  scheme  just  considered  is  somewhat  inefficient  since  seven 
units  of  time  are  required  to  radiate  three  units  of  energy.  If  in  the 
transmitted  signal  all  the  spaces  with  zero  energy  output  in  the  pulse 
envelope  had  also  had  an  output,  seven  units  of  energy  would  have  been 
radiated  in  seven  units  of  time.  One  might  use  some  sort  of  frequency 
shift  or  phase  modulation  to  accomplish  this  purpose,  still  utilizing  of 
course  the  same  output  tube*  Consider  the  same  waveform  as  before 
with  the  units  of  energy  coded  as  plus  and  minus  by  the  modulator  by 
means  of  phase  shift  modulation  of  either  zero  or  180°.  The  waveform 
is  as  shown  in  Figure  1. 3-a.  Upon  reception  this  waveform  is  now  fed 
into  a  delay  line  in  the  receiver  as  shown  in  Figure  1*  3-b.  MAM  repre¬ 
sents  a  straight- through  amplifier*  Figure  1*  3-c  shows  the  output  of 
the  delay  line  summer  with  respect  to  time.  The  output  energy  is  now 
seven  times  that  of  a  single  pulse  and  has  the  frequency  content,  and 
therefore  the  range  resolution,  of  one  of  the  single  pulses.  A  time  delay 
in  the  output  is  again  apparent,  as  are  a  number  of  noisy  sub- peaks. 

This  was  a  code  picked  at  random  and  the  results  were  good.  There  are 
certain  classes  of  codes  which  are  optimal,  having  even  smaller  noisy 

sub-peaks.  These  codes  are  of  length  2n- 1  and  have  been  studied 
5  6  T 

extensively.  “  *  *  They  are  known  as  L  codes  or  pseudo- rando m  noise 

codes.  An  example  of  one  of  these  pseudo- random  noise  codes  is 
-+--+  +  +,  where  the  minuses  and  pluses  indicate  the  phase  of  the  RF 
signal.  When  this  particular  noise  code  is  impressed  upon  its  matched 
filter,  Figure  1.4-b,  the  output  of  this  summed  delay  line  is  as  shown 
in  Figure  1. 4-c.  It  is  to  be  noted  that  the  maximum  height  of  the 
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Matched  filter  detection  of  a  complementary 
code  of  length  8. 

FIGURE  1.5 
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"hash"  (sub-peaks)  is  one,  with  a  main  peak  of  seven,  whereas  in  the 
previous  example  picked  at  random  the  hash  had  a  height  of  two. 

Elaspas  has  made  a  study  and  design  of  a  radar  utilizing  this  form  of 

g 

modulation  and  matched  filter  detection.  He  infers  that  codes  of 
reasonable  length  should  be  used.  For  example,  a  code  of  length  1Z7 
would  be  a  reasonable  length  to  obtain  the  benefits  of  a  ps eudo- random 
matched  filte  r  radar.  An  optimum  code  of  this  length  has  a  maximum 
subpeak  of  13,  which  gives  a  good  coding  signal-to-noise  ratio. 

There  is  a  way,  however,  to  completely  eliminate  the  hash,  al¬ 
though  this  required  the  transmission  of  two  separate  sequences.  Both 
of  these  sequences  could  be  transmitted  through  one  output  tube  with  a 

moderate  band  width.  These  codes  which  eliminate  the  hash  are  known 

9 

as  complementary  codes.  An  example  of  a  complementary  code  pair 
being  detected  by  a  matched  filter  pair  is  shown  in  Figure  1.  5.  The 
two  codes  being  transmitted  will  be  called  A  and  B  for  convenience, 
with  their  respective  matched  filters  also  being  labeled  in  a  similar 
fashion.  The  receiver  output,  after  the  complementary  codes  are  detec¬ 
ted  by  the  matched  filter  and  their  outputs  summed,  is  shown  in  Figure 
l.5-g.  This  demonstrates,  as  was  mentioned  before  ,  that  the  hash  level 
is  zero.  This  lack  of  hash  is  indicative  of  an  infinite  coding  signal-to- 
noise  ratio.  This  particular  feature  is  the  prime  advantage  of  complemen¬ 
tary  sequence  pairs  over  other  forms  of  noise  modulation  Figure  l„5-g 

shows  the  same  amount  of  delay  as  did  a.ll  the  other  forms  of  matched 
filter  detectors,  and  shows  that  the  sum  of  all  the  energy  sub-pulses  is 
incorporated  into  the  main  pulse  at  exact  match. 

Each  of  the  systems  described  in  this  introduction  has  had  the  char¬ 
acteristic  that  the  transmitted  signal  is  spread  over  a  wide  frequency 
spectrum.  This  has  in  itself  two  advantages:  One  is  that  enemy  detection 
of  the  signal  becomes  more  difficult,  and  the  other  is  that  even  after 
detection  the  jamming  of  a  wide  spectrum  noise- like  signal  becomes 
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extremely  difficult.  Since  the  characteristics  just  described  are  de¬ 
sirable  for  a  secure  communication  system  as  well  as  for  a  radar 
system,  all  of  the  advantages  that  accrue  for  the  complementary  se¬ 
quence  modulation  in  radar  would  apply  in  a  communications  scheme. 

The  initial  major  objectives  of  this  investigation  were  two  in 
numbe  r. 

1.  The  formulation  of  a,  set  of  operations  on  complementary 
sequences  which  form  a  group,  along  with  a  theoretical 
study  of  the  invariant  properties  of  such  a  group. 

2.  The  search  for  new  complementary  codes*  An  exhaustive 

computer  search  for  a  new  code  of  length  26  was  made 
as  well  as  a.  partial  search  for  codes  of  length  34. 

During  the  investigation  two  more  major  objectives  were  added 
to  the  work. 

These  we  re : 


3.  The  application  of  Hamming  distance  to  both  the  comple¬ 
mentary  sequences  and  their  group  formulation,,  ^ 

40  The  application  of  Hamming  vectors  to  complementary 

sequence  pairs  and  also  to  their  group  formulation* 

The  original  objectives  were  oriented  toward  the  generation  of 

new  codes,  while  the  last  two  objectives  were  oriented  toward  the 

decomposition  of  codes*  All  of  these  objectives  were  attained  and 

consititute  a  large  portion  of  the  rest  of  this  work*  The  next  chapter 

is  a  study  of  complementary  sequences  designed  to  bring  the  reader 

to  the  state  of  the  art,  and  contains  some  proofs  which  have  not  ap- 

1  2911  12 

peared  elsewhere*  1,1  1  Chapter  3  is  the  formulation  of  the 

operations  group  on  complementary  sequences  and,  along  with  several 
theorems  and  proofs  it  also  contains  the  identification  of  this  operations 


II 


group  with  a  known  group.  Chapter  3  as  written  requires  only  a  min¬ 
imum  of  group  theory  knowledge.  Chapters  4  and  5  have  to  do  with 
the  Hamming  distance  and  Hamming  vectors  of  complementary  sequen¬ 
ces  both  in  standard  form  and  in  group  form.  Chapter  6  is  a  short 
chapter  and  is  concerned  with  the  pmof  of  two  theorems  used  m  the 
searching  for  new  codes.  Chapters  7  a,nd  8  describe  the  computer 
search  for  new  kernels  of  length  26  and  34  respectively,  the  actual 
computer  programs  used  being  listed  in  the  appendicie  s*  The  last 
chapter,  Number  9,  contains  the  summary,  conclusions,  and  suggestions 
for  future  research  in  this  interesting  field. 
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CHAPTER  II 

COMPLEMENTARY  SEQUENCES 
This  chapter  is  presented  to  provide  a  background  in  complemen¬ 
tary  sequences  to  enable  the  reader  to  better  understand  the  material 
in  later  chapters.  Most  of  the  material  in  this  chapter  is  an  adaptation 

of  Golay1  s  “Complementary  Series"  with  expansions  and  deletions  to 

9 

fit  the  needs  of  this  paper.  The  formalization  of  some  proofs  and  the 
addition  of  new  proofs  given  here  have  not  appeared  in  print  elsewhere; 
however,  it  is  obvious  from  the  tone  of  Golayvs  work  that  he  was  aware 
of  these  proofs.  This  modification  of  Golay5  s  work  is  offered  for  com¬ 
pleteness  of  the  dissertation  and  not  as  new  work, 

In  the  introduction  a  pair  of  codes  of  length  n  were  discussed  which 
had  the  characteristic  that  when  detected  by  matched  filters,  the  sum  of 
the  two  filter  outputs  was  everywhere  zero,  except  where  the  signal 
patterns  had  zero  time  delay  with  reference  to  the  filter  patterns*  At 
the  zero  pattern  delay  time  the  output  was  2n,  leading  to  an  infinite 
coding  signal- to-noise  ratio*  A  study  of  the  necessary  and  sufficient 
conditions  for  this  to  be  true  will  now  be  considered. 

Assume  a  complementary  code  pair  (A,  B)  and  suppose  that  one  of 
the  pair,  i,  e.  A,  is  longer  than  the  other.  If  this  be  true  then  there  will 
be  an  output  from  the  A  filter  with  no  cancelling  output  from  the  B 
filter.  This  is  not  allowed;  therefore  the  A  and  B  codes  must  be  of 
equal  length.  Therefore  assume  a  code  A,  n  bits  in  length,  with 

elements  an  through  a  ,  where  each  a  is  either  +1  or  -1,  Also, 

I  °  n  i 

assume  a  second  code  B,  n  bits  in  length,  with  elements  b  similarly 
defined.  Detect  each  of  these  codes  with  a  matched  filter  (the  time 
inverse  of  the  code)  and  then  sum  the  outputs  of  the  A  filter  and  the 
B  filter.  At  time  1  (see  Figure  2*  1)  the  signals  and  b  have 
entered  the  filter  segments  a^  and  b  ^  respectively.  Since  at  this 
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a  . 

n 


a3a2al 


time 


bn'  '  '  *b3b2bl 


The  a  and  b  of  the  filter  are 
1  1 

either  ±  1  indicating  straight  thru 
amplifiers  or  inverters., 


Time  1  a  ana0a , 
4  o  2  1 


signal 


a  a  ,  a  filter 

n  n-  1  n- 2  1 


b3b2bl 


b  b 


signal 


n  n- 1  n- 2. 


b^  filter 


Time  2 


a4a3a2al 


a  a  ,  a  ^ 
n  n-J  n-2, 


b4b3b2bl 


b  b  ,  b  ^ 
n  n-  1  n-  2 


Time  r  a  .......  a , 

r  1 

an . *°an+l-r  nhi-(r4!)  .  w  * 


Matched  filter  detection  of  a  complementary 
pair  of  length  n0 

FIGURE  20  1 
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time  the  delay  between  signal  pattern  and  filter  pattern  is  not  zero,  the 

summed  output  voltage  must  be  zero.  This  gives  the  condition  a}an  + 

bjb^  =  0.  At  time  Z  (Figure  Z.  1)  the  signals  are  one  segment  farther 

into  the  filter  and  the  output  must  again  be  zero,  giving  aian  j  + 

a^a  +  b,  b  ,  +  b?b  =  0. 

Z  n  1  n-  1  ^  n 

At  time  3:  a,  a  0+a0a  ,+a0a  +b  b  ^+b^b  ,+b-b  =  0 

1  n-Z  Zn-1  3  n  1  n-Z  Zn-1  3  n 

At  time  r  (r  4  n) :  a,a  +a^a  ,  #  ,  *  .  , ,  ,  . 

x  1  1  n+l-r  Z  n+ 1  -  ( r  -1 }  r  „  .  .  +a  a  +b,b  ,,  +D  .  . 

r  n  1  nti - r 

+  b  b  =0 


r  n 

At  time  n:  a.a,  +  a0a0  +  ......  +  a  a  +b  ,  b  ,  +  ......+  b  b  =  2n 

1122  nnl  l  nn 

for  at  time  n  the  pattern  of  the  signal  must  exactly  match  the  pattern 
of  the  filter.  The  explicit  statement  of  the  conditions  given  above  for  all 
times  is 


F  = 
J 


i  =  n-j 

S 

i=  1 


a.a.  . 

i  i+J 


'f 


i=n-j 

E 

i=  1 


b.b  ,  • 
i  i+J 


=  0 

-  2n 


j  4  o 
j  =  o 


(2.1) 


An  examination  of  the  equation  at  time  1  (j=  n  -  1)  shows  the  following 
eight  possible  solutions* 


al 

a 

n 

bl 

b  n 

1 

1 

1 

-  1 

1 

1 

-  1 

1 

- 1 

-1 

1 

-  1 

- 1 

-1 

-  1 

] 

- 1 

1 

] 

1 

1 

-1 

1 

1 

- 1 

1 

-  1 

-  1 

1 

-1 

-1 

-  1 

These  possibilities  show  that  if  the  A  pair  are  alike  *  the  B  pair  must 
be  unlike,  and  if  the  A  pair  are  unlike,  the  B  pair  must  be  alike*  The 
equation  at  time  2(j~n-Z)  indicates  that  if  the  number  of  likes  (like  pairs, 
i.e.  a^and  a^  ^,)  of  the  A  code  and  the  number  of  unlikes  (unlike 
pairs)  of  the  B  code  are  equal,  and  if  the  number  of  unlikes  of  A 


equals  the  number  of  likes  of  B,  the  equation  is  satisfied.  Extending 
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this  reasoning  to  time  r  shows  again  that  if  the  likes  of  A  equal  the 
utilikes  of  B  and  if  the  unlikes  of  A  equal  the  likes  of  B  the  equation 
is  again  satisfied* 

It  will  now  be  proved  that  the  last  half  of  the  previous  statement  is 
not  necessary,  since  the  first  half  of  the  statement  contains  the  necessary 
and  sufficient  requirement. 

Let  =  number  of  unlike  pairs  in  A  for  the  spacing  specified. 

(j  in  equation  Zc  1) 

Let  L  =  number  of  like  pairs  in  A  for  the  spacing  specified, 

(j  in  equation  2e  1) 

Let  =  number  of  unlike  pairs  in  B  for  the  spacing  specified* 

(j  in  equation  Zu  1) 

Let  L^  =  number  of  like  pairs  in  B  for  the  spacing  specified. 

(j  in  equation  Z\  1) 

Assume  that  L  =  U  ,  we  want  to  prove  that  U  =  L,  0 
a  d  a  d 

Since  the  sum  of  the  number  of  like  pairs  at  each  spacing  with  the  num¬ 
ber  of  unlike  pairs  at  the  same  spacing  must  equal  the  total  number  of 
pos  sibilities  , 

+  L^  =  n-(n-r)  for  r  4  n  where  n  is  the  length  of  the  code 

and  (n-r)  is  the  number  of  units  of  time  delay 
from  exact  match* 

For  the  same  reason 

Ub  +  1^  r  n-(n-r)* 

But  since  U,  =  L  , 
b  a 

U  +  U,  =  n-  (n-  r  )o 

a  b  '  9 

Therefore  U  =  L.  . 

a  b 

The  definition  of  a  complementary  pair  of  sequences  will  now  be 

given.  A  pair  of  binary  sequences  of  equal  length  with  the  number  of 

like  pairs  of  one  sequence  equal  to  the  number  of  unlike  pairs  of  the 

other  sequence  for  each  possible  spacing  is  said  to  be  a  complementary 
code. 
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Two  schemes  of  representation  will  be  used  in  this  paper.  These 
schemes  will  be  used  interchangeably  at  the  convenience  of  the  author. 


The  first  representation  is  the  one  already  presented  where  each  a 
and  b^  is  either  +1  or  -le  The  operation  used  in  this  case  is 
ordinary  multiplication.  The  second  representation  will  use  the  opera¬ 
tion  of  modulo  two  addition  denoted  by  where  the  elements  a.  and 

1 

b^  are  either  0  or  1„  That  these  two  groups  (1,-1;*  )  and  (0,1;  ®  ) 
are  insomorphic  is  shown  in  the  following  tables. 


• 

1 

- 1 

© 

0 

1 

1 

1 

- 1 

0 

0 

1 

- 1 

-1 

1 

l 

.1 

0 

where  1  —  0 

% 

-i  _  i 


The  first  form  wiU  normally  be  used  when  correlating  of  the  codes 
with  their  matched  filter  is  under  discussion  while  the  second  form  will 
normally  be  used  when  discussing  the  intrinsic  properties  of  the  indi¬ 
vidual  codes. 

The  necessary  and  sufficient  condition  for  a  pair  of  sequences  to  be 

complementary,  is  in  terms  of  the  modulo  two  representation 
i=n-j  i~n- j 

F.  =  E  {a.$a.  )  =  E  (b.©b.,.®l)  for  all  j,  l<j<n-l  (2.2) 

J  i  =  1  1  13  i  =  1  1  1+J 

or  an  equivalent  formulation  by  change  of  summation  index 


f . 
J 


i=j  i=j 

£  y 

.  .  (a.0  a  (b.®  b  ...  ®  1) 

i=l  l  n-j+i  i=l  1  i  n-j+i  ’ 


(2.3) 


The  above  E  is  standard  summation.  If  modulo  two  summation  is 
desired  the  symbol  E  will  be  utilized.  (The  equivalent  conditions  are 
given  since  they  are  both  used  in  formal  proofs  in  a  later  portion  of 
the  paper.  Due  to  the  possibility  of  confusion  it  seems  more  appropriate 
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to  introduce  the  second  form  now  rather  than  later  when  needed,  ) 

has  been  used  to  represent  these  same  equations,  and  will  be 
used  elsewhere  in  the  paper  when  it  is  more  convenient  to  use  it. 

An  interesting  property  of  any  system  of  codes  is  the  possible  in- 

variances  it  may  have  under  various  types  of  transformations.  It 

therefore  seems  appropriate  to  consider  the  possible  transformations 

on  a  complementary  sequence  pair  which  leave  the  pa,ir  complementary 

and  of  the  same  length  after  the  transformation,,  A  necessary  tool  for 

this  study  is  the  property  that  a^©a  =  a^  ©  where  the  symbol 

a.  means  the  complement  of  a.,  0  =  1 1  1  =  Cb 
1  l 

Theorem  2e  1 

The  modulo  two  sum  of  a  pair  of  binary  numbers  is  equal  to 
the  modulo  two  sum  of  their  complements. 

The  truth  of  the  theorem  is  obvious  from  the  following  exhaustive 
table  of  four  possibilities: 


a.  W  a.  ,  . 
i  i+j 

1  ©  1  =  0 

1  ©  0  =  1 

0  ©  1  =  1 

0  ©  0  =  0 


Therefore  a.  ©a.  ,  .  =  a 
i  i+j  i 


a.  w  a. , . 

i  i+j 

I©1  =  0©0  =  0 

I  ©  0  =  0  ©  1  =  1 

0©  1  =  1©  0=1 

0  ©  6  =  1  «  l  =  0 


a. 


i+j 


In  each  of  the  following  proofs  it  will  be  considered  that  A  and  B 
are  a  complementary  sequence  pair  of  length  nc  Complementing  each 
element  in  a  code  is  called  complementing  the  code0 
Theorem  2*  2 

Complementing  the  A  code,  o^  the  B  code,  or  both  the 
A  and  B  codes,  results  in  a  pair  of  complementary  codes,  (The 
proof  is  given  just  for  the  A  code  but  is  identical  for  the  B  code  and 


18 


from  both  of  these  for  the  A  and  B  codes,  ) 

By  hypothe  sis 

i-n-j  i=n-j 

F.  =  E  (a. ©  a.  .)  =  E  (b  «  b  «  1). 
J  i  i+j  1  i+J 

i=  1  i=l 


Complementing  the  A  code  transforms 


a.  —  a. , 
i  i 


but  by  Theorem  2U  1 


a.  a...  =  a,  ®  a...  * 

i  i+j  i  i+j 

Upon  substitution 


i=n-j  i=n-j 

F.  =  S  (5.  «  a.  ,  .)  =  £  (b  «  b  .  .  •  1)  . 

J  i«l  ‘  1+J  iml  i  >+J 

Therefore  complementing  the  A  code  has  no  effect  on  complementarity, 
and  similarly  for  the  B  code  or  for  both  codes* 

#  >!< 

Interchanging  the  first  and  last  bits  of  a  code,  the  second  and  next 
to  last  bits,  the  third  and  third  from  last  bits,  and  so  on  is  called  time 
inversing  or  time  reversing  a  code, 

Theorem  20  3 

Time  inversing  the  A  code,  or  the  B  code,  or  both  the  A 
and  B  codes,  results  in  a  complementary  pair*  (This  proof  is  again 
just  for  the  A  code,  but  also  applies  to  the  B  code  and  therefore  to 
both  the  A  and  B  codes,  ) 

By  hypothesis 


i=n-j 

f .  =  £  (a  la.,.) 

j  i=1  i  i+j 


i~n-  j 

s  (b  «  b  ,«  1). 
i 


The  time  inverse  of  A  causes  each 
i  n+1 -l 

Applying  this  time  reversal  just  to  the  A  code  shows  that 
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i=n-j  i=n-j 

F.  =  S  (a.  ®  a...)  -  £  (a 
J  ; _ i  1  1+J 


n+l-i  ^  ?n+l-(i+j) 


i=l  ‘  i=l 

Expanding  both  expressions  for  a  few  terms  and  remembering  that 
modulo  2  addition  is  commutative  as  is  regular  addition  we  see  that 
(on  the  left)  (on  the  right) 


for 

j=n-  1 , 

Fn-1 

=  aj 

«  a. 

Fn-1 

-  a  ©  a. 

n 

n 

for 

J=n- 2 , 

Fn-2  = 

an-  1  +  a2®  an 

Fn-  2 

=  a  ©  a 
n 

These  expressions  indicate  that  by  a  change  in  summation  index  the  two 
forms  of  Fj  are  equal,  as  is  the  case,  Therefore  time  reversing  the 
A  code  has  no  effect  on  the  complementary  property,  and  similarly  for 
the  B  code  and  for  both  the  A  and  B  codesB 


The  altering  of  a  code  is  a  transformation  where  every  other  bit  of 
both  codes  is  complemented.  In  a  different  portion  of  the  paper  a  dis¬ 
tinction  is  made  between  altering  odd  bits  and  altering  even  bits.  In  the 
following  proof,  however,  just  the  even  bits  are  altered,  although  a 
slightly  modified  proof  would  hold  for  altering  the  odd  bits* 

Theorem  2e  4 

The  result  of  altering  a  complementary  pair  of  sequences  is 
again  a  complementary  sequence  pair,, 

Since  the  codes  are  a  complementary  pair, 
i-n- j  i-n-j 


F.  =  E  (a  ,  ©  a .  .  .) 

1  i  i+J 


i=  1 


2  (b.  ®  b.A.  ®  1). 

i=l 


i  i+j 


Divide  the  F.  into  two  groups,  F,  for  j=  an  even  number  and  F 
J  K  r 

for  j=  an  odd  number*  The  portion  of  the  above  equations  for  j=  an 


even  number  is 
i=n-  k 


i=n-k 


fr=  ,s,  =  s,  <bi#b,+kfi>- 

1=1  1=1 
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For  the  theorem  to  be  t  rue,  when  j  is  an  even  number  the  above  set 
of  equations  F^  must  be  satisfied  for  all  k  when  the  pair  is  altered, 
(The  proof  is  given  only  in  terms  of  A,  but  duplicate  steps  must  also 
be  applied  to  B*  ) 

When  i=an  even  number,  then  lfk-an  even  number,  In  the  alter¬ 
ing  process  since  all  even  bits  are  complemented  a.  ®  a  . ,  a.®a.  , ,  , 

1  l+k  i  l+k 

but  by  Theorem  2.  1  these  are  equal.  When  i=an  odd  number,  then 
l+k  =  an  odd  number  and  these  are  not  affected  by  the  altering  process, 
since  only  even  bits  are  complemented.  Therefore  F  is  invariant 

K 

under  altering,  or  symbolically  F^  =  F^  ^  ^ 

Examination  of  F  where  r  is  an  odd  value  of  i,  shows  that  of 

r  J 

the  pair  a^,  a^r  °ne  of  the  bits  is  complemented  and  the  other  is  un- 
chang  ed, 

a.  ®  a.  .  (a.®  a.  ,  or  a.®  a,.  )  =  a.  I  a.. 

i  l+r  -*■  i  l  +  r  i  l  +  r  i  l+r 

Therefore  F  .  =  (n-l)-F  for  both  the  A  code  and  the  B  code, 

r  alt  r 

This  signifies  the  likes  of  the  A  code  for  odd  spacings  are  changed  into 
unlikes,  and  the  unlikes  of  the  A  code  for  odd  spacings  are  changed  into 
likes;  however,  for  each  change  in  A  there  is  an  opposite  change  in 
B  so  the  total  likes  in  A  still  equals  the  total  unlikes  in  B  for  each 
spacing.  It  is  to  be  noted  that  altering  both  the  even  and  the  odd  bits  is 
the  same  as  complementing  both  codes, 

3jC  #  ♦ 

The  operation  of  interchanging  the  A  and  B  codes,  although 
trivial,  is  also  a  transformation.  There  are  a  total  of  sixty  four  pos¬ 
sible  transformations.  These  can  all  be  generated  by  combining  the 
above  listed  transformations.  They  will  be  discussed  much  more  thor¬ 
oughly  in  Chapter  3, 
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Whenever  a  new  facet  of  a  pair  of  sequences  such  as  the  comple¬ 
mentary  property  is  discovered,  it  is  natural  to  wonder  just  how  general 
this  property  might  be.  One  question  of  interest  would  be  the  possible 
length  limitations  of  the  code  pairs  that  might  be  forced  by  the  comple¬ 
mentary  property.  The  following  proof  due  to  Golay  shows  that 

Z  2 

complementary  sequences  must  have  a  length  n  of  the  form  n=x‘ +y  , 

where  n  is  an  even  integer  and  x  and  y  are  integers*  The  first 

theorem  to  be  proved  demonstrates  that  n  must  be  an  even  integer 

2  2 

and  the  second  theorem  shows  the  x  +  y  form. 

Theorem  2  .  5 

A  necessary  condition  for  a  sequence  pair  to  be  complementary 
is  that  their  length  be  an  even  number. 


By  hypothesis 
i=n-j 

F.  =  £  a. ©a.  . 

j  i=1  i  i+j 


i~n“  j 

S  b  I 
i  =  l  1 


i+j 


Be t  G.  =  F.  modulo  2,  so  that 

j  ^  i=n- j  i=n- j 

G.  =  §  a  ©  a. ,  .  =  E  b  ©b.A 

J  4-1  1  1+J  4-1  1  1  +  J 


1  . 


i=l  A  i=l 

An  expansion  of  G  for  a  few  different  values  of  j  will  suggest  the 

proof.  Start  with  j^n-1* 

G  ,  =a,#a  sbJb  t  1, 
n- 1  In  In 

G  =  a,  ©a  ..©a^Sa  =  b ,  ©b  ^ltb^ftb  91, 
n-  2  1  n  - 1  2  n  1  n  -  1  In 

G  0  =a,®a.  ^  ©a  ,  ©a  1ffia.,©a  “b.©b 

n-3  1  n-2  n-  1  n-1  3  n  1  n-2 


)  -919b-9b  ©  1 

n-1  3  n 


Change  the  order  of  the  sums  in  G  >  and  G^  ^  to  achieve  the  following 

2 


G  0=ai  ©a^®a  i®a  =b1©b0®b  ,©b  ©  §  1  , 

n-2  1  2  n-1  n  1  2  n-1  n 


G  0=a1®a^©aQ®a  9®a  ,  ©a  -b  n  ©b->©b^®b  ->®b  ,  ©b  8 

n-3  1  2  3  n-2  n-1  n  1  2  3  n-2  n-1  n 


In  general  for  j-n-r 

G  =a,  ffia^©u  .  „  ©a  ©a  t  ,nla  ,  ,7A 

n-r  1  2  r  n-(r+l)  n-(r+Z) 


©a  , ©a  =b, ©b0( 
n-1  n  1  2 


c  •  o  c 
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>  ®b  ,  .  ,,®b  .  ®b  ,®b  «  £  1  , 

r  n-(4+l  n-( r+2)  n- 1  n  \  ’ 


and  for  j=n-(r  +  l) 


Gn-(r  +  l)"al9a2®  •-fer,ar  +  ]»an-r®an-(r*n®---9a„-l®a„=bl»b2*--- 


.  ,»b  ffib  .  ,»b  »b  ,  ,  ,«b  ,  ,,.l 

r  r+1  n-r  n-(r+z)  n-(r+2) 


r+1 

?  '• 


Forming  the  sum  G  ©  G  .  ,,  will  leave  only  five  terms  in  the 

n-r  n-(r+l)  7 

equation  since  G  and  G  t  ,  ,  x  differ  only  in  a  ,  a  .  ,  b  ,  n  , 

n-r  n-  (r  +  1)  y  n-r  r  +  1'  r+1' 

and  a,  bcause 
r  r+1 

Si®  S  1  =  1. 

1  1 

G  ©G  /  ,  v  =  a  , ,  ©a  =  b  .  .  ©b  0  1. 
n-r  n-(r+a)  r+1  n-r  r+1  n-r 


Adding  b^+i©b  modulo  2  to  both  sides  of  the  equation  gives 

a  ,  ©a  ©b  .  t  ®b  =*1.  (2.4) 

r+1  n-r  r+1  n-r  x  ' 

Let  n=2s-l  and  r=s-l;  then  to  satisfy  the  equation  just  derived 


a  ©a->  i  .  i  ©b  ©b ^  t  =  a  ©a  ®b  ®b  =  1  ,  which 

s  Zs-l-s+1  s  Zs-l-s+1  s  s  s  s 

is  obviously  false. 

Therefore  n  cannot  have  the  form  2s-  1  and  must  be  an  even  number 


* 


* 


Equation  2.  4  which  was  a  step  in  the  previous  proof,  is  a  very 
important  necessary  condition  for  a  pair  of  codes  to  be  complementary 
and  will  be  referred  to  in  the  future  as  the  parity  check.  The  next 
theorem  which  shows  the  n=x  +y  form,  starts  with  the  basic  assump¬ 
tion  that  A  and  B  are  a  complementary  pair  of  length  n. 

Theorem  2,  6 

A  necessary  condition  for  a  pair  of  codes  to  be  complementary 
is  that  their  length  be  the  sum  of  the  squares  of  two  integers. 

By  hypothesis 


i-n-j 


i  =  n- 


F.  =  S  (a  ©a  ^.01)  = 

j  i=i  1  i+j 

Assume  that  the  A  code  satisfies 


E 

i=  1 

F. 

J 


<bi®bitj>  • 

and  that  the  A  sequence  contains 
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p  ones  and  n- p  zeros.  The  A  code  has  a  weight  oi  p.  Similarly 

assume  that  the  B  code  satisfies  F\  and  that  its  sequence  contains 

q  ones  and  n-q  zeros.  The  B  code  has  a  weight  of  q,  For  the  total 

likes  in  the  A  code,  each  bit  is  matched  against  each  other  bit  exactly 

once.  For  the  total  unlikes  in  the  B  code  each  bit  is  again  matched 
against  each  other  bit  once.  There  are  two  possibilities  of  likes  in  the 
A  code,  1  matched  against  1,  and  0  matched  against  0*  The  total 
number  of  like  pairs  for  the  A  code  is  therefore  the  number  of  com¬ 
binations  of  p  ones  taken  2  at  a  time  plus  the  number  of  combinations 
of  n-p  zeros  taken  2  at  a  time,  for  the  unlike  pairs  in  the  B  code, 
q  ones  are  matched  against  n-q  zeros. 

This  gives  the  equation 

=  q(n-q) 

expanding  and  simplifying 
n=n^-  2qn+2q"'  +  2p  -  Znp, 
adding  and  subtracting  Zpq  gives 

n=n^- 2qn+2q^+2p  - Znp+Zqp- Zqp, 

collecting  terms  gives 

2  2  2  2  2 

n=  n  +q  +p  +2qp-2np-2nq  +  p  -2pq  +q 

combining  terms  gives 

n=  (n-p-q)2  +  (p-q)2  5) 

Therefore  complementary  sequence  lengths  are  permissible  only  in 
lengths  which  are  even  numbers  and  formed  by  the  sum  of  two  squared 
integers. 

*  *  * 

A  list  of  all  possible  code  lengths  up  to  200  with  the  number  of 
ones  allowed  in  the  A  and  B  codes  is  given  in  Table  2.  lv 
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^Indicates  possible  kernel 


Code 

Uno  rde  red 

Code 

Unordered 

.Lengths 

W  eights  of(A,  B) 

L»e  ngths 

W eights  of  (A,  B) 

*2 

(Z,  1 )( 1 , 0) 

104 

(5  8,  561(58, 481(48, 46)( 56,  46) 

4 

(3 , 3)(  1 ,  1)( 3 ,  1) 

*106 

(60,  55) (5 1 , 46)(60,  5 1 )(  55, 46) 

8 

(6 , 4)(4,  2) 

1 16 

(65,  61  )(55 , 51  )(65,  55)(  61,51) 

*10 

( 7 , 6 )(4,  3)(  7,  4)(6,  3) 

*122 

(67,  66)(56, 55)(6 7,  56)(66,  55) 

16 

(10,  10)(6,  6 )( 1 0 ,  6) 

128 

(72,  64)(64,  56) 

*18 

(12,  9)(9 , 6) 

*130 

(72,  69) (72, 611(61, 58)(64,  5  7) 

20 

(13,  1 1 ) ( 9 ,  7){  11,  7) 

(73,  66)( 73 , 64)(64,  57)(66,  57) 

*26 

(16,  1  5)  ( 1 1  ,  10) 

1  36 

(76,  701(66, 60)( 76,  66)(70,  60) 

(15,  10)(  16  ,  11) 

144 

(78,  78)(66, 661(78,  66) 

32 

(20,  1 6)  ( 1 6 ,  12) 

*  1  46 

(81,  76)(70, 65)(81 , 70)(76,  65) 

*34 

(21,  1 8)  (16,  13) 

148 

(81,  791(69, 67)(81, 691(79,67) 

(21,  1  6)  ( 1 8 ,  13) 

160 

(88,  84) (76 , 72) (8 8,  76)(84,  72) 

36 

(21,  21)(1  5,  1  5)(2  1  ,  1  5) 

*  1  62 

(90,  81)(81 , 72) 

40 

(24,  22)(18,  16) 

(24,  181(22,  16) 

1  64 

(91, 83)(81, 73)(9  1,81  )(83,  73) 

*50 

(30,  251(25, 20)(29, 28) 
(29,  22)(28, 21 ) ( 2 2 , 21 ) 

*1  70 

(94,  87)(83 ,  761(94, 83)(87,  76) 
(92,  911(79,  78)(92,  79)(91, 78) 

52 

(31, 2 7) (31 , 25) 

(25,  2 1 ) ( 2 7 , 2 1 ) 

*  1  78 

(97, 94)(97, 84)(84, 81)(94, 81) 

*58 

(34,  3 1 ) (2 7 ,  24) 

(34,  27)  (31 , 24) 

1  80 

(99, 93) (87, 81)(99,87)(93,  81) 

64 

(36,  28)  (36 , 36)  (2  8, 28) 

*194 

(106,  101)(93s  88) 

68 

(39,  37)(31 , 29) 

(106,93)(101,88) 

(39,  31)(37 ,  29) 

196 

(105, 105)(91, 91)(105, 91) 

72 

(42,  361(36, 30) 

200 

(110,  100)(108,  106)(  1 00 , 90) 

*74 

(43,  381(36, 31) 
(43,36)(38,31) 

(94,  92)(94,  108)(  106,  92) 

80 

(46,  421(38, 34) 

All  pos  sible 

complements  ry 

*82 

*90 

*98 

100 


(46,  38){42,34) 

(46,451(37,36) 

(46,  3 7) (45 , 36) 

(51, 48)  ((51,42) 

(42,  39)(48 , 39) 

(56,  49)(49 , 42) 

(55,  55) (55 , 45)(45, 45) 
(57,  51 )  (57 , 49)( 5  1 , 43) 
(49,43) 


code  lengths  up  to  200  with 
unordered  possible  weights. 


TABLE  2.1 
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The  shortest  possible  complementary  pair  is  A=ll,  B-iO.  This  pair  or 
any  of  its  transformations  is  called  a  kernel  of  length  two,  or  the  quad. 

A  kernel  is  a  basic  length  code  which  cannot  be  decomposed  into  shorter 
length  codes  by  an  inversion  of  the  standard  generating  methods  to  be 
explained  later.  Some  possible  kernel  lengths  are  2,  10,  1  8,  26,  34, 

50,  etc.  although  codes  for  all  of  these  do  not  exist.  Complementary 
sequences  which  are  not  kernels  are  called  composite  complementary 
sequence  s . 

The  above  list  of  kernel  lengths  did  not  include  n=“4  or  n  =  8 
which  are  possible  complementary  sequence  lengths*  These  if  they 
exist  must  therefore  be  composite.  All  possible  codes  of  length  four 
can  be  generated  from  the  exhaustive  list  of  possibilities  for  four  binary 
digits  given  in  Table  2e  2. 


0000 

1  000 

0001 

1001 

0010 

1010 

0011 

101  1 

0100 

1  100 

0101 

1  101 

0110 

1  110 

0111 

mi 

Table  2e  2W  All  Possible  Binary  Numbers  of  Length  Four. 

The  unordered  possible  pairs  of  ones  in  the  codes  are  (1,1),  (3,3) 
and  (3,  1).  Thus  all  possible  codes  of  length  four  have  either  one  1  or 
three  1 1  s  because  of  this  limitation*  This  reduces  the  table  to  the 


eight  following  numbers  of  length  four: 


0001 

0111 

0010 

1  101 

0100 

101  1 

1000 

1110 

The  possible  unordered  (1,  1) 

code 

pairs  are  the  following  four: 

1000 

1000 

0001 

0001 

0100 

0010 

0100 

001  0 
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The  possible  unordered  (3,  1)  pairs  of  codes  are  the  following  eight: 

0111  0111  1110  1110  1011  1011  1101  1101 

0010  0100  0100  0010  0001  1000  0001  1000 

The  possible  unordered  (3,3)  code  pairs  are 

0111  0111  1110  1110 

1011  1101  101.1  1101 


Comparing  the  last  code  to  the  kernel  of  length  two,  A=1  1  ,  B=10, 
shows  that  it  might  have  been  constructed  by  writing  in  time  sequence 

S  i  =  AB=a  ^a^b  i  b-,  =  1110, 

=  AB=a^a^>b,b-)  =  1101. 

The  next  to  last  sequence  pair  might  have  been  formed  by  inter¬ 
lacing  A  and  B  in  the  following  manner: 

Ti  =  albla2b2  =  1U0* 

Tg  =  >  =  1011, 

All  the  rest  of  the  16  pairs  could  be  considered  either  as  transforma¬ 
tions  of  the  last  pair,  or  as  being  formed  by  the  same  two  composite 
generating  operations  used  above  on  the  8  transformations  of  the  ker¬ 
nel  of  length  two,  A  general  proof  will  now  be  given  to  show  that  both 
the  time  sequence  scheme  (Si  ,  S2)  and  the  interlace  scheme  (Ti,  T2) 
will  always  form  complementary  sequence  pairs  providing  that  A  and 
B  are  a  complementary  sequence  pair. 

Theorem  20  7 


If  A— a^a^9-^. 


a  -1  a 
n-  1  n 


B_blb2b3’  ‘  *  “bn-  lbn 

are  a  complementary  sequence  pair,  then 

..  .  .  b 


C  =  a.  a, ...  a  b 

12  n  1  2 


n 


D  =  a,  *  a  b,bv 

12  3  n  1  2 

are  a  complementary  sequence  pair. 


n 
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Theorem  2,  8 

If  (A,  B)  are  a  complementary  sequence  pair,  then 

C  -  a  i  b  i  3.  -)e4,,UBo»3.  b 

1  1  < t  c  n  n 

D  =  a1b1a2b2.......anbn 

form  a  complementary  sequence  pair0 

The  following  notation  will  be  used  in  the  proof: 
i=n- j 

U  =  E  a  fa.  .  for  all  j  1  <  j  <  n-i 

L.  =  1~n~ Ja.$a.,  .  f  1  for  all  j  a  <  j  <  n-1  . 

a  Ei  i+j  J  J 

i=  1 

It  is  understood  that  the  two  CD  sequences  in  the  two  theorems 
above  are  not  the  same  but  because  the  proofs  are  identical  two  sets  of 
symbols  will  not  be  used.  For.each  spacing  j  theie  are  three  possible 
ways  for  bits  to  match,  within  A,  within  B ,  or  from  a  bit  of  A  to  a 
bit  in  B*  These  possibilities  are  denoted  by  I  and  Labor  Ua’ 

U^,  U  ^  re spectively* 

The  necessary  and  sufficient  condition  for  C  and  D  to  be  a 

complementary  sequence  pair  is  given  by  L  =  La+  L^+L^  and 

U  =  U  +Ur  +  U  -  where  L  = 
d  a  b  ab  c  d 

Now  L  =  U,  and  L,  =  U  ,  ,  ...  r  , 

a  b  jq  a  by  deiinition  ot  complementary. 


Adding 

L.  +  L,  =  U  +UU, 
a  d  a.  b 

b  ut  t  -l 

i=n-l  l = n  -  1_ 

U,  =  Ur  since  E  b.fb  .  .  =  E  b  ®  ^i+i 

b  b  ,i  l+l  ,  1  J 

i=l  J  i“l 

so  substituting 

L  +  D  =  U  +  Ur 
a  o  a  b 

All  bit  matches  from  A  to  B  in  the  C  sequence  are  paired  in  the  D 

sequence  with  bit  matches  from  A  to  B  therefore. 
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and  adding 


L  + 
a 


L  ,  =  U 
ab  a 


+  uB  +  u 


ab 


or  L  =  U  , 
c  d 


* 


* 


* 


Golay  describes  two  other  methods  of  generating  composite  sequen¬ 
ces  from  shorter  complementary  sequences,.  These  methods  will  not  be 
used  too  extensively  in  the  rest  of  the  paper  and  will  therefore  just  be 
mentioned  rather  than  proved*  Given  two  complementary  sequences 
pairs  (A,  B)  (C,D),  A  of  length  m  and  C  of  length  n. 

Let 


Ux  =  ACl  AC? . .  A°n  Bdl  Bd? .........  Bdn 

U2  =  Adn.  ...........  Adl  ’BCn . .  .  .  BCl , 


(2.6) 


and 

Vi  =  Aci  Bdl . . . .  Acn  Bdn 

,  ,  _  (2.  7) 

V3  =  Adn  BCn . .  . . . . . . .  Adl  BC\ 

It  can  be  shown  that  the  pairs  of  U  are  complementary  as  are  the  pairs 
of  V.  Where,  if  an  exponent  is  a  one  the  A  or  B  code  is  left  un¬ 
changed,  and  if  the  exponent  is  a  zero  the  A  or  B  code  is  comple¬ 
mented.  The  lengths  of  codes  which  can  be  generated  from  these  two 
methods  are  Zmn  where  m  is  the  length  of  the  (A,B)  sequence  pair 
and  n  is  the  length  of  the  (C  ,  D)  sequence  pair*  There  is  also  another 
special  method  of  generation  which  applies  only  to  codes  of  length  2ro 
This  method  is  explained  in  Chapter  5  where  it  is  necessary  for  the 
completeness  of  a  proof. 

The  kernel  of  length  2,  or  the  quad,  has  been  discussed  in  the 
previous  pages.  There  are  two  kernels  of  length  10.  These  are 
1001010001  and  10000001  10,  and  a  second  kernel  0101000011  and 
00001  001  10.  The  next  possible  kernel  size  is  18  and  Golay  has  proved 


by  exhaustive  search  that  no  kernel  of  this  length  exists,  Kruskal  has 

1 3 

since  then  completed  the  proof  analytically*  One  of  the  goals  of  the 
writer  in  his  study  was  to  make  an  exhaustive  search  for  codes  of  kernel 
length  26.  This  search  disclosed  there  was  only  one  kernel  of  length 
26.  Chapter  7  of  this  paper  describes  this  search  in  detail  and 
Chapter  8  describes  attempts  to  find  kernels  of  length  34. 
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CHAPTER  III 

THE  GROUP  OF  OPERATIONS 

A  large  portion  of  this  investigation  is  devoted  to  the  set  of  allow¬ 
able  transformations  on  complementary  sequences.  These  allowable 
transformations  possess  the  characteristic  that  the  length  of  the  sequen¬ 
ces  and  the  complementary  property  are  invariant  under  transformation. 
To  be  useful  in  searching  for  new  kernels  it  is  also  desirable  that  these 
transformations  have  the  following  features: 

1.  That  they  form  a  group  for  all  general  transformations. 
Where  a  general  transformation  is  one  which  can  be 
applied  to  any  complementary  code  pair. 

2.  The  transformations  as  represented  be  simple  manipula¬ 
tions  on  the  elements  of  the  group. 

3.  The  supplementary  characteristic  (to  be  explained  in 
Chapter  6)  be  a  property  of  the  representation. 

a 

4.  The  parity  check  be  easily  made. 

Some  of  these  transformation  operations  were  discussed,  and  the 
proof  of  invariance  given  for  all  group  generators,  in  the  background 
chapter,  Chapter  2,  on  complementary  sequences.  As  a  continuation, 
a  brief  review  of  the  transformation  operations  and  a  tagging  of  symbols 
to  these  operations  will  prove  useful.  Each  basic  operation  symbol  and 
its  definition  is  given  in  Table  3.  10  More  definitions  are  given  than  are 
necessary  to  generate  the  group;  however,  the  redundance  seems  to  the 
writer  to  be  an  aid  to  understanding. 

A  check  through  Table  3.  1  shows  that  the  transformations  affect 

the  complementary  sequence  pair  in  four  different  ways.  Ci  ,  C2  , 

Tx  ,  T2  affect  all  the  bits  of  one  of  the  code  pair.  EST,C  affect  all  bits 

of  both  codes,  while  Ai  and  A2  affect  half  of  the  bits  of  both  codes.  The 

identity,  I,  of  course  has  no  affect  on  either  code.  A  reasonable  way 

to  form  the  elements  of  the  group  to  account  for  the  above  affects  would 

be  to  divide  each  code  of  the  pair  into  two  units.  These  would  be  A 

(even  bits)  and  A^  (odd  bits),  also  (even  bits)  and  Bq  (odd  bits). 
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I 


The  Identity 


a.  -*  a. 

1  1 

b  -»  b 
1  i 


Ti 


A  code  time  reversal 


T3 


a.  “* a  . 

i  n+l-i 

b.  -b. 
i  i 

B  code  time  reversal 


a. 

i 

b. 

i 


a. 

i 


n+ 1  -  i 

Time  reversal  of  both  codes 


a.  -*■  a  ,  ,  . 

i  n+ 1  -  i 


Ai 


b.  -  b  X1  ; 
i  n+ 1  - 1 


Alte  r  odd  bits 


a2i-  1  "*  a2i-  1 

a,.  -*  a-,. 

2i  2i 

b2i-  1  "  b2i-  1 

b2i  -  b2i 

A3 

Alte  r  even  bits 

a2i-  l-*  a2i-  1 

a2i~*  a2i 

b2i-  1  b2i-  1 

b-,  b? 

2i  2i 

Ci 

Complement  the  A  code 

a.  -*  a . 
i  i 

b  -  b 

i  i 

c3 

Complement  the  B  code 

a.  -*  a. 

i  i 

b.  —  b 

i  a 

c 

Complement  both  codes 

a .  -*  a . 

i  i 

b.  -  b. 

i  i 

E 

Exchange  the  A  code  with 

a.  -»  b. 

the  B  co 

,  i  i 

b  -*  a. 

i  i 
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Some 


operations. 


TABLE  30  1 


Ao  ala3a5 . *  ”  '  an-  1 

Ae  =  a2a4a6‘  an 

B  =  b,  b~bK.  e  .  * . b  , 

o  13  5  n- 1 

B  =  b->b Ab, . .  .  .  .  b 

e  Z  4  6  n 

However,  this  particular  division  for  representing  the  complemen¬ 
tary  sequences  proves  to  be  very  awkward  for  the  time  reversal  opera¬ 
tions  (T}  ,T2  >  T)  and  for  the  ease  of  the  parity  check  condition. 

Fortunately,  one  simple  change,  that  of  time  reversing  the  second  and 

14 

fourth  row  (even  bits  of  A  and  B),  clears  up  these  difficulties. 

The  method  for  converting  from  a  complementary  sequence  pair 

into  the  standard  group  form  is  therefore 

<,.  I=a,  a.a,. •  . .  e  a  , 

135  n- 1 

TT 

II=anan-2--—a2  (3-D 

,111=  b  .  b,b  r.  ...  b  . 

1-  U  /  135  n- 1 

B=b  b->b, . b  .b  - ( 

12  3  n-  1  n  \  T»r  ,  ,  , 

MV=b  b  . b-> 

n  n-  2  2 

If  the  matrix  (I  II  III  IV)  is  taken  as  the  symbol  for  the  above  grouping 
the  Identity  operation,  I,  should  yield  (I  II  III  IV).  Post  multiplying 
(I  II  III  IV) 

10  0  0 

by  0  0^0  8^ves  (I  II  HI  IV),  therefore  I  can  be  taken 

0  0  0  1 


as  the  identity  matrix  of  rank  4,  using  the  operation  matrix  multipli¬ 
cation  on  the  right.  Similarly  Tn  has  the  effect  of  exchanging  the  first 
and  second  columns  yielding  (II  I  III  IV).  Now  post  multiplying 


(I  II  III  IV)  by 

0  10  0. 

I  1  0  yields  (II  I  III  IV)  the  desired  result  so  the 

0  0  0  1 

matrix  just  used  could  represent  Ti  .  Similarly, 
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1000  0100  1000  1000 

T  0  1  0  0  T  1  00  0  r  _  0  1  o  or_o  1  00 

2  0  0  0  1  000  1  ^100l0v'a0010 

0010  0010  0001  oool 


Ai= 


1  0  0  0  1  0  0  0 

0  100  01  0  0 

ooioA8  ooio 

0  0  0  1  0  0  0  1 


0  0  10 
0  0  0  1 
10  0  0 
0  10  0 


Matricies  are  often  a  very  convenient  method  of  expressing  trans¬ 
formations.  However,  when  the  elements  are  few  and  the  numbers  in 
them  simple  it  is  often  more  convenient  to  express  the  operation  dif¬ 
ferently.  For  instance  (I  II  III  IV)  operated  upon  by  Ti  yields 

(II  I  III  IV).  This  could  be  just  as  well  expressed  as 

_  I  II  III  IV 
11  II  I  III  IV  • 


This  indicates  that  the  top  row  i§  transformed  into  the  bottom  row  by  the 
operation,  or  in  even  more  shortened  form,  Ti  =  (II  I  III  IV)  with  the 
top  row  understood  which  means,  as  is  common  practice,  Ti  (I  II  III  IV) 
=  (II  I  III  IV).  Another  example  is  =  (I  II  III  IV),  [which  means 

Ai  (I  II  III  III)  =  ( I  II  III  IV)]  where  the  bar  across  I  and  III  indicates 
the  negative  element  in  the  matrix  or  in  actuallity  the  complement  of  the 
portion  of  the  code  contained  in  I  and  III.  Table  3,  2  is  a  list  of  all 
32  transforms  generated  from  Table  3.  1  both  by  symbol  and  group 
element  representation.  There  are  actually  64  transform  operations 
in  the  entire  group.  However,  the  exchange  operation  (III  IV  I  II)  is 
not  considered  along  with  all  the  32  elements  it  would  generate.  In 
this  paper  the  pair  (I  II)  and  (III  IV)  are  generally  considered  as 
being  an  unordered  pair,  which  eliminates  the  exchange  operation. 

The  second  feature  listed  under  desirable  characteristics,  that 

of  simple  manipulation  of  the  elements  of  the  group  is  satisfied. 

The  supplementary  property  which  is  discussed  in  Chapter  6  is 

also  satisfied  by  the  formulation  of  groups  from  complementary 
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sequences  by  this  method.  The  parity  check,  a  ®b  $a  ,  ,  <&b  .  ,  .  =  1, 

1  11  n+l-i  n+l-i 

is  also  satisfied  under  all  transforms  and  is  easily  checked  since  each 
column  of  equation  3.  1  is  one  solution  of  the  parity  check. 


I  =  a 

1  a3' 

a  c  a 

i 

II  =  a 

5 

i  *  *  #  n- 

a 

n  n 

-2  *  • 

• • *an+l-i  a2 

III  =  b 

lb3 

b  5  ■  • 

•  •  •  o  •  •  • 

i 

IV  =  b 

b 

n  n 

-2  *  * 

*  °  ’  bn+i-lb2 

I  =  ( I  II 

III 

IV) 

R 

= 

(i  IT 

IV 

ffi) 

A1=  (T  II 

Tn 

IV) 

Q 

= 

(ii  T 

IV 

HT) 

a2=  (I  7T 

III 

m 

P 

(IT  I 

Tv 

hi) 

T i=  (II  I 

III 

IV) 

O 

sr 

or  r 

m 

IV) 

T3=  (I  II 

IV 

III) 

N 

= 

(H  I 

m 

Tv) 

T  =  (II  I 

IV 

III) 

M 

= 

(T  n 

IV 

in) 

Ci=  (T  il 

III 

IV) 

L 

=s 

(ITT 

m 

Tv) 

c3=  (I  II 

m 

TV) 

K 

= 

(1  Ti 

Tv 

IH) 

c  =  (1  n 

TIT 

Tv) 

J 

= 

(IT  T 

Tv 

TTI) 

z  =  (ii  T 

Hi 

IV) 

H 

= 

(i  ii 

Tv 

HIT 

Y  =  (IT  I 

HI 

IV) 

G 

= 

(TTI 

IV 

III) 

x  =  (IT  i 

hi 

IV) 

F 

= 

(I  ~n 

m 

IV) 

w  =  (ii  T 

hi 

“IV) 

D 

= 

(J  ii 

m 

TV) 

V  =  (1  II 

IV 

HI) 

ct> 

(ii  i 

Tv 

TTI) 

U  =  (I  II 

Fv 

III) 

B 

= 

or  i 

IV 

Tn) 

s  =  (i  IT 

IT 

III) 

IT 

= 

(ii  i  ■ 

Tv 

iu) 

E  =  (III  IV  I  II) 


Elements  of  the  unordered  operations  group. 


TABLE  3.  Z 
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pq 

C  PH  3CQ  >  OSHffiWOJt5HUOm>-i><OilHp.'isI<;EHbHoW 

maPn  t'SuittJHP‘KW2H2;HO  juo!>MtsiD>flX<bE-iUH4 

joh  WhSW^bXsiMiHOiW^PMOHbE-ia 
h<4POOtHD5CQOH2;hJ2WffiHe;'Q.^X!>5:p)  cmaHU^bF4 
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J!30H>H04WR<fe<tafflPHDDiU  U'S.h  Hec5SXNm>2 
ZHWKMk-CJOOPh  COifQ^h^filxJtJH^H'Q.OgUEHDCtiSH  J 
55  hJ  H  O  x  03  O'®.  S3  -Sp4^fl|CqoifM>UHTOto2HM^lHlDDiH 
OEHjlz;Ki5<!oag<ft<fiCst|iicq>mHUh  Jtsi.Ii;  b  H  ^  (S  t)  W 

CM  CM  rH  ■ — I  rH  CM 

Ph  tfflaHbWl^ilE-iSMffiE-iZhJOHOtDlS^^^^Mfi  [5'Q.O  U  ft 

<-<03  r-H  0}  C\)  rH 

Oi PQ  tP4bE-|>isi<lMffiE-'2;(JOE-i^UHDiD!>Hmi><<<>h'S.c5ot-)fl 
Wt)WH(fl&SH4'Q.E-u5  0OHO  JHOi|lc4  tfi^^Siag^OlxJ 

cot3Dil>SiilPH4KcbH^s.bouoHo^  c  ca  oa-<  c>3fe  >  h  W  h  j  m 

rH  03  H  03  03  H 

>(UDMCSa<!M'S.b!5H0HUU20ffl 

t£iOOHUbt5  E-i^S-W  H^<MiHCQCria«20!I12!> 

CM  CM  rH  rH  rH  CM 

X  S=  H^^oclPHmuuuH'siHb^KSRh  t  <5  a>  ffl  5  J  h  M  ^ 
HNX>Htj^mtfHOUOH'Q.bi5HM'?4ao  fDPLi!>OS2ffim 

i — i  < — <  03  03  rH  03 

siH>XO!?<ioi5»UHUutb  b'Q.H  aw&cfifQ^PimptfHhJHWD 

03  rH  03  rH  rH  03 

Hej'Q.bP4atJHHm5D|>tS  Jo|z;u  cmfiiiU 

CO  H  rl  03  rH  03 

H21l!Wt3tfHHEHta(CQOIfi<lft<l>|H  O'Q.O  bOJUOWHX^a 

i — (  03  rH  03  rH  03 

HOghJ|H>HHE-ipC(<!fi<IC9(ic  taoSDOOcSObW^aScqX^MBH 
^fih^OH^mfMJSHOMSHffib^^iHlDlsI^OJDi  tOO^eJCQ 
-^Fnp<?HUiSl|>0)^0(JghffiW2hh3i>c>tf  XmfilDmuuej'Q.  t 

03  rH  03  rH  rH  0i 

l>q>P40)^0SHEHe2'S.Rft  t  PQ  O 

03  rH  rH  03  03  rH 

U  U  H  O  p  ^  M  N1  P  >  Di  W 

rH  03  rH  03  rH  03 

O  i— i  00^R02c5^N>X(tiMt3>ffl  CfHjE-|z;!il1SMH<tl'i!OJ(ii4 

HUOO^-liH^frcsllHixlls^pWK  0)Ph  O^^cq^bMtt&HPfq  insw 


O  PQ  Eh 
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TABLE  3.3 

Operations  Multiplication  Table 
(Side)  (Top)  =  Operation 


Another  feature  inherent  in  this  formulation  of  the  group,  which  is 
very  useful  in  searching  for  new  kernels,  is  the  interlace  breakdown 
into  the  group.  This  is  useful  because  a  kernel  may  be  considered  as  the 
interlace  of  two  half  length  codes  where,  as  will  be  shown  in  Chapter  7, 
the  number  of  zeros  and  ones  in  each  half  length  code  are  easily  pre- 
dete  rmined. 

Although  this  group  of  operations  is  indeed  a  group,  it  has  never 

been  proved  formally.  In  order  for  a  collection  of  elements  to  be  a 

1  5 

group,  they  must  satisfy  the  following  four  conditions: 

1.  A  group  has  closure. 

2.  A  group  is  associative  under  the  group  operation. 

3.  A  group  has  an  identity  element. 

4.  Each  element  in  a  group  has  an  inverse. 

That  this  operation  group  has  closure  is  demonstrated  by  Table  3.  3 
which  is  the  full  multiplication  table  of  this  group  of  order  32.  The 
elements  under  the  operation  multiplication  are  associative  since  they 
are  expressible  as  matricies  and  matrix  multiplication  is  associative. 
The  element  I  is  the  identity,  and  as  seen  from  Table  3.  3  all  elements 
have  an  inverse.  Therefore  these  32  elements  have  satisfied  all  the 
requirements  for  a  group. 

There  are  three  special  subgroups  of  order  4  which  within  the 
subgroup  deal  with  only  one  type  of  operation,  Ti  T2  T I,  Ai  A2  Cl, 

Ci  C2  Cl.  All  of  these  subgroups  have  the  unusual  property  that  each 
element  in  the  subgroup  is  its  own  inverse.  This  group  of  Order  32 
has  been  identified  as  isomorphic  with  Senior's  group  number  44.  ^ 
This  will  be  discussed  in  some  detail  at  the  end  of  this  chapter. 

A  study  of  complementary  codes  in  (I  II  III  IV)  form,  or  as  we 

will  call  it  sequence  quadruple  form,  reveals  many  interesting  proper¬ 
ties.  Some  of  these  properties  and  their  formal  proofs  will  follow. 
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In  these  proofs  the  complementary  pair  (A  B)  will  be  assumed  to  be  of 

i=j 

length  n.  The  symbol  L.  will  mean  E  (a  #a  .  .  I  1)  for  all 

.  ,  1  m*  i  —  i 

. i=l  J 

j,  1  <  j  <  n-  1  and  similarly  Ug  =  ^  ^ 


i=l 


b.«b  ,  . 

i  n+i-j 


Theorem  3.  1 

If  (A,B)  are  a  complementary  sequence  pair  and  if  either  (I,  II) 
or  (III,  IV)  in  sequence  quadruple  form  are  complementary,  then  the 
other  pair  must  also  be  complementary. 

1.  Assume  (I,  II)  are  comlementary. 


2. 

3. 

4. 

5. 


6. 

7. 

8. 

9. 

1  0. 

1  1. 
1  2. 

1  3. 


Lj  =  UjT  by  the  definition  of  complementarity. 

UT  =  Ljt  by  the  definition  of  complementarity 
Adding  the  equations  in  steps  2  and  3  gives 
LI  +  LII  "  UI  +  UII* 

Since  only  even  spacing  s  are  concerned  in  the  sequence 
quadruple  form,  when  related  back  to  the  original  (A,B) 
pair,  the  likes  of  A  are  equal  to  the  unlikes  of  A  for  all 
even  spacings, 
even  space. 

and  by  the  definition  of  complementarity,, 

Therefore  L  =U  for  even  spacings. 

±j  Jd 


This  will  be  expressed  as  for 


Lm+U m=k.  since  the  total  number  of  matches  at  any  spacing 
is  equal  to  the  likes  plus  the  unlikes,  and  is  dependent  upon  j. 


LIV+UIV  "kj* 

Adding  8  and  9  gives  By  1  jy  ”  By 

(uin+uiv)  55  2k; 


J 


L'IN+Liy  -  k^  ,  by  step  A 

Lj.^r-Ujyj=  0  by  subtracting  the  equation  in  step  8  from  that 
in  step  1  10 

Therefore  B  v  =U^j  and  the  (III,  IV)  pair  is  complementary. 
Since  the  proof  would  have  been  the  same  if  (IIIS  IV)  were 
chosen  complementary  rather  that  (I,  II),  the  general  state¬ 
ment  is  true* 
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Theore m  3  „  2 


If  (A,B)  are  a  complementary  sequence  pair  in  sequence  quadruple 
form  and  if  one  of  the  pair  (I,  II)  or  (III,  IV)  is  not  complementary 
neither  is  the  other. 

1.  Assume  (I,  II)  not  complementary. 

2.  There  are  two  possibilities  for  (III,  IV)  either  complementary 
or  not  complementary^ 

3.  Assume  (III,  IV)  are  complementary,  then  by  the  Theorem  3.1 
the  (I,  II)  pair  must  also  be  complementary,  but  this  is  con¬ 
trary  to  assumption, 

4.  Therefore  (III,  IV)  cannot  be  a  complementary  pair. 


* 


* 


Theorem  3.  3 

A  necessary  and  sufficient  condition  for  (I,  II)  and  (III,  IV)  to  be 

) 

two  complementary  pairs  when  a  complementary  sequence  pair  is 
written  in  the  sequence  quadruple  form  is  if  the  likes  for  even  spacing  of 
one  of  the  code  pair  is  equal  to  the  unlikes  at  the  same  spacing  for  the 
same  code  of  the  pair.  First  the  sufficient  portion: 

1.  L.  =  U  for  even  spacing  s. 

A  A 

2.  L, i  +  L«T7  =  L  for  even  s pacings* 

3.  Uj  +  Up  =  UA  =i  La  for  even  spacmgs, 

4.  L.  +  L.tj  +  Uj  +  Up  =  2L^  for  even  spacings  by  adding  the 
equation  in  step  2  and  step  3. 

5.  L.t  +  UT  =  k.  since  the  total  of  likes  and  unlikes  must  equal  the 

i  J 


pos  sibilities 


6. 

7. 

8. 
9. 


LII  +  UII  =  kj  ' 

Lt  +  UT  +  L.tt  +  UTT  =  2k.  =  2L.  by  adding  the  equations  in  steps 

I  I  II  ii  J  A 

5  and  6  and  comparing  to  those  in  step  4„ 

L,t  +  L,tt  =  k.  rewriting  step  2  with  k.  substituted  for  L,  . 

I  II  J  J  A 

L  -  Uj  =  0  subtracting  the  equation  in  step  5  from  that  in 

s  te  p  8 , 
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10.  Therefore  Lu  =  Uj  and  (I,  II)  are  a  complementary  pair  and  by 
Theorem  3.  1  (III,  IV)  are  also  a  complementary  pair. 


11. 

12. 

13. 


For  the  proof  of  the  necessary  condition  assume  that  (I,  II)  are 
complementary  and  that  U  <  L ^  for  some  even  spacing. 

Therefore  Uj.  +  <  L  +  L  for  some  even  spacing. 


Howeve  r 

=  U, 


I 


II 


,  by  assumption  (I,  II) 
and  =  L)T 


are  complementary  and 


14.  Adding  the  two  equations  in  step  13  gives  T(  +  L^  (  =  Uj+Up  . 

15.  This  is  incompatible  with  step  12,  therefore  it  is  necessary  that 

UA=  LA  or  Ug  =  for  all  even  spacing  s  in  order  to  have 

(I,  II)  and  (III,  IV)  complementary. 


* 


Theorem  3.4 

*> 

If  the  sequence  quadruple  form  of  a  complementary  pair  can  be 
transformed  into  the  pattern  (I  II  I  II),  either  the  original  pair  is  the 
quad  or  (I,  II)  and  (II,  IV)  are  complementary  pairs. 

1.  Given  that  (I  II  III  IV)  is  transformable  into  (I  II  I  II),  and 
is  so  transformed. 

2.  Lj+LTj+Ljp+Ljy=Ui+Up+UpT+Ujy=  2k^  by  the  supplementary 
property,  Theorem  6.1. 

3.  and  since  ai®a1+j=  ai®^i+j  by  Theorem  2.1. 

4.  Therefore  2Lj+2L,p=  2Uj*f2U p=  2k  ^  by  substituting  into  step  2, 
and  using  the  condition  of  step  1Q 

5.  Lj+Lp=  kj  - 

6.  VUI=  kj  since  the  total  of  likes  and  unlikes  must  equal  all 
pos  sibilitie  s. 

7.  L,jj-Uj=  0  subtracting  the  equation  in  step  6  from  that  in  step  5. 

8.  Therefore  L  j=  Uj  and  (I,  II)  are  complementary. 

9.  Since  (I,  II)  are  complementary,  so  are  (1,7!)  and  any  general 
transformations  of  { 1 , 1 1)  and  (III,  IV). 

=;=  *  * 
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This  last  property  should  be  extremely  useful  in  decomposing  long 
codes,  since  it  avoids  the  tedious  task  of  checking  to  see  if  (I  r  II)  are 
a  complementary  pair.  Whenever  a  complementary  pair  expressed  in 
sequence  quadruple  form  has  the  pattern  (I  II  I  71) ,  certain  additional 
transformations  are  allowed  which  are  not  true  in  general.  These 
allowable  operations  are  the  time  reversal  of  I  and  III,  II  and  IV, 
or  both.  Also  allowed  as  an  operation  is  the  altering  of  either  even  bits 
or  odd  bits  for  all  four  sequences.  The  proof  of  all  these  operations 
being  allowable  is  quite  simple,  but  does  require  an  addition  to  the  no¬ 
tation  to  cover  these  particular  operations.  Double  or  quadruple  sub¬ 
scripts  will  be  used  to  denote  these  operations,  ^13*^24  anc^  ^1^34 
represent  respectively  the  time  reversal  of  (1*111),  (II*  IV)  and  the 
time  reversal  of  both  pairs.  A ^  represents  the  altering  of  the  odd 
bits  in  all  four  sequences  while  A?4  stands  for  the  altering  of  even  bits 
in  all  four  sequences.  Actually  complementing  any  two,  or  all  four  of 
the  four  sequences  will  also  leave  a  complementary  pair  when  reassem¬ 
bled,  but  each  of  these  operations  is  already  accomplished  by  the 
general  transformations. 

Theorem  3.  5 

Whenever  a  complementary  sequence  pair  is  put  into  sequence 
guadruple  and  then  transformed  so  that  the  resultant  group  pattern  is 
(I  II  I  71),  I  and  III,  or  II  and  IV,  or  both,  can  be  time  reversed 
and  the  sequence  pairs  formed  on  reassembly  from  the  quadruple ,  after  the 
operations  (TJ3  T  »4  or  T1  7^)*  will  still  be  complementary  sequence 
pairs. 

1.  Given,  (A,B)  is  a  complementary  pair  of  length  n  transform¬ 
able  in  quadruple  form  to  (I  II  I  II),  and  is  so  transformed. 

2.  (I,  II)  is  a  complementary  sequence  pair  by  Theorem  3.4. 
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3.  Therefore  ( 1^  II) ,  (I,  II)  and  (I,  II)  are  complementary  pairs 

by  Theorem  2e  3. 

4.  From  this  it  is  seen  that  the  original  sequence  quadruple  after 
the  operations  T^3»T24  and  Ti  >34  are  (I  11  J  TIT, 

(I  H  I  II)  and  (I_  12  I_  U) .  When  interlaced  in  the  customary 
manner  to  form  pairs,  the  pairs  they  form  will  be  complementary 
sequence  pairs  by  Theorem  2C  8C 

>}C 

Theorem  3.  6 

Whenever  a  complementary  sequence  pair  is  put  into  sequence  quad¬ 
ruple  and  then  transformed  so  that  the  resultant  group  pattern  is  (I  II 
I  TT),  the  operation  A,  the  altering  of  odd  bits,  or  the  operation 

the  altering  of  even  bits,  will  allow  the  quadruple  upon  reassembly 
to  still  be  a  complementary  pair. 

1.  By  hypothesis  A  and  B  are  a  complementary  pair  which 
when  expressed  in  sequence  quadruple  form  are  transformable 
to  (I  II  I  TT),  and  is  so  transformed. 

2.  (I,  II)  is  a  complementary  sequence  pair  by  Theorem  30  4. 

3.  (I,  II)  when  operated  upon  by  either  AT  or  A7  is  still  a 
complementary  pair  by  Theorem  2.  4U 

4.  If  the  symbol  (  1,  2)  is  used  for  the  pair  (I,  II)  altered  it  is 
seen  that  the  pair  formed  from  the  interlace  of  (12  12)  is 
a  complementary  pair  by  Theorem  2.  8C 

>*e  *  * 

These  last  two  theorems  should  prove  quite  useful  when  trying  to 
estimate  an  upper  limit  on  the  possible  number  of  composite  comple¬ 
mentary  sequences  of  any  particular  length.  An  example  to  illustrate 
the  use  of  Theorem  3C  5  might  prove  helpful*  Taking  pair  number  one 
from  Appendix  II, 
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A=1 101 0001  1101  1110 
B=1000010010001011 
1=10001011=1 
11=01  111011  =  11 
111=1000101 1=1 
IV  =  1  0000100  =.ir 

Since  this  is  already  in  (I  II  I  II)  form  without  transforming  it, 

can  be  applied  directly  giving 

1=1 1010001 
11=01 11 1011 
111=11010001 
IV  =  10000100. 

Reassembling  gives 

A1= 11 1 1001101 0101 10 
B:= 101001 1000000011 

which  is  identical  with  number  three  in  Appendix  II. 

As  an  example,  to  illustrate  the  use  of  Theorem  3.  6,  applying  A^ 

to  the  same  (A,B)  pair  gives 

1. =11 01 11  10 
2. =00101  1  10 
3.  =  110111  10 
2\  =  1 1  01 0001 , 

reassembling  give 

A2=101101 1110111000 
B2=1110001011 101101. 

This  is  the  same  as  number  6  in  Appendix  II  with  the  transformation 
TCs=£)  applied. 
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T  ransform 

Unorde  red 
Equivale  nt 

T  ransform 

Unordered 

Equivalent 

1. 

i  =  (i  ii  i  ny 

a2 

17. 

R  = 

(l  Ti  Ti  T) 

O 

2. 

a1=(1  ii  r  ny 

C 

18. 

Q  = 

(ii  T  n  T) 

J 

3. 

a2=(i  IT  l  ii) 

1 

19. 

P  = 

(TT  I  ii  l) 

T 

4. 

T!=(ii  l  l  TT) 

s 

20. 

O  = 

(ITT  1  TT) 

R 

5. 

t2=(i  ii  IT  I) 

X 

21. 

N  = 

(ii  l  T  ii) 

U 

6. 

T  =(II  I  TT  1) 

p 

22. 

M  = 

(I  n  Ti  i) 

Y 

7. 

C!=(I  IT  i  IT) 

F 

23. 

L  = 

(Ti  T  T  ii) 

V 

8. 

C2=(l  11  1  11) 

D 

24. 

K  = 

(T  Ti  n  T) 

Z 

9. 

C  =(T  TT  1  ii) 

Ax 

25. 

J  = 

Til  T  ii  T) 

Q 

10. 

z  =(ii  T  T  IT)) 

K 

26. 

H  = 

(l  n  ill 

W 

1  1. 

Y  =(TT  i  1  in 

M 

2  7. 

G  = 

in  ITT  l) 

B 

12. 

x  =(TT  l  i  U) 

t3 

28. 

F  = 

(i  n  T  Tij 

Ci 

13. 

W  =(I1  1  I  11) 

H 

29. 

D  = 

(T  II  l  II) 

c3 

14. 

v  =(T  n  TT  I) 

L 

30. 

IT  = 

(II  I  II  T) 

0 

15. 

u  =(I  ii  ii  l) 

N 

31. 

B  = 

(TT  I  TT  T) 

G 

16. 

s  =(i  TT  n  i) 

Ti 

32. 

0  = 

(II  T  ii  l) 

IT 

E  =  (I  II  I  ID  a2,i 


Redundance  due  to  (1  II  I  II). 
TABLE  3.4 
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The  (I  II  I  II)  form,  or  one  of  its  transforms,  appears  quite 
often  in  the  study  of  composite  complementary  sequences.  Since  two 
of  the  components  are  the  same,  and  the  other  two  are  complements, 
it  is  reasonable  to  assume  that  some  of  the  transforms  are  redundant, 
that  is  yield  indistinct  results.  Table  3C  4  lists  all  the  transforms 
with  the  elements  I,  II  and  their  complements,  If  (A9  B)  are  taken 
as  an  ordered  pair,  the  number  of  different  codes  formed  from  the 
transformation  group  are  32  instead  of  64  because  A2(I  II  I  II) 
and  E(I  II  I  TIj  are  the  same.  If  (A,B)  are  taken  as  an  unordered 
pair  the  number  of  transforms  forming  different  codes  is  decreased 
from  32  to  16.  Table  3.4  lists  each  transform  with  its  unordered 
mate. 

The  only  possibility  of  having  less  than  16  unordered  pairs  of 
codes  which  are  distinct  would  be  if  the  group  form  could  be  transformed 
into  the  pattern  (I  1  I  I).  The  only  case  known  to  the  writer  where 
this  is  true  is  the  quad.  For  one  version  of  the  quad  this  is 

A=  1 1  B= 10 

1=1 
11=1 
111=1 
IV =0. 

Table  3.  6  shows  the  eight  ordered  pairs,  or  if  the  cross  index  column 
is  used  the  four  unordered  pairs,  for  the  quad.  The  four  possible  code 
transformations  for  the  quad  are: 

12  3  4 

A=  1  1  1  1  00  00 

B=  10  01  10  01. 
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Patte  rn 

Transforms  with  pattern 
for  ordered  pairs 

Cross 
for  ui 

1. 

(i  i 

I 

1) 

I.Tx.H,  it 

3. 

2. 

(i  i 

T 

T) 

AlP  Y,V,  B 

6. 

3. 

(i  1 

i 

i) 

Aa,W,S,0,  E 

1. 

4. 

(i  i 

1 

i) 

t3,t,c3,n 

8. 

5. 

(T1 

i 

I) 

C1pO,K,  J 

7. 

6. 

(T  T 

T 

i) 

C,M,  L,G 

2. 

7. 

(i  1 

T 

i) 

Z,R,Q,F 

f. 

8. 

(T  i 

i 

i) 

X, U,P,D 

4. 

Redundance  due  to  (I  I  I  7)  . 
TABLE  3.5 


This  exhausts  the  transformations  of  special  types  and  the  redun¬ 
dancies  caused  by  (I  II  I  II).  The  operations  group  will  now  be 
compared  to  Senior’s  group  number  44  of  length  32  as  was  mentioned 
earlier  in  this  chapter. 

The  identification  of  any  group  with  a  group  in  Senior’s  list,  is 

accomplished  through  establishing  a  one  to  one  correspondence  between 

16  As  will  be  demonstrated,  Senior’s  group  num- 

group  generators. 

ber  44  of  order  32  is  in  one  to  one  correspondence  with  the  unordered 
operations  group  of  complementary  sequences  expressed  in  sequence 
quadruple  form.  Therefore  the  two  groups  are  isomorphic. 

Senior’s  generators  will  be  given  in  lower  case  letters,  while  those 
of  the  operations  group  will  be  in  their  standard  upper  case  symbols. 
Five  generators  are  required  for  Senior’s  number  44;  let  these  be 
a  ,  b ,  c  ,  d ,  e .  The  following  relationships  between  generators  are  re¬ 
quired  for  Senior’s  44, 

a2=b2=c2=d2=e2  =  1  where  1  is  the  identity. 

b~kb=a,  c-1ac=a,  d~1ad=a  where  x  1  is  the  inverse  of  x. 


1. 

2. 
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3. 

c^bcsb, 

d  1bd=b , 

4. 

d  1cd=c, 

5. 

e  1ae=a, 

e  1be=b , 

6. 

e  1ce=ca 

,  e  1de=bd 

A  check  of  the  above  requirements  shows  that  a,b,c,d  form  an 
abelian  group  of  order  16  while  the  Merl  generator  is  commutative 
only  with  a,b  and  is  non- commutative  with  c,d. 

Ltet  a=Ci,  b=C2,  c=Tx,  d=Ts,  e  =  Ax  »  From  Table  3.3  it  is  seen 
that  the  following  relations  are  true. 

1.  C?  =  C 1=  T?  =  Tf=  Af  =  l 

2 .  C 2C 1  C3  =  Ci,  T 1  Ci  T:  =  Ci  ,  T3  Ci  Tg  =  Ci  . 

3.  Tx  C2  Tj.  =  Ca  . 

4.  T3  To.  T3  =  Ti  . 

5.  A]_Ci  Ai  “  Gi  ,  Ax  C2  Ax  =  C2  „ 

6.  Ax  Tx  Ax  — •  Tx  Cx  ~  Ax  T 2  Ax  ~  C2  T 2  —  H  ,, 

This  one  to  one  correspondence  between  the  generators  of  Senior's 
group  number  44  and  a  set  of  generators  from  the  unordered  operations 
group  on  complementary  sequences  is  necessary  and  sufficient  to  prove 
the  two  groups  are  isomorphic* 

The  diagrams  for  Senior's  group  number  44  of  order  32  and 
those  for  its  subgroups  were  used  as  a  guide  line  to  check  the  generation 
of  subgroups  from  the  operations  group  of  unordered  complementary 
sequence  pairs.  These  subgroups  of  order  1,2, 4, 8  and  16  are  shown 
in  Appendix  IV*  All  groups  that  were  shown  in  Senior  s  diagrams  are 
included  in  these  charts;  however,  Senior’s  diagrams  included  only 
normal  subgroups  and,  although  much  care  was  taken  in  the  generations 
of  the  subgroups  of  the  operations  group,  there  is  a  possibility  that 
Appendix  IV  is  not  exhaustive  due  to  some  undetected  non- normal  sub¬ 


groups. 
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Senior  in  conjunction  with  others  has  rewritten  his  work  and 

17 

changed  his  designating  method.  Under  this  new  method  his 

number  44  of  length  32  is  now  j  •  The  ordered  operations  group 

of  length  64  is  also  identified  as 

An  identification  of  the  group  with  Millerfs  number  23  of  order 

32  was  also  accomplished  for  those  who  prefer  the  substitution  group 

1  8 

method  of  identification* 
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CHAPTER  IV 


HAMMING  DISTANCES  OF  COMPLEMENTARY  SEQUENCES 

The  Hamming  distance  is  a  well  known  property  of  certain  classes 
of  codes.  ^  This  very  useful  property  has  had  most  of  its  applica¬ 

tions  in  the  fjelds  of  binary  error  detecting  and  error  correcting  codes. 
This  chapter  will  demonstrate  the  usefulness  of  Hamming  distances  in 
the  field  of  complementary  sequences,  as  both  the  code  pairs  and  their 
sequence  quadruple  formulation  have  certain  invariances  in  Hamming 
distance  under  various  conditions. 

The  Hamming  distance,  D(US  V),  of  two  binary  vectors  or  sequences 
U  and  V  is  defined  as  the  number  of  positions  in  which  these  two  binary 
vectors  differ. 

i=n 

D(U,  V)  =  E  u.  9  v.  (4.  1) 

i=i  1 

For  example,  the  Hamming  distance  of  the  two  binary  vectors 

U=1011001 

V=01  11001 
7 

is  D(  U ,  V)  =  S  u.  $  v  =  (1  +  1+0+0+0+0+0)  =  2 
i=  1  1  1 

The  Hamming  weight  or  just  weight  of  a  binary  vector  is  its  distance 
from  the  null  vector. 

Theorem  4  1 

All  complementary  sequence  pairs  of  length  n  have  a  Hamming 
distance  of  n/Z. 

By  hypothesis 

A^=ala2a3c  .  •  •  •  •  an 

and 

B=blb2b3’  "  •  '  *  '  bn 

A  and  B  form  a.  complementary  sequence  pair. 


49 


1. 


z. 


3. 


All  complementary  sequence  pairs  satisfy  the  parity  check 

ai§an-i+l  ®bi  ®bn-i  +  l  =  1'£ora11 

ii  1  <  i  <  n/2  by  equation  2.  4„ 

Since  mod  2  addition  is  commutative 

a.  »b.  ©a  .  ©b  .X1  =  1  . 
l  l  n-i+1  n-i-fl 

Therefore  either  (a.  ©b  )  =*  l  or  (a  M  ®b  ,  , )  ,  but  not 

ii  n-i  +  1  n-i-fT' 

both. 

For  each  i,  1  <  i  <  n/2  one  of  the  two  vector  positions  in 
step  3  will  have  a  one  and  the  other  vector  position  a  zero. 

Summing  up  over  all  possible  i,  gives  n/2  for  the 
Hamming  distance. 


*  *  * 

When  a  complementary  sequence  pair  is  decomposed  into  a  sequence 

quadruple,  each  pair  of  sequences  in  the  quadruple,  which  form  one  of 

the  complementary  pairs,  has  a  Hamming  distance  which  will  be  shown 

to  be  characteristic  of  the  code  lengthc  The  sum  of  the  two  pairs  of 

Hamming  distances  from  the  quadruple  is  an  invariant  and  is  equal  to 

n/2  the  Hamming  distance  of  the  complementary  sequence  pair. 

For  example, 

A=  1001 010001 
B=  1000000 1  10 

are  a  complementary  sequence  pair  of  length  10. 

Decomposing  them  into  the  standard  sequence  quadruple,  (I  II  III  IV), 
give  s 

I  =  10000 

II  =  10110 

III  =  10001 

=  01000 


IV 
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I  ®  II  *  00110  Dll'll)  =  2 

III  $  IV  =  11001  D(TII,  IV)  =  3 

i  ®  ii  «  hi  ®  iv  =11111  d(i«ii,  iii®iv)  =  5  . 

This  distance  is  half  the  length  of  the  complementary  sequence  pair, 
and  is  required  because  of  the  necessary  condition  of  the  parity  check. 

If  the  distances  (2,  3)  for  the  pairs  given  in  the  above  example  are 
taken  as  an  unordered  rathe  r  than  an  ordered  pair  ,  this  (2,3)  pair 
will  be  invariant  through  all  general  transformation  operations. 

All  of  the  general  transformations  on  the  sequence  quadruples 
have  one  or  more  of  the  following  properties: 

1.  Change  the  order  of  (I,  II)  or  (III9IV)  or  both, 
for  example,  T1  ,  T2,  T  . 

2.  Complement  the  pairs  (I,  II)  or  III 9  IV)  or  both, 
for  example,  C\  ,  C2  *  C  . 

3.  Complement  one  from  the  pair  (I,  II)  the  other  f  rom 
the  pair  (III,  IV),  for  example,  A1  ,  A2  . 

4.  Exchange  (I,  II)  for  (III,  IV),  for  example,  E 

Since  all  operations  in  the  group  can  be  generated  by  multiplication 
from  the  set  listed  above,  an  examination  of  the  effect  on  Hamming 
distance  by  these  transformations  will  be  adequate  to  prove  the  mvar- 
ianc  e. 

Theore  m  4.2 

A  change  in  order  of  (I,  II)  or  (III,  IV)  or  both  will  not  change 

the  Hamming  distances  of  the  pairs  (I, II)  or  (II1,IV)0 

n 

1.  By  definition  D{U,V)  =  E  u  ©  w 

i“  1 

2.  u.  ©v.  =  v.  ®u.  modulo  2  addition  is  commutative 

1111 

i~n  i  =  n 

3.  DfU.V)  =  E  u  ®v.  =  E  v.  ®u  =  D( V  ,U) 

i=l  i=l 

*  *  * 
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Complementing  (I,  II)  or  (III,  IV)  or  both  will  not  change  the 
Hamming  distances  of  the  pairs  (I,  II)  or  (III,  IV). 

1.  u.  9  v.  =  u.  ®  v.  by  Theorem  2.  1 

1111  y 

2.  Therefore  E  u.  #  v.  =  E  u.  ©  v. 

ii  ii 

3.  and  D(U  ,  V)  =  D(U,V) 

❖  *  * 

If  it  is  as  sumed  that  the  length  of  the  sequence  quadruples  formed 
from  the  complementary  sequence  pairs  have  a  length  K  and  that  the 
distance  of  the  first  pair  is  R,  then  the  following  theorem  can  be  proved: 
Theorem  4.  4 

a 

Complementing  one  of  the  pair  (I,  II)  and  one  of  the  pair  (III,  IV) 
results  in  the  exchange  of  their  Hamming  distances,  or  expressed  dif¬ 
ferently,  gives  the  same  distance  pair  in  opposite  order. 

1.  Let  D(I,II)  =  R,  therefore 

D(III,  IV)  =  K-R 

as  the  parity  check  requires  D(I®TI,  III  ®IV)=K 

2.  Complement  one  of  the  first  pair  in  the  ( I  II  III  IV)  form, 
and  also  complement  one  of  the  second  pair. 

Each  vector  position  that  was  formerly  alike  and  summed 
to  zero  is  now  different  and  sums  to  one. 

Similarly  each  vector  position  that  was  formerly  different 
and  summed  to  one  is  now  alike  and  sums  to  zero. 

3.  D(I,II)  =  D(I,II)  =  K-R  from  2. 

4.  D(7lLIV)  =  D(III,IV)  =  K-(K-R)  =R  from  2. 

*  ❖ 

Theorem  4.  5 

The  exchange  operation  has  the  effect  of  changing  the  order  of 
Hamming  distances  of  the  pairs  (I, II)  and  (III, IV)  when  expressed 
in  sequence  quadruple  form. 


1  .  The  pairs  (i,  II)  and  (III,  IV)  are  not  changed  in  any  way 
except  by  position  in  the  exchange  operation* 

2.  The  only  way  to  effect  the  Hamming  distances  of  a  pair  of 
vectors  is  to  change  some  of  the  components  of  the  vectors. 

3.  Therefore  the  two  weights  are  the  same  but  reversed  in 
order. 


*  *  * 

An  examination  of  the  Hamming  distances  for  the  known  kernels  of 
2 ,  10  and  26  is  now  in  order  as  these  distances  are  a  characteristic 
of  the  kernel,  as  are  the  vectors  formed  in  the  process  of  obtaining  the 
distances  by  the  modulo  2  sums. 

The  first  case  to  examine  is  the  trivial  case  for  the  kernel  of 

■5 

length  2. 

A=  1 1  B= 1 0 

Broken  down  into  sequence  quadruple  representation 

1=1  D(1,H)  =  0 

11=1 

111=1  D(  III,  IV)  =  1 

IV  =  0 

There  are  two  kernels  for  length  10  to  be  examined. 

A=1001  010001 
B=  10000001  10 


1=10000 

I  © 

11=001  10 

D(I,  II)=  2 

11=101  10 

111=10001 

III  © 

IV=11001 

D(III,  IV)=  3 

IV=01000 

A=01 01 0000 1 1 
B  =  0000  1001  10 

1=00001 

11=10011 

I  © 

11=10010 

D(I,  II)=  2 

111=00101 

III  © 

IV=01  101 

D( III,  IV)=  3 
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Both  kernels  of  length  10  have  the  unordered  Hamming  distances  of 
(2,3)  when  expressed  in  sequence  quadruple  representation.  However, 
their  "Hamir  ing  vectors'1,  the  unordered  pairs  (00110,11001)  and 
(10  010,01101)  are  different  and  are  characteristic  of  the  kernel  of 
their  origin  under  all  possible  transformation  operations.  Hamming 
vectors  will  be  considered  at  length  in  the  next  chapter  and  proof  of  the 
above  statement  will  be  deferred  until  then. 

There  is  only  one  kernel  of  length  26  (Chapter  7)  and  it  is 
A=0100 1101 1110101 11100111010 

b=ioi  iooiooooi  mi  nooi  noio 

1=0010111110111  I  ©  11=0000001  101010  D(  1,  II)  =  4 

11=00  101  10011  101 

111=1101001110111  III  ©  I V=  1111110010101  D(III,  IV)=  9 

IV=001011 1100010 

This  is  a  complete  listing  of  the  Hamming  distances  for  all  known 
kerne  Is  and  it  is  exhaustive  for  n=2,10,18  and  26,  It  is  worthwhile 
to  note  that  the  sum  of  the  unordered  pairs  of  distances  is  always  equal 
to  one  half  the  length  of  the  code  and  that  the  Hamming  vectors  of  the 
pairs  are  complements.  Both  of  these  conditions  are  due  to  the  parity 
check  being  a  necessary  condition,, 

A  natural  extension  of  the  Hamming  distances  of  kernels  is  the 
consideration  of  Hamming  distances  of  composite  codes  when  put  into 
sequence  quadruple  form.  Six  different  code  pairs  of  length  eight  are 
listed  in  Table  4.  1,  These  pairs  will  be  examined  for  a  possible  in¬ 
variance.  After  this  invariance  is  noted  a  general  proof  will  be  given 
for  the  characteristic  distance  of  composite  codes. 

The  first  four  examples  in  the  table  when  decomposed  into  sequence 

quadruple  representation  have  (I, II)  and  (III,  IV)  as  code  pairs  which 

are  still  complementary.  These  codes  are  half  the  length  of  the  original 
pair,  and  being  complementary,  have  a  Hamming  distance  one  half  of 
their  length,  or  n /4  as  compared  to  the  original  codes  of  length  eight, 
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(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


A=llll 1001 
B=  11001010 

1=11 10 
11=101 1 

111=101  1 
IV=0001 

A=  11110110 
B=1 1000101 

1=  1 1  0 1 
11=01  11 

111=1000 

IV=1101 

A= 11101101 
B=1 1 100010 

1=11  10 
11=1101 

111=  1101 
IV=0001 

A= 1101  1110 
B=  1101  000  1 

1=101  1 
11=01  11 

111=1000 

IV=1011 

A=1 1101011 
B=ll 100100 

1=  1 1  1 1 
11=10  11 

111=1100 

IV=0101 

A=1011 1110 
B=101 10001 

1=  1 1  1 1 
11=01  10 

111=1100 

IV=1010 


H=(00  110011) 
I®  11=0101 

III  ®  IV=1010 

H=  (00  110011) 

I  ®  11=1010 

III  ®  IV=0101 

H=(000011  11) 

I  ®  II  =001  1 

III  ®  IV  =  1 100 

H=(000011  11) 

I  ®  II  =  1  100 

III  ®  IV  =  001 1 

H=(000011  11) 

I  ®  II  =  0110 

III  «  IV  =  1001 

H=(000011  11) 

I  ®  11=1001 

III  ®  IV=0110 


D( I,  II)  =  2 
D(III,  IV)  =  2 


D(I,  II)  =  2 
D(III,  IV)  =  .  2 


D(I,  II)  =  2 
D( III,  IV)  =  2 


D(I,  II)  =  2 
D(  III ,  IV)  =  2 


D(  I,  II)  =  2 
D( III,  IV)  =  2 


D(I,  II)  =  2 
D(III,  IV)  =  2 


Some  Hamming  distances  and  vectors 
of  codes  of  length  8. 


TABLE  4„  1 


tt 
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An  examination  of  the  last  two  examples  in  the  table  shows  that 
although  the  (I,  II)  and  (III,  IV)  pairs  are  not  complementary  their 
Hamming  distance  is  still  n/4  when  compared  to  the  original  codes 
of  length  eight.  That  this  feature  is  true  in  general  will  now  be  proven 
by  examining  and  exhausting  all  the  known  ways  of  generating  composite 
code  s . 

Theo re m  4,6 

Every  composite  code  of  length  n  (generated  by  one  of  the  standard 
methods)  when  expressed  in  the  sequence  quadruple  representation  will 
have  D(I,  II)  =  D(III,  IV)  =  n/4. 

There  are  presently  known  four  general  ways  of  generating  com¬ 
posite  codes  from  shorter  cqdes,  (there  is  also  the  method  applicable 
only  to  codes  of  length  2r  given  at  the  end  of  the  proof).  These  are: 

1.  Time  Sequence 
Sx  =  AB 
S2  =  AB 


2. 


Inte  rlace 

Tx  =  a  t  bn  a-jb.-,.  .  „  a_b 


112  2* 


n  n 


Ta  =  a  ,’E,  a^bp.  .  .  a  b 
1  1  c*  c.  n  n 

3.  Time  Sequence  Exponential 


C1  .  C2 


Ux  =  A  A 


A 


cmAl  d2 


B 


m 


d  d  , 

U3  =  A  mA  m:...  A  "B  lliB 


dl  T-,Cm„cm-  1 


.  .B 


4.  Interlace  Exponential 

C1  dl  C2  d2 
V1  =  A  B  A  B  .  .  _ 

d  c  d  .  c  , 

V2  =  A  mB  mA  m" 1 B  m_1 


c  d 
A  mB  m 


dl  Cl 
...A  B 


Throughout  this  proof  it  will  be  considered  that  (A,B)  are  a  comple¬ 
mentary  pair  of  length  r  and  (C,D)  are  a  complementary  pair  of 
length  m. 


=  1  is  a  necessary 


The  parity  check  a.  ©  a  ,  ,  ©  b  ®  b  ,  . 

i  n+l-i  1  n+l-i 

condition  for  a  sequence  pair  to  be  complementary.  One  of  the  simplest 

methods  for  making  the  parity  check  is  to  fold  both  codes  in  the  middle 

9 

and  double  them  back  on  themselves.  This  procedure  is  shown  in 
Figure  4.1.  Each 


a  ^a^a^.  . 

l 

- an/2-  ian/2 

a  a  !  .  .  .  . 

n  n-  1 

n+l-l 

*  *  *  ’  an/2+lan/2+l 

blb2b3*  *  *  * 

.  .  b. 

1 

*  ‘  bn/2-  lbn/2 

b  b  -  1  , 
n  n 

••  bn+l-i 

*  '  *  ‘  bn/2  +  2bn/2+l 

Foldover  Parity  Check  Method 
Figure  4.M 

column  in  Figure  4.  1  satisfies  the  parity  check.  If  code  pairs  are 
written  in  the  sequence  quadruple  form  such  as  in  Figure  4.  2  ,  it  is 
seen  that  the  columns  of  this  array  are  the  same  as  the  column  in  the 
first  array  with  rearrangements  in  order  provided  it  is  allowed  that 
the  columns 


a  and  a 

u  r 


b  b 

r  u 


b  b 

u  r 


are  the  same.  Since  modulo  2  addition  is  the  operation  used  in 
making  the  parity  check  and  in  determining  Hamming  distance,  and  the 
use  of  each  column  involves  only  modulo  2  addition  which  is  commu¬ 
tative,  the  two  columns  can  be  considered  the  same.  Therefore  the 
fold-over  method  of  grouping  and  the  sequence  quadruple  form  of  grouping 

have  invariant  columns,  order  being  ignored. 

57 


la3a5... 

an-  1 

n  n-  2  *  * 

a2 

lb3 

. -bn-3 

bn-l 

b  ?  . 

n  n-  2 

. b4 

b  2 

Sequence  Quadruple  Form 
of  Parity  Check 

FIGURE  4.  2 


Now  consider  the  composite  generating  methods. 

1.  Time  Sequence 

51  =  AB 

52  =AS  writing  in  the  fold- over  parity  check  form 

a  i  a9a"3 . 3,  a  i  a. 

1  2  3  r-2  r- I  r  A 

b  b  . b  b?  b,  D<A»B)  " 

r  r-  1  3  2  1 

a,a.a, .  a  0a  ,a 

1  ^  3  r-2  r-  1  r  T ^ /  A  lo 

bS  .  b,  b,  E.  D(A.B)  =  r/2 

r  r-  1  3  2  1 

The  top  pair  are  complementary  and  therefore  have  a  distance  of 
r/2,  similarly  the  bottom  two  rows  are  complementary  and  have  a  dis¬ 
tance  of  r/2  by  Theorem  40  1.  To  write  the  time  sequence  code  in 
sequence  quadruple  form  would  not  change  the  distance  since  it  would  be 
purely  a  reordering  of  the  parity  check  columns*  Therefore  D (I ,  II) 
would  equal  r/2  as  would  D(UI,>IV)B  In  terms  of  the  composite  code 
of  length  n  D(I,  II)=D(III,  IV)=  r/2  =  n/4  whether  or  not  (I,  II)  and 
(III,  rV)  happen  to  form  complementary  pairsB 

2.  For  the  interlace  scheme  the  (I,  II)  pair  and  the  (III,  IV)  pair 
are  obviously  complementary  sequences  by  construction  and  will  have  a 
distance  of  r/2  =  n/4. 

3.  Writing  the  time  sequence  exponential  pair  in  fold-over  form 

§ives  U1=AClAC2.....  ACm 

d.  d  j  where  the  symbol  B  means  B  time 

r  mR  m"  Ri  reversed. 
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u2=a  mA 

( 

.  .  .  A 

c  1  C  2 

B  B  . 

\ 

.  .  .  B 

The  B's  are  all  time  reversed  but  the  time  reversal  of  one  of  the 
complementary  pairs  does  not  change  the  complementary  property 
D(A ,  B  )=D(A,  B)  =  r/2.  For  each  bit  of  the  C  code  and  D  code  there 
is  a  distance  r/2  since  every  A  is  matched  with  some  transformation 
of  B  for  each  of  these  bits.  There  are  a  total  of  m  bits  in  the  C 
code  and  in  the  D  code,  therefore  the  total  distance  of  folded  back 

on  itself  is  m  x  r/Z  =  mr/2.  Similar  reasoning  gives  U3  folded  back 
on  itself  a  distance  of  mr/2.  A  rearrangement  of  the  columns  will  not 
change  the  Hamming  distances  Therefore  D(J,J1)=mr/2  and 
D(III,  IV)  =  mr /2 ,  since  the  total  length  of  the  code  is  n=2mr,  D(I,  II)  = 
D(III,  IV)=n/4. 

The  interlace  exponential  scheme  is  the  last  general  method  of 
code  generation  to  be  considered.  Written  in  the  fold-over  form  it  gives 


C1  dl  c2  d2 
V1=A  B  A  B  .  .  . 

^Cm/2gdm/2 

d  c  d  , 

B  mA  mB  m'A  .  . 

g^m/  2+1  ^cm/  2+1 

d  i  d  . 

v2=a  mB  mA 

^  m/  2+lgCm/  2  +  1 

c 1  dl  c  2 

B  A  B 

BCm/2  Adm/2 

Again  each  A  is  matched  with  some  transform  of  B,  D(  A,  B)  =  D( A,  B)= 
D(A,B),  etc.  =  r/2.  As  before,  rearranging  columns  into  sequence 
quadruple  form  will  not  effect  the  Hamming  distances.  The  distance  of 
the  top  two  rows  is  mr/2  and  the  bottom  two  rows  (V2)  also  is  equal 
to  mr/2.  The  total  length  of  the  code  is  2mr,  therefore 
D(I,  II)=D(III,  IV)=mr/2=n/4. 

Y 

5.  When  composite  codes  are  of  length  2  ,  they  are  sometimes 
formed  in  a  special  way  from  the  interlaced  sections  of  a  complementary 
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pair.  These  sections  or  pieces  of  the  codes  must  be  of  length  p=2 
where  j<r.  Using  the  symbols  A  =a ^ a^a^a^.  .  .  .  a 


A2p~ap+lap+2* 


a^p  and  so,  the  following  pair  would  be  complementary  if  A  and  B 


formed  a  complementary  pair. 


S t  =A  B  A^ 

lp  p  p  2p  2p 


=A  B  A.  B. 

Zp  p  p  2p  2p 


Anp/pBnp/p 

Anp/pBnp/p 


Rather  than  try  to  use  the  general  term  to  prove  that  D(1, 1I)=D (III,  IV ) 

=  n/4  as  was  done  in  the  previous  cases,  a  short  general  example  of 
length  8  will  be  used.  These  two  codes  (A,B)  of  length  eight  will  be 
interlaced  four  bits  at  a  time  (p  =  2)  to  form  a  code  of  length  16. 

S14=ala2a3a4blbZb3b4a5a6a7a8b5b6b7b8 
S24=ala2a3a4blb2b3b4a5afea7a8b556B7b8 
Writing  this  composite  pair  in  fold-over  form  yields 

S14=ala2a3a4blb2b3b4 

b8b7b6b5a8a7a6a5 

S24=ala2a3a4blb2b3b4 

b8b766b5a8a7a6a5 

D(A,B)=4  by  Theorem  4.1.  folded  over  is  just  a  column 

rearrangement  of  the  complementary  pair  (A,B),  and  since  a  re¬ 
arrangement  of  columns  does  not  effect  Hamming  distance,  the  Hamming 
distance  of  folded-over  is  also  4.  Similar  reasoning  holds  for  S-,^ 

Since  the  standard  sequence  quadruple  form  is  again  just  a  re¬ 
ordering  of  the  columns,  the  distance  of  the  (1,11)  pair  and  the 
(III,  IV)  pair  is  n/4.  For  this  example, 


1=ala3blb3a5a7b5b7 

II=b8b6a8a6b4b2a4a2 


D( I,  H)=4 
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IIl=a1a3b1b3a5a7b5b7 


D(1II,  IV)=4 


IV=S8S6a8a6B482a4a2 


D(I©II ,  III®IV)=8. 


Therefore  any  composite  complementary  sequence  pair  when  written  in 
sequence  quadruple  form  has  D(  1, II)=D(III,  IV)=n /4. 


*  *  * 

Although  there  is  no  proof  that  the  4  standard  generating  methods, 

j* 

plus  the  special  one  applicable  only  to  2  codes,  are  the  only  way  to 
form  complementary  codes  which  are  composite,  an  exhaustive  search 
of  code  lengths  16  and  20  revealed  no  codes  which  were  not  formed 
by  these  standard  generating  methods. 

The  full  usefulness  of  the'J  Hamming  distance  property  is,  of  course, 
not  yet  known.  The  property  of  the  distance  of  the  complementary  pair 
being  equal  to  half  the  code  length  as  shown  in  Theorem  4.  1  ,  is  used 
in  many  of  the  proofs  in  the  chapter  on  Hamming  vectors.  This  property, 
D(A,B)=n/2,  indicates  that  the  code  pairs  are  orthogonal,  Because  of 
this  orthogonality  property  the  leakage  of  an  A  code  carrier  into  the 
B  code  filter,  or  vice  versa,  will  have  no  effect  on  the  output  at  exact 
match. 

The  primary  use  the  writer  has  found  for  the  Hamming  distance  of 

the  (I,  II)  and  (III,  IV)  pairs,  has  been  the  error  check  that  this  invariant 

property  offered  in  decomposing  codes  by  hand.  One  possible  use  for  the 

invariance  of  D(I,II)  and  D(II1,IV)  could  be  in  the  search  for  kernels  of 

new  lengths.  If  a  characteristic  Hamming  distance  for  unknown  kernels 
could  in  some  way  be  deduced,  based  only  on  length,  this  invariant  prop¬ 
erty  of  distance  would  act  as  a  most  powerful  screen  in  the  search.  The 
(I,  II)  and  (III,  IV)  pairs  are  themselves  orthogonal,  unless  the  code  is  a 
kernel.  This  may  provide  other  possible  methods  for  searching  for  new 
codes.  It  is  difficult  to  say  in  advance  what  utilization  will  be  made  of  any 

invariant  property,  but  it  is  fairly  safe  to  say  that  the  property  will  find  use. 
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CHAPTER  V 

THE  HAMMING  VECTORS  OF  COMPLEMENTARY  SEQUENCES 

The  Hamming  distance  between  the  code  pairs  in  a  complementary 
sequence  was  shown  by  Theorem  4.  1  to  be  one  half  the  length  of  the 
code.  This  Hamming  distance  was  obtained  by  first  finding  what  I  shall 
call  the  Hamming  vector,  H(A,B)  =  A  ®  B  or  for  convenience  just 
H=A®B,  and  then  counting  the  number  of  ones  H  contains.  The  Ham¬ 
ming  vector  is  itself  useful  in  recognizing  the  composition  of  a  composite 
code  when  it  is  desired  to  break  the  code  down  into  shorter  lengths.  This 
decomposition  might  be  used  to  find  previously  unknown  kernels,  or  to 
set  an  upper  bound  on  the  possible  number  of  pairs  of  codes  of  any  length. 

An  examination  of  each  of  the  methods  for  generating  longer  codes 
from  shorter  ones  as  explained  in  Chapter  2  will  show  that  each  method 
has  a  characteristic  Hamming  vector  form,  although  it  may  be  necessary 
to  time  reverse  one  of  the  codes  to  obtain  this  form.  The  method  of 
generating  codes  of  length  2n,  which  does  not  apply  to  complementary 
codes  in  general,  will  also  be  studied  with  its  characteristic  Hamming 
vector. 

In  the  following  proofs  (A,B)  form  a  complementary  code  pair  of 
length  n  while  (C,D)  form  a  complementary  code  pair  of  length  m. 

The  symbol  (0)n  or  just  0n  will  mean  a  string  of  n  zeros,  similarly 
(01  )n  will  symbolize  a  string  of  n  zero  ones,  and  (a^®b^)n  will  be  a 
string  of  n  elements  equal  to  (a  $b_.)  # 

T  heo  re  m  5,1 

The  time  sequence  form  of  generating  a  composite  code  has  the 

characteristic  Hamming  vector  H=(0)n(l)  4 

Use  the  time  sequence  generating  method  on  (A,  B)  to  form 
Si=  AB 
Sg=  AB  . 
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Take  the  Hamming  vector  of  (Si,  Ss),  H=Si  ©  Ss= 
seen  to  be  true  by  inspection  and  definition  of  modulo 

*  *  * 


0n  1 n ,  thi sis 

( 5.  1) 

2  sum. 


Theorem  5.  2 

The  interlace  form  of  generating  a  composite  code  has  the  charac¬ 
teristic  Hamming  vector  H=(01)n. 

Use  the  interlace  generating  method  on  (A,B)  to  form 

T ^  • a i  bT  a.b^a^b^.  .  .  .  .  .  .  .  .  a  b 
^  1  I  Z  Z  3  3  nn 

T 2=  a1b1a^b-}a~b~........„a  S  * 

I1ZZ33  nn 

Take  the  Hamming  vector  of  (Ti,  T2),  H=Ti  ©  T2=(01)n,  (5.  Z) 

this  is  seen  to  be  true  by  inspection  and  definition  of  modulo  Z  sum. 

>\<  >jc 


Theorem  5.  3 

The  time  sequence  exponential  form  of  generating  a  composite  code 

has  the  characteristic  Hamming  vector  H=(c,©d  )n  (c^®d  ,  )n  ... 

1  m  Z  m- 1 


(c  ®d1)n  (d,  ®c  )n.  ...d  ®c,  )n . 

mil  m  ml 

Use  the  time  sequence  exponential  method  on  (A,B)  and  (C,D) 
to  form 


C 1  C2 

Ui=A  A  .... 

d  d  , 
TT  .  m  .  m- 1 
US=A  A  ... 


c  d,  d9 
.A  mB  lB  Z. 


B 


d,  c  c  , 

1-n  m,.,  m- 1 


A  *B  B  / . .  B 


m 


1 


Take  the  Hamming  vector  of  (Ui  ,  U2  ), 

H=Ui®  U3  =(c,  ®d  )n  (c,®d  ,  )n.  ...  .  (c  ®d.  )n(d  ®c  )n.  .  . 

1  v  1  m  Z  m-1  m  11  m 


(d  ®  c,  )n  . 
m  1 


(5.3) 

This  is  seen  to  be  true  by  inspection  and  definition  of  modulo  Z  sum. 


It  is  to  be  noted  that  each  of  theTheorems  5. 1  ,  5C  2,  5.  3  have 

Hamming  vectors  which  are  anti  symmetric  about  their  center.  This  is 

required  to  satisfy  the  fold- over  method  of  making  the  parity  check. 
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Theorem  5.  3  also  shows  that  along  with  the  anti-s  ymmetrj.c  characteris tic 

the  Hamming  vector  for  time  sequence  exponential  composition  has  unit 

clusters  of  n  zeros  or  n  ones  and  is  also  anti-s ymmet ric  in  both  halves. 

This  is  easily  seen  because  the  first  sum  is  (c  .  ©d  )  while  the  last  sum 

to  operate  on  the  A  code  is  (c^M^),  but  it  is  known  through  the  parity 

check  that  c,©d  ®c  ®d  =  1  .  Therefore  one  of  the  sums  (c,©d  )  or 
1  m  m  1  1  m 

(c  flkL  )  must  be  zero  and  the  other  must  be  one.  To  summarize,  for 
'ml 

the  time  sequence  exponential  form  of  composition  the  first  half  of  the 
Hamming  vector  is  anti -symmetric  about  its  center  (the  quarter  length 

point)  and  the  first  half  and  the  last  half  of  the  vector  are  also  anti- sym¬ 
metric. 

Theorem  5.4  will  concern  the  interlace  exponential  method  of 
forming  composite  codes  and  it  will  be  seen  that  this  method  will  again 
form  a  Hamming  vector  anti-s  ymmet  ric  about  its  center,  but  with  a  dif¬ 
ferent  arrangement  of  the  clusters  of  zeros  and  ones. 

Theorem  5.  4 

The  interlace  exponential  method  of  generating  a  composite  code  has 

the  characteristic  Hamming  vector  H-(c  ©d  )n(d,©c  l^c^ld  ,  )n,  .  .  . 

1  m  1  m  c  m- 1 

(c  ®d,  )n(d  9c,  f. 
m  1  m  1 

Use  the  interlace  exponential  method  on  (A,B)  and  (C  3  D)  to  form 


C1  ^1  c2  ^2 
V1=A  B  A  B 

c  d 

,  .  A  mB  m 

d  c  d  . 

Va=A  mB  mA  m-.1. 

dl  ~cl 

A  B  l. 

Take  the  Hamming  vector  of  (VT,V2) 

H=V!®  V2=  (c.fd  )n(d.§£  )n(c?©d  ©d,)n(d  ®c,  )n  (5.4) 

1  1  m  1  m  2  m-  1  m  I'm  1 

This  is  seen  to  be  true  by  inspection  and  definition  of  modulo  2  sum. 

s|<  # 

The  interlace  exponential  method  again  gives  a  Hamming  vector 

which  is  anti-symmetric  about  its  center  due  to  the  condition  of  satisfying 
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the  fold-over  form  of  the  parity  check,  but  its  clusters  of  zeros  or  ones 

are  now  2n  in  length  rather  than  just  n  as  was  the  case  for  the  time 

sequence  exponential  method  of  generation.  This  characteristic  of 

clusters  of  2n  in  length  is  caused  by  the  necessary  condition  of  the 

parity  check  for  the  (C,D)  code  pair.  For  example,  the  first  two  clusters 

are  (c-.©d  )n  and  (d,©c  )r  and  the  second  two  clusters  are  (c~>®d  ,  )n 

1  m  1  m  2  m-  1 

and  (^2®Cm  1^*  t^le  Parity  check  to  the  i=l  and  i  =  2  bits 

in  the  (C,D)  pair  gives  the  following  two  equations: 

cn®c  ©d ,  ®d  =1  and  c^ffic  1©d:>©d  ,  =1  . 

1  m  1  m  2  m-1  2  m-  1 

Adding  one  modulo  2  to  both  sides  of  both  equations  and  rearranging 

gives  (c  ,  ©d  )©(d,©c  ©1)  =0  and  (c9©d  1  )  ©(d7®c  1®1)  =  0.  It  is 

1  m  1  m  2  m-1  2  m-1  9 

true  in  gene  ral  that  x  ©  1  =  x:  the  re  fore  (c  ,  ©d  )  ©  ( d,  #c  )=0  and 

lml  m 

{c0©d  ,  )  ©  (d0®c  ,  )  =  0  6  The  sums  (c ,  ©  d  )  and  (d,  ®  c  )  must 

2  m-1  2  m-1  1  m  1  m 

be  alike  to  satisfy  this  constraint  with  each  of  the  sums  giving  n  ones 

if  unlike,  or  n  zeros  if  alike,  and  since  the  pairs  are  in  time  sequence 

the  clusters  must  be  in  groups  of  2n.  A  similar  reasoning  will  show 

that  the  parity  check  for  i  =  2  gives  identical  results  for  the  third  and 

fourth  terms  in  the  Hamming  vector  for  the  interlaced  exponential  form 

of  generating  composite  codes*  Therefore  equation  5e  4  can  be  written 

HIV,,  V2)  =  (Cl  f  dm)2n  (c2  ®  dm_i>2n--*- (dm  .  (5.5) 

Before  going  into  the  special  case  applicable  only  to  codes  of  length 

2n  it  might  prove  helpful  to  give  an  example  of  generating  a  code  by  one 

method  and  then  decomposing  it,  to  show  that  it  could  have  been  composed 

by  a  second  method.  The  time  sequences  exponential  scheme  will  be  the 

generating  method  used. 

A=  1 0  C=1001010001 

B= 1 1  D=  10000001  10 

U1=100101  1001  10010101  1011000000000000111100 

U2=01 1010010101010101 100011111 1001100111100 

H=U!®U2  =1111111100110000000011111111001 100000000 
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In  passing  it  should  be  noted  that  the  shortest  clusters  of  zeros  and  ones 
in  H  are  of  length  2  and  that  H  is  anti-  symmetric  about  its  center. 

If  Ui  is  now  time  reversed,  indicated  by  ,  the  complementary 

pair  is  now: 

Ui=001  1  1  10000000000001  101  10101001  1001  101001 
U3=01 1010010101010101 1000111 111001100111100 
H  =  Uj  «U2 =01 01 0101010101010101010101 01 0101 01 010101  =  (01)'° 

This  form  of  the  Hamming  vector  indicates  that  H  is  composed  of  an 
interlace  pair  of  length  20.  Decomposing  both  and  U2  by  assum¬ 

ing  an  interlace  pair  gives 

Uj=01  10000001011  10101  10 
01100000011000101001 
H  =  00000000001  11  11  1111  1  =  o10i10 

u3=  01100000010111010110 
10011  11  110011  10101  10 
h  = 111 1111 1110000000000  =  i10o10 

The  Hamming  vector  of  indicates  that  this  code  of  length  20 

was  composed  by  the  time  sequence  method  from  codes  of  length  10. 

If  the  form  AB  is  assumed  for  the  codes  of  length  10  then 

A=01 10000001 
B=01 110101 10 

This  (A,  B)  pair  is  seen  to  be  the  transformation  ETC^-  EG  on  the 
original  kernel  of  length  10  which  was  used  to  generate  the  code  by  a 
different  method.  Similarly  the  pair  from  U2  is  also  a  transformation 
of  the  original  (C  ,  D)  pair*  Therefore  the  given  code  of  length  40 
could  have  been  generated  by  two  different  methods* 

A  second  example  will  now  be  given  using  the  same  (A,  B)  and 
(C,D)  kernels  as  were  used  in  the  first  example,  but  his  time  utilizing 
the  exponential  interlace  scheme  of  generation  rather  than  the  time  se¬ 
quence  exponential  form. 
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A=  1 0 
B=  1 1 


C=1001010001 

D=  10000001  10 


Vi=101 101000100100001 001 000010001  1101  111000 
V8=01 0010 11 101101 11010001 110100011101111000 
H=V!®V2  =111111111111111100001 1 110000000000000000 

Note  that  the  smallest  cluster  is  four,  as  would  be  expected  of  an  inter¬ 
lace  exponential  scheme  with  the  (A,B)  pair  of  length  two.  The  H 
vector  is  again  anti- symmetric  about  its  center.  Time  reversing  Vi 
give  s 

V, =00011 11011 100010000100100001001000101101 

v7=0100101 1101 10111010001110100011101111000 

H=Vlt®V2  =0101010101010101010101010101010101010101  =  (01)20 

As  before  this  form  of  H  indicates  an  interlace  code  pair.  Decomposing 

both  Vj  and  V-,  into  interlace  code  pairs  gives 

Vj=001 11 1010001000101 10 
01101000010001000011 
H  =01010101010101010101  =  (T  0  1 ) 1 0 
V2=  001  11  1010001000  101  10 
10110111101110111100 
H  =  10101010101010101010  =  (10) 10 

Vj  and  V2  both  break  down  into  interlaced  pairs  of  length  10w  The  top 
row  of  breaks  into  A=01  1 0000001 ,  B  =  0111010110  which  is  the 

transformation  EG  operating  on  the  original  (C,D)  pair  as  in  the 
previous  example.  Each  of  the  codes  for  V2  would  break  down  in 
similar  fashion  to  some  transformation  of  the  original  (C5  D)  code  pair. 
One  use  for  the  method  just  shown  could  be  to  search  for  larger 
kernels  than  those  currently  known.  An  example  might  be  kernels  of 
length  fifty  found  through  the  use  of  time  sequence  exponential  or  inter¬ 
lace  exponential,  where  both  (A;  B)  and  (CjD)  are  codes  of  length  10, 
These  would  form  codes  of  length  200  0  From  this  as  a  starting  point, 

a  decomposition  to  codes  of  length  100  and  then  to  length  50  might 
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be  possible  through  this  technique.  This  technique  will  be  discussed 
much  more  thoroughly  towards  the  end  of  this  chapter,  but  it  is  neces¬ 
sary  to  consider  the  special  case  of  2n  code  composition  and  the 
Hamming  vectors  of  the  sequence  quadruple  form  before  continuing  with 
this  interesting  possibility. 

The  special  generating  method  which  is  applicable  only  to  codes  of 
r 

length  2  is  actually  a  generalization  of  both  the  time  sequence  and  the 

interlace  methods.  As  mentioned  in  the  previous  chapter,  this  method 

takes  a  code  pair  of  length  2  and  from  them  forms  a  code  pair  of 
r  T 1 

length  2  .  The  allowable  combinations  are  formed  by  interlacing 

pieces  of  the  A  code  and  the  B  code  of  length  2m,  whe  re  m= ,  0,1,2,.  .r. 

As  an  example  a.  code  of  ^length  8  can  be  used  to  generate  codes  of 
length  16.  In  this  case  r=3  and  m  has  the  possible  values  0,  1,2,3. 
These  values  of  m  give  section  lengths  of  1,2,  4,  8  respectively. 

A=ala2a3a4a5a6a7a8 

B=blb2b3b4b5b6b7b8 
for  m=0, 

S 1 1 =a  1 b  1 a  2b  2a  3b3a4b4a  5b 5 a  6b  6a  7b  7a  8b  8 
S21=albla2b2a3b3a4b4a5b5a6b6a 7b7a8b8 

n  118 

H=S11®S21  =  01010101010101 01.  =(01)  =  (0  1  ) 

This  is  exactly  the  same  as  the  interlace  procedure  which  is  a  standard 
generating  method* 

The  section  lengths  are  2  bits  long  for  m-1. 

S12=ala2blb2a3a4b3b4a5a6b5b6a7a8b7b8 

S22=ala2blb2a3a4b3B4a  5a  bb5b6a  7a8b7B8 
H=S12»S22  =  00110011001  10011  =  (0212)4, 
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m=2  gives 


S14=ala2a3a4blb2b3b4a5a6a7a8b5b6b7b8 

S24=ala2a3a4blb2b3b4a5a6a7a8b5b6b7b8 

h=s14»s24  =00001  11100001  11  1  =  (04lV, 

For  m=3  the  lengths  are  eight,  which  is  the  entire  code. 

S18=AB 

s28=Ab, 

This  is  the  same  as  the  time  sequence  form  giving  a  Hamming  vector 
H=00000000111  1  11  11  =  (0818). 


The  Hamming  vectors  from  this  special  construction  are  altering 
clusters  of  zeros  and  ones,  each  cluster  being  the  length  of  the  segment 
used  in  the  generation. 

Theorem  5.  5 

r+ 1 

Complementary  sequences  of  length  2  formed  from  an  (A,  B) 
r  . 

pair  of  length  n=2  using  the  special  generating  method  with  segments 

2  n/m 

of  length  m  have  a  characteristic  Hamming  vector  H=(0  1  ) 

Use  the  special  generating  method  utilizing  segments  of  length  m 
of  the  pair  (A,  B)  to  form 

b  1  “■  a  i  a  ^  a  ...  a  b  i  b  ^  e  ®  .  b  a.  _i  a  .  ^  a  b  >  ^  b  >  o  •  •  *  b 

lm  12  3  m  1  2  m  m+1  m+2„  „  2m  m  +  1  m+2  n 

S2m=ala2a3*  *  *  am  blb2“  *  "  bm  am+la  m+2*  *  *  a2m  bm+lbm+2*  *  *  * 


m 


Take  the  Hamming  vector  of  (S  , 

H=S ,  VS,  =  (0mlm)2n/m. 
lm  2m 

This  is  seen  to  be  true  by  inspection  and  definition  of  modulo  2  addition. 

;|c  *  >|c 


It  now  seems  suitable  for  completeness  to  include  a  theorem  which 
has  been  alluded  to  many  times  but  not  proven.  This  is  the  anti- sym¬ 
metric  property  of  all  Hamming  vectors  of  complementary  sequences. 
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The  Hamming  vector  of  a  complementary  sequence  pair  has  a  pat¬ 
tern  which  is  anti- symmetric  about  the  center  of  the  Hamming  vector. 


1. 

A=a1a-)a0ay1  .........  a  na 

1  2  3  4  n-ln 

B  blb2b3b4  . . .  bn- lbn 

is  a  complementary  pair  of 
length  n 

2. 

H(A,B)=  h1h2h3h4.  . .  •  «  hn 

by  definition  of  Hamming  vector 

3. 

a.  ®  a  .  .®  b.  <B  b  ,,  .  =  1 

i  n+l-i  i  n+l-i 

by  equation  2C  4 

4. 

a .  ®  b .  =  h 

1  1  1  an  d 

both  by  definition  of  h 

5. 

a  i-i  .  <8  b  ,  T  .  =  h  ,  n 
n+l-i  n+I-i  n+l-i 

1 

6. 

h.  ®  h  .  ,  =1 

l  n+ 1 - i 

7. 

Add  (hn_j_i  ^  ®  1 )  to  both  sides  of  the  equation 

8. 

h  , .  .  =  h. 

n+ I - l  l  s 

Make  a  change  of  reference  to  the  center  of  the  vector,  which  lies 
between  m/Z  and  m/Z  +  1,  rather  than  the  end  by  having  i=n/Z  -r  then 
n+l-i=n/2+l  +  r 

9.  h  =  h  /9  +l+r  for  all  0  <  r  <n/2  -1  by  step  8, 

n/Z-rn/Z  ~  ~  j  r 

1  0.  Therefore  the  Hamming  vector  is  anti- symmetric  about  its 
center  when  it  is  formed  from  a  complementary  pair. 

*  *  * 

The  Hamming  vectors  of  (1,11)  and  (III,  IV)  in  the  sequence  quad¬ 
ruple  form  are  invariant  under  transformation  if  HfI,II)  and  Hf  III,  IV) 
are  taken  as  an  unordered  pair.  The  following  four  theorems  which 
parallel  Theorems  4.  2 , 4.  3 , 4.  4 , 4.  5  on  Hamming  di  stance  s  will  dem¬ 
onstrate  this  characteristic  of  unordered  Hamming  vector  pairs. 

Theorem  5.  7 

A  change  in  the  order  of  (I,  II)  or  (III,  IV)  or  both  will  not  change 


H(I,  II)  or  H(III,  IV). 
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1.  H(U,  V)=u,  ®  by  definition. 

2.  u.  ®  v.  =  v.  §  u.  since  modulo  2  addition  is  commutative, 

1111 

3.  therefore  H(U,V)  =  H(V,U) 

sje 

Theore  m  5.8 

Complementing  the  pair  (I,  II)  or  (111,1V)  or  both  will  not  change 
H(I,  II)  or  H(III,IV). 

1.  H(U,V)  =  u  $  by  definition, 

2.  u.  ®  v.  =  u.  ®  v.  by  Theorem  2.  1, 

1111  J 

3.  therefore  H(U,V)  =  H(U,V), 

#  Jfc 

In  the  following  proof  it  is  assumed  that  each  vector  I,  II,  III,  IV  is 
of  length  m.  It  is  also  necessary  to  show  that  x  =  x.  This  is  easily 
seen  sincex®  l  =  x  and  x  ®  l=x,  therefore  x  =  x  $  (1  ®  1 )  =  x. 

Theore  m  5.9 

Complementing  one  of  the  pair  (1,11)  and  one  of  the  pair  (III,  IV) 

exchanges  the  Hamming  vectors  H(I,II)  and  H(IIISIV)# 

I=a,  a.a, . a  , 

13  5  n-  1 

II=a  a  o  .  •  •  .  *  a~> 
n  n-  2  2 

III=b1b_bc;.._...b  . 

13  5  m-  1 

IV=b  b  -2 . b7 

mm  2 

1.  H(I,  II)=h.h7h0h>1  h 

1  2  3  4 .  m 

2.  H(III,  IV)-h ^h^h^h^.  .  .  .  .  v  .  0  h^  because  of  the  pa rity  c heck. 

3.  H(1,II)=H(I,  Uj  =h1h2h3h4. 

4.  H(III,  IV)=H(III,IV)  =  R1E2fi3R4R5.  .  .  Rm  =h1hzh3h4.  . .  .  hm  . 

5.  Therefore  H(I,  II)=H(III,  IV)  =H(IIIjIV), 

6.  and  H(III,  IV)  =  H(I ,  II)  =H(I,nT. 

«  *  * 
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Theorem  5„  10 


Exchanging  the  (1,11)  pair  with  the  (IU,IV)  pair  exchanges  their 
Hamming  vectors. 

1.  A  change  in  the  Hamming  vector  can  only  be  caused  by  a  bit 
change, 

2.  Exchanging  the  (I,  II)  pair  with  the  (III,  IV)  pair  changes 
position  but  not  bitsc 


* 


Earlier  in  this  chapter  codes  of  length  40  generated  from  the  quad 
were  decomposed  into  kernels  of  length  10  by  use  of  the  Hamming 
vector.  The  conjecture  was  made  at  that  point  that  perhaps  codes  of 
length  200  manufactured  from  kernels  of  length  10  could  be  decomposed 
into  kernels  of  length  50.  In  order  to  prove  this  decomposition  is  not 
possible  several  more  theorems  are  necessary.  These  theorems  will 
be  proved  before  proceding  to  the  decomposition  problem. 

Theo  re  m  5.11 

Under  all  possible  general  transformations  of  a  complementary  pair 
where  n  >  2  there  are  only  four  possible  Hamming  vectors,  two  of 
which  are  complements  of  the  other  twoe 

1.  (A,B)  are  a  complementary  pair  of  length  n  where  n  >  2, 

2.  H(A,  B)  =  H  interlaced  with  H  ,  where  H  are  the  odd  bits 

o  e  o 

of  the  Hamming  vector  and  H  are  the  even  bits  of  the  Hamming  vector 


in  reverse  order.  The  symbol  H~H^  *  H  will  be  used  for  H^  inter¬ 
laced  with  H  in  exactly  the  same  manner  that  (1,11)  and  (III,  IV)  are 
e 

interlaced  to  form  A,  Bc  Also  h  =  a.^®b^  as  before, 

3.  Let  H1=(l®lll)=a1©b1,  a3§b3.  .  .  aR_  j©bn_  1  =  h^h^  .  .  hn_  j , 
H2=(U»lV)=an#bn,an_2»bn2...a2»b2  =  hnhn_  2. . .  h2, 

H3=»»  IV)=ai«lbn,a3®bn_2 . an_1*b2=KIK3H5...tln_l, 

H4=(Il®UI)=an«bran  2«b3. ....  a2®bn_]=Mn«ri  2. . .  «... 
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by  Theorem  2.  1, 


4.  u,©  v  =  u.  ffi  v. 

ii  ii 

therefore 

(BTT!)=h1  (TT«iv)=h2  (1©iv)=h3  (U®m)=H4. 

5.  If  h.=u.®v. ,  then  h.®  1=  h  =  u.®v.  ©  1  =  u.®  v.  =  u.©  v- 

iii  i  ill  1111 

the  refo  re 

H(U.V)  =  H(U  ,  V)  =(U  ,  V) 

and  Hj  =(T®III)  =  ( I©TTT)  H3={T©IV)  =  (I«lv) 

h2=(TT©iv)  =  (ii«Tv)  H4=(mni)  =  (ii©m). 

6.  H,  *  H  <,  =  h ,  h~  h  -  e  .  .  .  .  h  . 

1  u  1  b  <3  n 

H-  *  H,  =  h  h  ,  . hlf 

2  1  n  n-  1  1 

7.  But  h  =h  .  t  .  by  Theorem  5.6. 

n/2-r  n  /  2  +  1  +  r  J 

8.  H,  *  Hn  =  h ,  h  ~  . . h  /nh  /  ^  e  .  .  h  t 

1  2  12  n/2  n/2  2  1 

H2  *  H1  =  hi  h2 - «-”Rn/2  hn/2”-*h2hr 

9.  Therefore  H]  *  H2  -  H  ,~~H^  =  H,  *  H  . 

10.  H.  *  H.  =  K.K,  .....  U  . 

3  4  I  2  n 

11.  H.*H,=U  ,  _ K,. 

4  3  n  n-  1  1 

12.  ii  /7  =R  j?ti . 

n/2-  r  n/2  +  l+r* 

13.  Therefore  HQ  *  H  .  =  =  H.  *  fiL. 

3  4  4  3  4  3 

14.  Table  5*  1  is  an  exhaustive  list  of  transformations  on  a 
complementary  pair  in  the  sequence  quadruple  form  with  the  Hamming 
vector  for  each  transformation, 

H ^  *  H2>  H  *  G1  ,  H3*H4>  H^H^  and  their  complements  are  the 
only  possibilities  which  exist*  Therefore  only  four  Hamming  vectors 
are  possible  and  by  steps  9  and  13  they  occur  in  pairs  which  are 
complements  of  one  another. 

#  #  # 
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The  E 
of  (I 

would  take  place  in  the 
Hamming  vectors 


HAMMING  VECTOR 
H(I)=H]*H2 

h(c1)=h1*h2 
h(c2=Hi*h2 
H(C)  =H1*H2 
H{A1)=H1*H2 

h(a2)=h^h2 
h(t1)=h4*h3 
H(T2)=H3*H4 
H(T)  =H?«H1 
H(Z )  =H4*H3 
H(Y)  =H4-H3 
H(X)  =H4*H3 
H(W)  =H4-H3 
H(  V)  =H3*H4 
H(U)  =H3^'H4 
H(S)  =H3^H4 
H(R)  =H3*H4 
H(Q)  =H?*H1 
H(P)  =H2^H1 
H(O)  =H4*H3 
H(N)  =H4^H3 
H(M)-  H3*H4 
H(L)  =H4*H3 
H(K)  =H3«H4 
H(J)  =  H2*H. 

H(H)  =H3*H4 
H(G)  =H?*H1 

h(f)  =Hj*h2 

H(D) 

H(it)  =H  *H 


COMMON  VECTOR 
H(I) 

H(I) 

H(I) 

H(I) 

H(I) 

H(I) 

H(T  2 ) 

H(T  ^ ) 

H(I) 

HtTj) 

H(T  x  ) 

h(tx) 

H(T  j  ) 

HIT^ 

H(T  ^ ) 

H(Tl) 

h(t  x ) 

H(I) 

H(I) 

h(tx) 

h(tx) 

H(T1) 

H(TX) 

H{TX) 

H(I) 

H(T  x ) 

H{  I) 

H(I) 

H(I) 

H(  I) 

H(  I) 

H(I) 


H(B)  =H?;!'H1 
H(O)  =H2*H1 

transformation  just  permutes  the  first  two  and  last  two  pairs 
II  III  IV)  and  since  modulo  2  addition  is  commutative  no  change 

H  vector.  The  last  column  identifies  all 
in  terms  of  two  vectors  and  their  complements. 

TABLE  5.1 
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Theore  m  5.12 


If  a  complementary  pair  (A,  B)  of  length  n  is  operated  upon  by  a 
complementary  pair  (C  ,  D)  of  length  m  to  form  codes  U1?  U2  by  the 
time  sequence  exponential  method,  or  to  form  codes  Vx  ,  V2  by  the  inter¬ 
lace  exponential  method  and  either  Ui  or  U2  is  time  reversed,  symbo¬ 
lized  by  U  or  if  Vi  or  V2  is  time  reversed  symbolized  by  V  then 

H(U)  =  H(V)  =  I-^A1  ,  B1  )"m  where  the  (A1  ,  B1)  pair  indicates  some 
general  transformation  of  the  original  (A,  B)  pair. 

1.  A  ®  B  =  H(M)  .  where  B  indicates  the  time  reverse  of 

the  original  B 

2.  A  «  B  =  H(H). 


3. 

B  0  A  =  A  ® 

B  =H(TX). 

4. 

B  «  A  =  A  • 

B  «H(L). 

5. 

Referring  to  Table  5,  1, 

H(M)=H(T1), 

H(H)=H(T1) , 

H(L)  =  H(T1). 

6. 

C1  c2 

U  i=A  A  .... 

c  d, 

...  A  m  B  ’.  ... 

.  Bdm 

d  d  . 

u2=a  mA 

d ,  c 

,  .  .  A  1  B  .... 

B^  1 

and 

VX=A  *  B 


C1  dl  .  C2 


c  d 
A  m  B  m 


d  c  d  ,  d,  c, 
V3=A  m  B  mA  “:?.A  1  B 


Time  reversing  one  of  the  pair  in  each,  for  example,  the  use 
of  T2  yields 


and 


C1  C2 

Ui  =  A  A 

c  1  c2 
U3  =  B  B 


C1  dl  c2 
Vi  =  A  BA. 

8  1  dl  '2 
V2  =  B  A  B  . 


c  d, 
A_mB  1 
c  d , 
B  mA  l. 


B 


m 


m 


c  d 
A  mB  m 


B  mA  m. 
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7. 


c.  c. 

H(U)  and  H(V)  consist  of  clusters  of  (A  ©  B  ±)  and 
d.  d. 

clusters  of  (B  1  #  A  1). 

8.  The  Hamming  vectors  for  these  clusters  have  the  following 
four  possible  forms: 

A©  B  ,  A  ©  B  ,  B  ©  A,  and  B  ©  A. 

These  four  possible  vectors  were  seen  in  steps  1,2,  3,4,  5  to  all  be 

equal  to  H(T1).  If  Ux  and  V1  had  been  time  reversed,  using  the 

operation  T1?  the  result  would  have  been  H  (Tx), 

1  0.  Since  all  the  clusters  are  identical  and  there  is  a  cluster  for 

each  of  the  m  bits  in  C  and  each  of  the  m  bits  in  D, 

H(U)  =  HV)  =  [H( T]_)  ]  for  U2  or  V2  time  reversed  or 

in  general  s 

2m  1  x 

H(U)  =  H(V)  =  H(A1,B1)  where  ( A  ,  B  )  signifies  some  transfor¬ 

mation  of  (A  ,  B). 

>;<  >\< 

Theorem  5e  1  3 

The  kernel  of  length  2  when  operated  on  by  a  complementary  pair 
by  either  the  time  sequence  exponential  or  interlace  exponential  method 
can  always  be  decomposed  at  least  once  by  the  interlace  method. 

1.  If  one  of  the  code  pair  formed  by  either  exponential  method 

is  time  reversed,  H(U)  =  H(V)  =  H(A1,B1)  by  Theorem  5.  14. 

2.  There  are  only  two  possible  Hamming  vectors  for  the 
quad,  01  and  10. 

3.  Therefore  H(U-)  =  H(V-  )  =  (01)^m  or  (10)^m  both  of  which 
are  decomposable  by  the  interlace  method. 

*  m 
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Theo  re  m  5,14 


A  complementary  pair  (Ui,  U2)  formed  by  the  time  sequence 
exponential  method  cannot  be  decomposed  by  the  time  sequence  method. 

1.  Let  (Ux  ,  U2  )  be  a  complementary  pair  of  length  Zmn  formed 
from  the  complementary  pair  (A,  B)  of  length  n  operated 
upon  exponentially  by  the  complementary  pair  (C,D)  of 
length  m,  formed  in  the  same  manner  as  in  Theorem  5.  3. 
H(U)=(c  t  9  dm)n(c2»dm_1)n.  .  .(cm®d1)n(d1®cm)n...(dm®ai)n. 
The  first  mn  bits  in  H(U)  are  m  clusters  composed  of  n 


2. 

3, 


zeros  or  n  ones  dependent  upon  the  bits  of  H(C,D), 

,  -r-v  f-f-i ,  ci  tt  r. mn ,  mn  .  mn  ~mn  r  . 

4.  By  theorem  5.  1  H=U  1  or  L  U  ,  tor  decomposition 

of  a  code  of  length.  Zmn  by  the  time  sequence  method. 

5.  This  requires  the  distance  of  (C,D)  be  zero  or  m,  but  by 
Theorem  4.1  D(C,D)  =  m/2,  therefore  H(U)  is  not  decom¬ 
posable  by  the  time  sequence  method. 

6.  If  Ui  or  U2  is  time  reversed  H(U)  =  H(A1,B1)2m  by 
Theorem  5.12. 

_  .  .  .  .  r  T  t  in*  mn  _  mn  .  .  .  . 

7.  This  is  incompatible  with  the  form  H=U  1  since  it  would 

require  the  first  m/2  clusters  of  H(A1,B1)  to  have  a  distance 
of  zero  and  the  second  m/2  clusters  of  H^A1,  B1  )  to  have  a 
distance  of  n.  Both  are  impossible  by  Theorem  4.  1  which 
requires  a  distance  of  n/2. 


Theo  rem  5.15 

A  complementary  code  pair  (V^Vs)  formed  by  the  interlace  ex¬ 
potential  method  (using  the  same  definitions  and  symbols  as  Theorem  5.  4) 

can  be  decomposed  by  the  time  sequence  method  only  if  H(C,D)  = 

nm/2  m/2  ,m/2  m/2. 

U  1  or  1  0 
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1.  H(V)  =  (c  ©d  )Zn(c,©d  .  }2n _ (c  «d,)2n  by  Theorem  5.  4 

'  '  l  m  Z  m-  1  'ml' 

and  equation  5,  5. 

Z.  The  clusters  of  zeros  and  ones  in  this  Hamming  vector  are  in 

lengths  of  Zn  and  are  dependent  in  value  upon  the  bits  of  H(C,D). 

3.  H=Omnlmn  is  the  requirement  for  decomposition  of  this  code 
by  the  time  sequence  method. 

4.  If  H(C,D)  =  Om^  2  or  complement,  it  will  satisfy  step  3 

/nZnvm/Z  ,,  Znvm/Z  ^mn,  mn  ,T/, TK  ,  r- 

since  (0  )  (1  )  =0  1  =H(V),  and  it  also  satisfies 

Theorem  4,  1 . 

5.  If  Vi  or  V3  is  time  reversed  H(V )  =H(A1 ,  B1  m  by 
Theorem  5.1  Z. 

6.  This  is  incompatible  with  the  form  H=0mnlmn  since  it  would 
require  the  first  m/Z  clusters  of  H(AX,  B1  )  to  have  a  distance 
of  zero,  and  the  second  m/Z  clusters  to  have  a  distance  of  n. 
Both  are  impossible  by  Theorem  4.  1  since  it  requires  a 
distanc e  of  n /Z. 

7.  Step  4  showed  that  H(C,D)  =  0m/  Jlm  .  This  satisfies  the  re¬ 
quirement  for  the  time  sequence  decomposition  by  inspection. 

Any  other  arrangement  of  H(CPD)  would  fail  to  allow  decom¬ 
position, 

*  >:<  sj< 


Theo  rem  5 .  1  6 

Any  complementary  sequence  pair  (Vi,  Vs)  formed  by  the  interlace 
exponential  method  as  in  Theorem  5,4  (using  the  same  definitions  and 
symbols),  the  quad  as  the  (C ,  D)  pair,  is  time  sequence  decomposable. 

1.  The  quad  has  only  two  possible  Hamming  vectors,  01  and  10. 

Z.  This  satisfies  the  form  H=0m  where  m=Z,  the  quad 

length. 

3,  Vi  and  V2  are  therefore  time  sequence  decomposable  by 
Theorem  5e  1  5. 


* 


* 


* 
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Theorem  5.17 


A  complementary  pair  formed  by  either  exponential  method  is  de¬ 
n/2 

composable  by  the  interlace  method  if  and  only  if  H(A"  ,  B1)  =  ( 0 ]  )  or 

n  /  2 

(10)  .  All  definitions  and  symbols  are  assumed  to  be  the  same  as 

those  in  Theorems  5*  3,  5.  4  and  5.  12. 

1.  By  Theorem  5.  2  these  codes  are  decomposable  by  the  inter¬ 
lace  method  only  if  they  are  of  the  form  H=(01)nm  or  (10)nm. 

2.  H(U)  and  H(V)  are  clusters  of  at  least  n  ones  or  n  zeros 
and  therefore  cannot  have  the  form  (01)nm„ 

3.  H(U)  and  H(V)  =  H(A1,B1)'  m  by  Theorem  5.12. 

4.  If  and  only  if  H(AX,  B1)  has  the  form  (01  )n  or  (10)n^  can 

H(U)  or  H(V)=  [(0.1)n/?']  ?m=(01)  nm,  or  using  10  instead 

of  01,  (10)nm.  This  satisfies  the  requirement  of  Theorem 

r  r 

5.  2  that  the  Hamming  vector  be  of  the  form  (01)  or  (10)  . 

# 

A  brief  description  of  the  method  used  to  prove  that  codes  of  length 
200  formed  from  kernels  of  length  10,  cannot  be  decomposed  into  ker¬ 
nels  of  length  50,  will  now  be  given  as  a  guide  to  the  actual  proof  of  the 
theorem.  The  method  used  is  to  exhaust  all  possible  methods  of  de¬ 
composition  by  examining  all  possible  Hamming  vectors  for  the  (AfB) 
and  (C,D)  pairs  Table  5.2  lists  all  possible  Hamming  vectors  for 
kernels  of  length  10c 

Kernel  1  Kernel  2 

A= 10000001 10  A= 1001 1011 11 

B= 1001 010001  B=1 100001010 

The  4  possible  Hamming  vectors.  The  4  possible  Hamming  vectors. 

0001010111  0101100101 

1110101000  101001 1010 

0000101  1  11  1100101100 

1111010000  001 1010011 

All  possible  Hamming  vectors  of  kernels  of  length  10o 

TABLE  5.2 
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The  two  methods  of  generating  codes  of  length  ZOO  from  kernels 
of  length  10  are  the  time  sequence  exponential  method  and  the  interlace 


exponential  method  given  respectively  by 


C1  .  c2 


.A  1  °B  1 . Bd10 


1 


U1=A  A 

T1  A  d  1  0  A  .  d  1  T3  C  1  0  ' 

U2=A  A  ...A  B  . B 

V1=ACl  B^1  A"2  B^ . AC'°Bdl0 

V2=Adl°B"10Ad9  B"9 . Adl0BJl 


(5.6) 


(5.  7) 


where  (A,B)  and  (C,D)  are  kernels  of  length  10. 

The  Hamming  vectors  of  these  codes  of  length  ZOO  are  anti- sym¬ 
metric  about  their  centers  and  are  in  clusters  of  10  zeros  or  10  ones 
for  H(U),  and  in  clusters  of  Z0  zeros  or  Z0  ones  for  the  H(V)  vector. 
Time  reversing  either  Ui  of  U2  and  similarly  either  V1  or  V2  » 
gives  Z0  clusters  of  Hamming  vectors  of  the  forms  listed  in  Ihble  5.  2, 
since  each  of  the  20  A  codes  is  now  matched  with  some  transform  of  B. 

The  known  possible  ways  for  kernels  of  length  50  to  generate  com¬ 
posite  codes  of  length  20  are: 

1.  Either  interlace  or  time  sequence  kernels  of  length  50  to 
form  composite  codes  of  length  100  and  then  either  inter¬ 
lace  or  time  sequence  these  n=100  codes  to  form  composite 
codes  of  length  ZOO. 

Take  kernels  of  length  50  and  either  time  sequence  expo¬ 
nential  or  interlace  exponential  with  the  quad  to  form  codes 
of  length  ZOO. 


2. 


3.  Take  kernels  of  length  Z  and  either  time  sequence  exponen¬ 
tial  or  interlace  exponential  with  kernels  of  length  50. 

An  examination  of  all  possible  Hamming  vectors  of  length  ZOO  from 

each  of  these  methods  will  show  the  incompatability  of  these  Hamming 

vectors  with  those  generated  by  the  kernels  of  length  10. 

Theorem  5.18 

No  codes  of  length  ZOO  formed  from  kernels  of  length  10  can  be 

decomposed  into  kernels  of  length  50  by  standard  methods  of  decomposi¬ 
tion.  80 


Given  U1  U2  and  Vi  V2  both  complementary  pairs  of  length 
formed  by  the  time  sequence  exponential  and  interlace  expo¬ 
nential  methods  as  given  in  equations  5.  6  and  5.  7  respectively. 
The  time  sequence  method  of  decomposition  is  impossible  for 
each  of  the  following  Hamming  vectors  for  the  following  reasons. 

a.  H(U)  can  never  be  decomposed  by  time  sequence. 

Theorem  5.  14. 

5  5 

b.  H(V)  requires  a  form  0  1  by  Theorem  5.  15. 

Table  5.2  contains  no  vector  of  this  form. 

c.  H(V  )  and  H(U)  can  never  be  decomposed  by  time 
sequence.  Theorems  5.14,15. 

The  interlace  sequence  method  of  decomposition  is  impossible 

for  each  of  the  vector  forms  for  the  following  reasons. 

a.  H(U)  and  H(V)  can  never  be  decomposed  by  interlace. 
Theo  rem  5.  1  7. 

b.  H(V)  and  H(U  )  are  decomposable  by  interlace  only  if 
H(A1,B^‘)is  of  the  form  (01)  or  (10) ^  by  Theorem  5.17. 
Table  5.  2  contains  no  such  vectors. 

Assume  that  one  of  the  Hamming  vectors  H(U)  or  H(V)  has  a 

second  exponential  form,  first  with  the  kernel  of  length  2 

operated  on  exponentially  by  the  kernel  of  length  50,  and 

secondly  by  the  kernel  of  length  50  operated  on  exponentially 

by  the  quad.  L.et  H1  be  the  designator  for  the  Hamming 

vectors  of  the  2  and  50  combinations. 

a.  Consider  first  the  quad  acted  on  exponentially  by  the 

kernel  of  length  50.  H^U)=H"(V)=  H1(A1,BX)<I  But  since 

A^B1  is  the  quad,  by  Theorem  5.  13,  it  is  decomposable 
by  the  interlace  scheme.  This  is  not  possible  since 
step  3  above  exhausted  all  possible  interlace  decompo¬ 
sitions. 

b.  Next  consider  the  kernel  of  length  50  acted  on  exponen¬ 
tially  by  the  quad.  Considering  first  the  exponential 
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interlace  scheme,  by  Theorem  5.  15  if  Hx  (C  ,  D)=01  or  10, 
which  it  does  for  the  quad.  Therefore  the  code  is  time  se¬ 
quence  decomposable.  However,  in  step  2  all  of  these  time 
sequence  possibilities  were  exhau sted  and  therefo re  this  de¬ 
composition  is  not  possible,  H1(U)  would  have  a  form  con¬ 
taining  either  50  zeros  or  50  ones  in  clusters,,  For  the 

original  10  length  codes  to  form  these  would  require  a 
5  5 

0  l  Hamming  vector  which  is  not  listed  in  Table  5.  2. 

This  exhausts  all  possible  Hamming  vectors  that  might  lead  to  de¬ 
composition,  therefore  codes  of  length  200  formed  from  kernels  of 
length  10  cannot  be  decomposed  into  kernels  of  length  50  by  standard 
methods. 

$  ❖  * 

Table  5.3  shows  the  possible  Hamming  vectors  for  codes  of 
length  26. 

A=01 00 1101 111010111100111010 

B= 10 1100100001 11111100111010 

The  4  possible  Hamming  vectors. 

Ill  1  111  1  11  1101000000000000 

000000000000101 11111111111 

1 1 101 11011 1010100010001000 

0001 000100010101 1101110111 

All  possible  Hamming  vectors  of  the  kernel 
of  length  26. 

Table  5e  3 

An  examination  of  Table  5.  3  which  is  exhaustive  for  the  Hamming 
vectors  of  length  26  and  Table  5.2  for  10,  shows  that  the  next  two 
theorems  can  be  proved  by  identical  methods  to  those  used  in  Theorem 
5.  18.  Since  these  proofs  exactly  parallel  that  of  Theorem  5.  18  they 
will  not  be  given. 

Theorem  5e  1  9 

Complementary  pairs  of  length  520  formed  from  the  time  sequence 


82 


exponential  or  interlace  exponential  from  kernels  of  length  26  and  10 
cannot  be  decomposed  by  standard  methods  into  kernels  of  length  130. 

*  *  * 


Theo  rem  5.20 

Complementary  pairs  of  length  1352  formed  from  the  time  sequence 
exponential  or  interlace  exponential  method  from  kernels  of  length  26 
cannot  be  decomposed  by  standard  methods  into  kernels  of  length  338. 

^  J1,: 

The  study  of  Hamming  vectors  thus  far  has  disclosed  that  a  pattern 

r  /ril  .n/2  _  n / 2  ,n/2  ,  _  .  .  _  ,  . 

ol  (01)  or  U  1  or  their  complements  in  either  the  (A,B) 

pair  or  (C,D)  pair,  is  required  in  order  that  a  code  formed  from  one 
kernel  length  by  an  exponential  method  be  decomposed  into  a  pair  from 
a  different  kernel  length.  The  quad  is  only  kernel  which  thus  far  satis¬ 
fies  this  criteria.  However,  the  other  kernels  come  close,  as  a  check 

* 

4  4  2  2 

shows  that  0  101  and  (01)  10(01)  are  among  the  Hamming  vectors  of 

12  12 

length  10,  and  0  101  is  a  Hamming  vector  of  length  26. 

If  the  quad  is  involved  in  the  formation  of  composite  codes  it  makes 
their  decomposition  more  likely.  This  is  shown  in  Theorem  5.  12  and 
Theorem  5.  16.  Appendix  II  demonstrates  this  fact  quite  clearly  since 
it  contains  all  codes  of  length  16,  less  some  operational  redundancies. 
Each  of  the  code  pairs  and  its  time  inverse  is  decomposable  by  some 
method  into  shorter  length  codes.  Appendix  III  lists  all  the  pairs  of 
length  20,  less  some  operational  redundancies.  Each  code  pair  or  its 
time  inverse  can  be  broken  into  its  originating  pair  through  the  use  of 
Hamming  vectors,  but  no  code  and  its  time  inverse  can  both  be  decom¬ 
posed  by  the  standard  Hamming  vector  methods. 


The  writer  would  like  to  end  this  chapter  with  a  conjecture  that 
kernels  of  length  50  do  not  exist.  This  conjecture  is  based  upon  three 
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bits  of  evidence,  the  proof  of  Theorem  50  18  being  the  first.  The  second 
piece  of  evidence  is  based  upon  the  fact  that  50=7'  +1  =5^+5^.  From  this 
one  might  assume  that  codes  of  ZOO  formed  from  lO^s  might  cover  just 
one  set  of  the  ones  determined  by  equation  2e  5.  This  proved  not  to  be 
true  however,  since  the  number  of  ones  in  both  the  exponential  formations 
of  codes  of  length  ZOO  from  kernels  of  length  10,  if  they  were  decompo¬ 
sable  into  kernels  of  length  50,  would  exhaust  all  possible  unordered 
pairs  of  ones  for  n=50  as  determined  by  equation  2fc  5e  The  last  bit  of 
evidence  to  base  this  conjecture  upon  is  that  the  code  length  is  not  twice 
a  prime  number.  The  only  previous  kernel  which  might  have  existed,  but 
was  not  twice  a  prime  number,  was  n=18,  and  it  did  not  exist*  Similar 
reasoning  extends  this  conjecture  to  kernels  of  length  130  and  338. 
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CHAPTER  VI 


SUPPLEMENTARY  AND  CYCLIC  COMPLEMENTARY  CODES 

This  chapter  is  concerned  with  two  classes  of  codes.  The  first  of 

these  classes  consists  of  quadruples  of  sequences  with  the  property  that 

the  total  number  of  likes  at  each  spacing  equals  the  total  number  of  un- 

2 

likes  at  the  same  spacing.  These  are  called  supplementary  codes.  The 
second  class  of  codes,  which  will  be  defined  later,  consists  of  the  cyclic 
complementary  codes.  Both  of  these  code  types  have  complementary 
sequence  pairs  as  a  possible  subset;  this  allows  the  use  of  their  prop¬ 
erties  as  screens  in  the  search  for  new  kernels. 

Supplementary  c ode  s  will  be  di scus  sed  f ir st ,  since  the  n=26  search 
involved  the  use  of  this  property  as  a  screen,  while  the  cyclic  comple¬ 
mentary  property  was  not  used  until  the  search  for  kernels  of  length  34, 
Only  one  theorem  will  be  proved  for  each  of  these  types  of  codes,  since 
these  two  theorems  were  the  only  ones  applied  in  the  actual  search  for 
new  kernels. 

If  U , U-U0 . u 

12  3  n 

v,  v9  v- . v 

12  3  n 

w  .  w~  w0 . w 

12  3  n 

x,x-x0 . x  are  sequences  of  zeros  and  ones  and 

12  3  n 

satisfy  the  constraint  of  equation  6C  1,  they  form  a  quadruple  which  is 
supplementary. 

i=j 

E  (u.©u  .  .  J&l)  +  (v.©v  .  .§1)  +  w  ©w  .  .©1)  +  x.©x  .  .©1)  = 

'  l  n+1-j  '  '  l  n+1-j  l  n+i-j  i  n+i-j 

i=l 

j 

E  (u.flhi  .)  +  v.©v  ,  ,  .)  +  (w.©w  .  .)  +  x.©x  .  .)  =  2j  (6.  1) 

1  n+i-j'  i  n+l-j  x  l  n+i-  j  l  n+i-jf  J  '  ' 

for  all  j,  1  _<  )<  n  -  1.  Note  that  there  are  2j  like  pairs  and  2j  unlike 

pairs  at  a  spacing  of  n-j. 
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For  simplicity  in  the  proof,  we  use  the  symbolism  that  was  used  in 


Chapter  2, 
let 


Ltt  =  £  u.  ®  u  ,  .  .©1 
U  .  ,  1  n+i-j 


J 

£ 

i=  1 

i 
£ 
i=  1 


for  all  j,  1  j<  j  <  n-  1 


UTT  =  £  u.  ©  u 
U  i  n+i-j 


A  restatement  of  equation  6.  1  in  this  symbolism  is 

^ W  +  =  UU  *  UV  UW  *  UX  ^ 

Theorem  6.  1 

Any  complementary  sequence  pair,  (A,B),  written  in  standard 
(I  II  III  IV)  form  has  the  property  that  I,  II,  III  and  IV  are  a 
supplementary  quadruple. 


1. 

2. 


la=ub  and  ua=lb 


by  the  definition  of  complementary. 
Considering  just  the  even  values  of  j  rather  than  all  values, 


3. 

4. 

5. 

6. 


la-li  +  Ln, 

ua=ui  +  uir 

LB=LIII  +  ^IV’ 

UB=UIII  +  UIV 

Li  +  =uiii  +  urv* 

L,m+L,iv=ui  +  uir 

Adding  the  equations  in  steps  3  and  4  gives 

Li  +  Ln  +  +  Liv  =  ui  +  TTn  +  um  +  urv  • 

This  is  the  same  as  equation  6.  2  and  I,  II,  III,  IV  formed 


in  the  sequence  quadruple  form  from  a  complementary  pair 
are  supplementary* 

Since  the  supplementary  property  applies  only  to  the  even  spacings 
in  a  code  pair,  the  converse,  that  all  interlaced  supplementary  quad¬ 
ruples  are  complementary,  is  not  true. 

*  * 
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The  cyclic  complementary  property  is  useful  in  itself  for  a  communi¬ 
cations  or  telemetering  system;  however,  at  this  time  our  concern  is 
for  its  property  as  a  necessary  condition  for  a  code  pair  to  be  complemen¬ 
tary.  A  cyclic  sequence  or  code,  as  the  name  suggests,  is  a  never 
ending  sequence  of  zeros  and  ones  which  has  a  period  of  n  bits.  A  cyclic 
complementary  sequence  pair,  is  a  pair  of  cyclic  codes,  each  of  period 
n,  where  the  number  of  likes  of  one  sequence  equals  the  number  of  un¬ 
likes  of  the  other  sequence  for  n  possible  matches  for  all  spacings 

from  1  to  n-1,  or  stated  in  terms  of  the  bits  of  the  A  and  B  codes 
n  i=n 

C .  =  s  a. a  a  ...  =  E  b.®  .,.®  1  1  <  j  <  n- 1 .  (6.3) 

J  i=i  1  n-  j+i  i=1  i  n-j+i  “  J  “  . 

where  a  =  a  and  b  ,  =  b  since  the  period  is  n. 

n+i  i  n+i  -  i 

Theorem  6.  2 

A  complementary  sequence  pair  (A,  B)  of  length  n,  if  written  in  a 
cyclic  fashion,  is  always  a  cyclic  complementary  pair  of  period  n. 

1.  Since  (A,  B)  are  a  complementary  pair  they  satisfy 
equation  2.3. 

j  i=j 

f.  =  S  a.®  a  ...  =  S  b.®  b  .®  1  for  j,  1= j=n- 1. 

J  i  =  1  i  n-j+i  i=1  i  n-j+i  J 

2.  Expanding  for  f ^  gives 

a,®  a  =  b ,  ffi  b  ©  1  which  is  equal  to  a  ©a,  =  b  ©  b,  ©  1 « 
lnln  nlnl 

3.  Expanding  for  f  ^  gives 

(a,®  a9)  +  (a^©  a0 )  +  a~©  a.)  +.  . . .  (a  n  ©  a  )=(b,  ©b^©l)  + 

'  1  2'  v  2  3  3  4  n-1  n  12' 

(b2©  b3©  1)  +  (b3©  b4©  1)  + _ _  .+(bn-1®  bn©  1). 

4.  Adding  the  equations  of  step  2  to  step  3  gives 

5.  Similarly  expanding  f^  gives 

(a^©  4  (a2®  an)  =  (b  j©  bR_ -^  ©  1)  +  (b2®  b^©  1) 

which  is  equal  to 

(an„i©  aL)  +  (an©  az)  =  (b^©  bj©  1)  +  (bn©  b2  ©  1). 
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6.  Expanding  ?  gives 

(a1«a3)  +  (a2®a4)  +  (a3®a5)  +  .  .  .  (an_2®an)  =  (b^b^  1)  + 

(bz®b4®l)  +  (b3®b5®l)  +....+  <bn  z®bn®l) 

7.  Adding  the  equations  of  step  5  to  step  6  gives 

8.  Continuing  this  same  procedure  through  all  possible  values  of 

j  would  show  that  C.  =  f .  +  f  Therefore  the  cyclic  com- 

J  J  n“J 

plementary  constraint  equations  are  based  on  the  sum  of  two 
restricted  portions  of  the  complementary  constraint  equations. 
Therefore  all  complementary  pairs  have  to  be  cyclically 
complementary.  The  converse  is  not  nece  ssar  y  true  since 
Cj  is  a  sum,  and  in  a  sum  the  addend  and  augend  are  not 
unique.  * 

sjtf  sjfi 


Since  C.  =  f.  +  f  ■,  each  constraint  of  the  complementary  property 
J  J  ^  J 

is  used  twice  as  j  varies  from  1  to  n-1;  therefore  is  symmetric 

and  is  centered  at  C_  /->.  In  the  application  of  the  cyclic  complementary 


n 


/2- 


property  in  Chapter  8,  half  of  the  characteristic  cyclic  number  is  de¬ 
leted  due  to  this  symmetry  property. 

Both  of  the  screens  described  in  this  chapter  were  suggested  by 
14  19 

Dr.  Golay;  *  however,  to  the  best  of  my  knowledge  the  theorems 
and  their  proofs  have  not  appeared  in  print. 
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AN  EXHAUSTIVE  SEARCH  FOR  KERNELS  OF  LENGTH  26 

CHAPTER  VII 


In  his  paper  Golay  emphasized  the  need  for  an  exhaustive  search  for 
complementary  code  pairs  of  length  26.  This  chapter  describes  such  a 
search.  The  purpose  of  this  search  was  to  determine  if  there  were  any 
kernels  of  length  26  and  if  any  of  these  exist,  the  number  of  such 
kernels.  It  would  prove  very  difficult  to  surmise  the  number  of  kernels 
in  a  code  of  this  length  predicated  on  the  known  kernels  of  other  lengths, 
since  n=2  had  one  kernel,  n=10  had  two,  and  n=  1 8  had  zero.  An  ex¬ 
haustive  search  for  kernels  of  length  26  might  throw  some  light  on  a 
possible  general  method  for  finding  kernels  of  longer  lengths,  or  might 
disclose  some  sort  of  pattern  showing  the  distribution  of  kernels  among 
the  possible  code  lengths.  * 

The  only  feasible  method  to  accomplish  this  search  was  with  a  high 
speed  digital  computer.  This  chapter  is  concerned  with  the  computer 

program  and  the  results  of  the  exhaustive  search  for  kernels  of  length  26. 

52 

At  first  glance  it  would  appear  that  there  are  2  possibilities  to 
be  screened  for  an  exhaustive  search  for  all  possible  kernels  of  length 
26.  Although  this  is  true,  the  application  of  some  of  the  theorems  de¬ 
veloped  in  the  earlier  chapters  immediately  e  liminate  s  from  considera- 

52 

tion  large  blocks  of  the  2  possibilities  which  are  either  redundant  or 
impossible.  This  is  a  very  necessary  procedure  since  even  with  the 
highest  speed  computers  of  today  it  would  be  impossible  to  investigate 
this  number  of  possibilities  in  a  life  time. 

The  program  was  developed  for  the  CDC  1604;  this  computer  ac¬ 
complishes  approximately  200,000  operations  per  second.  This  speed 
coupled  with  the  powerful  screens  used  in  programming  made  possible  a 
reduction  of  the  computer  run-time  to  approximately  75  hours.  A  brief 
outline  of  the  various  screens  used  and  their  reduction  factor  will  now  be 
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given  before  continuing  in  more  detail  with  the  actual  programming 
techniques  used.  The  word  code  or  code  pair  as  used  in  this  cha.pter 
and  also  in  the  next  chapter  will  be  understood  as  possible  complemen¬ 
tary  code  pair. 

The  first  screen  utilized  was  the  number  of  ones  which  must  appear 

Z  Z 

in  each  of  the  code  pairs.  This  was  given  by  n=(n-p-q)  +  ( p - q ) 

as  derived  in  chapter  Z,  where  n  is  the  length  of  the  code,  p  is  the 

weight  of  the  A  code  and  q  is  the  weight  of  the  B  code.  The  possible 

solutions  fo  r  n=Z  6  were  the  unordered  weight  pairs  (16,15),  (16,11), 

(11,  10),  (15,  10).  The  pair  (p=l6,  q=l5)  was  arbitrarily  chosen  to  be 

used  in  the  program.  This  reduced  the  total  number  of  possible  code 

pairs  from  Z  ^8x10  down  to  -q-gr-  —  x  ’qgy-j-y r~  ^  4x10  .  This 

was  a  reduction  by  a  factor  of  ZOO,  which  brought  the  life  time  search 

down  to  a  little  less  than  a  year,  but  this  amount  of  time  was  still  not 

feasible  for  a  computer  search.  Golay  remarked  in  a  footnote  that  a 

complementary  pair  might  be  thought  of  as  being  composed  of  two  inter- 

9 

laced  half  length  codes,  in  this  case  n=l3.  Since  it  was  obvious  the 
total  number  of  ones  in  both  of  the  interlaced  pair  must  equal  the  number 
of  ones  in  one  of  the  complementary  code  pair,  equation  Z.  5  was 
applied  again,  this  time  for  n=13.  The  unordered  weight  possibilities 
for  n=13  were  (9,  7),  (9,  6),  (6,4),  (7,4).  In  order  to  be  compatible  with 
the  (16,  15)  pair  for  n=Z6,  the  (9,  7)  pair  was  selected  for  the  code  con¬ 
taining  16  ones  and  the  (9,  6)  pair  was  chosen  for  the  other  code  which 
contained  15  ones. 

The  basic  problem  had  now  been  reduced  to  all  possible  combinations 

of  the  four  thirteen  bit  sequence  families.  This  was  equal  to 

~ x  ^r^jT  x  71  ^7 — -  x  ■  “gr '  ~  715x715x1716x1716  or  approximately 

1  Z 

1.5x10  possibilities.  Although  this  reduced  the  problem  by  one  more 


90 


order  of  magnitude  it  was  still  much  too  large  for  a  practical  search. 

Dr.  Golay  in  a  private  communication  with  the  author  pointed  out  the 

possibility  of  using  the  supplementary  characteristic  as  an  additional 
14 

screen.  Theorem  6.  1  gives  the  necessary  and  sufficient  condition 
that  complementary  codes  sequence  quadruple  form  be  supplementary 
sequence  quadruples*  The  codes  of  length  13  were  therefore  categor¬ 
ized  according  to  their  first  3  and  last  3  bits.  The  number  of  like 
pairs  at  spacing  twelve,  at  spacing  eleven,  and  at  spacing  ten  was  com¬ 
puted  for  each  of  the  13  bit  sequences.  For  example,  0111100111011 
has  0  likes  for  spacing  twelve,  1  like  for  spacing  eleven,  and  3  likes 
for  spacing  ten.  The  block  number  given  to  0111100111011  was  there¬ 
fore  013.  Out  of  sixty  four  possible  combinations  only  eighteen  different 
block  numbers  were  generated  for  each  of  the  twrf  number  sequences  of 
all  possibilities  of  thirteen  bits  with  nine  ones  and  with  seven  ones.  An 
examination  of  Figure  7.3  reveals  that  there  are  a  widely  varying  number 
of  members  within  each  of  the  blocks.  The  supplementary  characteristic 
states  that  the  total  number  of  like  pairs  for  each  spacing  must  equal  the 
total  number  of  unlike  pairs  for  the  same  spacing*  With  four  codes  at  a 
spacing  of  twelve  there  are  4  possibilities,  therefore  two  of  them  must 
be  likes;  at  a  spacing  of  eleven  there  are  eight  possibilities,  therefore 
four  of  them  must  be  likes;  and  likewise  at  a  spacing  of  ten  there  are 
twelve  possibilities  of  which  six  must  be  likes  to  satisfy  the  supplemen¬ 
tary  property.  Block  numbers  were  added  for  each  possible  combination 
to  see  if  the  total  were  246,  and  in  cases  where  the  total  was  246  those 
particular  codes  within  these  blocks  were  then  sent  to  the  next  necessary 
but  not  sufficient  condition  for  the  codes  to  be  complementary.  If  the 
total  was  not  246,  the  blocks  of  codes  were  rejected  The  power  of  this 
screen  is  shown  in  the  example  where  say  two  of  the  blocks  had  code 
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numbers  000,  000*  The  only  possible  combination  to  add  to  Z46 

would  be  123,  123  for  the  other  two  codes*  This  cuts  the  total  number 

of  possibilities  in  this  particular  block  combination  down  from 

106 

1716x1716  =  2*9 x  to  only  64  possibilities.  Although  this  admittedly 

was  the  most  extreme  example,  this  screen  reduced  the  number  of 

9 

possible  codes  checked  to  about  2x10  or  further  reduction  of  three 
orders  of  magnitude. 

The  next  necessary  condition  to  be  checked  was  the  parity  test, 
equation  2.  4.  All  code  quadruples  which  satisfied  the  parity  test  were 
then  sent  to  the  necessary  and  sufficient  like  pair,  unlike  pair  check  for 
each  spacing.  The  like,  unlike  check  subroutine  was  quite  long  and  also 
involved  word  unpacking.  This  subroutine  which  was  ignored  in  the 
calculation  of  run  time  increased  the  computer  search  time  from  a  calcu¬ 
lated  70  hours  to  an  actual  run  time  of  the  order  of  75  hours. 

The  half  length  codes  in  the  program  were  interlaced  in  the  same 
manner  as  was  done  in  the  operations  group  formulation*  This  method 
of  combination  made  the  parity  check  easy  to  calculate  since  the  bit 
positions  in  the  computer  words  were  the  same  as  the  sequence  posi¬ 
tions  of  the  codes*  Modulo  2  addition  was  directly  applied  to  the  code 
quadruple  to  form  a  vector  of  all  ones,  providing  the  parity  test  held. 

If  the  code  vector  was  not  all  ones  the  code  quadruple  was  rejected. 

It  is  to  be  noted  that  by  selecting  the  ordered  (16,15)  pair  for  the 
number  of  ones  in  the  A  and  B  code,  the  operations  Ci,  C2  »  and  C 
were  eliminated  from  possible  redundant  c  onsideration*  These  were 
eliminated  by  their  requirement  for  (10,15),  (16,11),  and  (10,11) 
ordered  pairs  of  ones  respectively,  based  on  equation  2.  5,  The  inter¬ 
lacing  of  the  half  length  codes  with  ordered  pairs  of  (9>  7)  ones  and 
(9,  6)  ones  respectively  eliminates  from  consideration  the  operations 
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Flow  Chart  for  Computer  Search 
For  Kernels  of  Length  26. 

FIGURE  70  1 
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Ti  ,  T3  ,  T  and  Ai  since  these  would  have  required  (7,9,)  (9,6), 
(9>7)(6,9)>  (7,9)(6,9)  and  (4,7)(4,6)  ones  respectively  The  only  basic 
transformation  not  deleted  by  controlling  the  weight  m  each  of  the  code 
quadruples  was  As.  This  is  seen  to  be  true  because  if  the  pair 
(9>7)(9,6)  has  the  operation  A3  performed  on  it,  the  result  is  (9,6)(9>7). 
If  this  code  is  now  transformed  by  the  trivial  exchange  operation,  E,  a 
code  pair  of  the  original  form  is  generated.  The  scheme  which  was  used 
to  eliminate  the  As  r edundanc  y  will  be  discussed  late  in  this  chapter. 

Figure  7.  1  is  a  rough  flow  diagram  for  the  program  used  in  the 
exhaustive  search  for  kernels  of  length  26.  With  modifications  in  data 
as  necessary,  this  program  was  also  used  for  exhaustive  searches  of 
composite  codes  of  16  and  20, 

A  brief  description  of  the  blocks  in  Figure  70  1  will  give  the  reader 
some  familiarity  with  the  programming  philosophy  used  in  this  search. 

The  actual  programs  as  written  for  the  CDC  1604  are  given  in  Appendix  V. 

The  number  generation  subroutine  is  able  to  generate  any  number  up 
to  48  bits  in  length  (word  size  of  the  1604)  with  any  number  of  ones  up 
to  a  maximum  of  48.  This  is  accomplished  in  the  1604  by  inserting  a 
number  of  the  form  00000.  .  .  00001  111. 111  into  the  accumulator, 
where  the  total  length  of  the  number  is  48  bits  and  it  contains  p  ones 
and  has  a  code  length  of  n  This  number  is  then  shifted  left  48-n  bits. 
This  shifted  result  is  then  checked  to  see  if  it  is  negative,  which  indicates 
a  one  in  the  most  significant  bit  position.  A  procedure  is  now  set  up  to 
count  a  one  or  zero  and  then  shift  left  one  bit  a,nd  repeat  the  process  for 
the  entire  n  bits.  If  p  ones  are  counted  the  number  is  stored,  if  not  it 
is  rejected.  A  number  one  larger  than  its  predecessor  is  inserted  into 
the  system  and  the  process  is  repeated  until  the  ones  are  all  in  the 
uppermost  bit  positionsa 
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Blocks  2,3,4  are  used  to  compute  a  supplementary  characteristic 
number  for  each  of  the  715  thirteen  bit  numbers  with  9  ones  and  also 
for  each  of  the  1716  thirteen  bit  code  numbers  with  7  ones;  then  in 
accordance  with  their  supplementary  numbers  these  are  stored  into 
group  blocks*  The  supplementary  number  was  derived  by  using  a  mask 
to  expose  for  consideration  only  the  first  3  and  last  3  bits  of  each  of 
the  13  bit  code  numbers.  These  masked  numbers  were  then  compared 
with  an  exhaustive  list  of  all  possible  combinations  of  zeros  and  ones  in 
these  6  bits.  These  64  possible  combinations  were  in  an  M  to  one 
correspondence  with  a  list  of  18  possible  supplementary  characteristic 
numbers  (see  Figure  7.2).  This  correspondence  allowed  the  temporary 
attachment  of  this  supplementary  number  to  each  of  the  codes.  The  code 
numbers  were  rearranged  into  blocks  according  to  supplementary  number 
and  at  the  same  time  had  their  temporarily  attached  supplementary 
number  deleted.  This  part  of  the  program  is  purely  for  putting  the  data 
into  a  useable  state  and  is  used  only  once  per  run0  The  rest  of  the  pro¬ 
gram  is  highly  iterative. 

Consideration  of  the  (I  II  III  IV)  form  of  the  possible  26  bit 
codes,  shows  that  each  of  the  18  blocks  of  I,  IIS  III,  and  IV  must  be 
cornpared  against  each  other  indicating  1  8x1  8x1  8x1  8  1  05,  000  possible 

block  comparisons.  This  was  cut  down  somewhat  by  use  of  a  slightly 
different  technique.  This  was  accomplished  by  loading  the  computer 
accumulator  with  246  and  from  this  subtracting  the  first  block  number 
of  I,  then  subtracting  the  first  block  number  of  III,  and  then  subtrac¬ 
ting  the  first  block  number  of  II0  Rather  than  subtracting  the  Block 
number  of  IV,  the  difference  obtained  was  instead  checked  against  all 
all  possible  block  numbers  of  IV  to  see  if  it  was  listed*  Since  if  the 

block  number  did  exist  only  one  could  exist*  This  procedure  cut  the 
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possible  number  of  block  comparisons  by  a  factor  of  eighteen  and  left 
only  5800  combinations  to  examine  rather  than  105,  000,  If  the  block 
did  exist  all  the  codes  within  this  (I  II  III  IV)  block  grouping  were 
then  given  the  parity  test.  If  the  IV  block  number  did  not  exist  II  was 
stepped  ahead  one  block  number  and  the  process  repeated.  After  II 
had  cycled  through  its  18  block  numbers,  III  was  then  stepped  ahead 
one  block  number  and  the  entire  process  repeated.  Similarly  when  III 
had  cycled  its  18  block  numbers,  I  was  stepped  ahead  and  the  process 
repeated.  To  avoid  the  A2  operation  redundancy,  whenever  I  was 
stepped  up  one  block  number  II I  was  started  from  this  same  block 
number  rather  than  at  the  first  block  number.  This  avoided  the  block 
number  combinations  (I  III) ,  and  (III  I)  both  appearing  except  when 
block  I  was  equal  to  block  III,  This  block  number  cycling  allowed  a 
vantage  breakpoint  for  partitioning  the  program  into  suitable  size  in¬ 
crements  for  computer  run  times. 

The  parity  check  was  made  by  adding(Modulo  2)  a  code  from  I 
to  a  code  from  III  to  a  code  from  II.  This  resultant  was  then  checked 
to  see  if  it  existed  as  a  code  in  the  list  from  the  IV  block.  If  it  did 
exist  it  was  complemented  and  then  sent  to  be  unpacked  for  the  necessary 
and  sufficient  like,  unlike  test*  If  it  did  not  exist  another  set  of  codes 
was  sent  in  and  the  process  repeated  until  all  possible  code  combinations 
in  (I  II  III)  were  exhausted. 

The  methods  used  in  the  unpacking  and  like-unlike  subroutines  are 
quite  straight  forward  and  will  not  be  amplified  here,  although  it  is  to  be 
noted  that  the  like-unlike  subroutine  is  not  limited  in  the  lengths  of  the 
codes  it  can  test,  up  to  the  machine  storage  size,  whereas  all  the  other 
subroutines,  since  they  use  packed  words,  do  have  quite  restrictive 
code  length  limitations. 
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This  program  was  first  used  to  make  an  exhaustive  search  for 
codes  of  length  10  for  test  purposes  since  all  kernels  of  this  length 
were  known.  After  this  satisfactory  checkout  the  exhaustive  search  for 
codes  of  length  26  was  undertaken. 

During  the  first  part  of  the  computer  run  for  n=26  Dr.  Golay  in¬ 
formed  the  writer  by  personal  correspondence  that  he  had  discovered  a 
kernel  of  length  26,  using  a  Mby-handn  technique.  ^  9  ^  The  exhaustive 
search  revealed  only  a  transformation  of  the  code  pair  that  Dr.  Golay 
had  discovered.  Therefore  it  was  proved  that  only  one  complementary 
kernel  of  length  26  existed  and  that  it  was  (ignoring  allowable  trans¬ 
formations) 

a=oioo noi  uioioii  noon  ioio 
B= 1011001000011111  1  100  11  10  10  . 
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Supplementary 
Characte  r  is  tic 
Number 

0 

20 
1  10 
1 

1  1 
21 
101 
111 
121 
2 

12 
22 
102 
1  12 
122 
13 
103 
123 

Total 


Number  of  13 
bit  codes  with 
nine  ones 

Number  of  13 
bit  codes  with 

seven  ones 

14 

70 

14 

70 

28 

140 

44 

112 

44 

1  12 

44 

1  12 

22 

56 

88 

224 

22 

56 

14 

70 

98 

196 

14 

70 

* 

28 

140 

70 

56 

70 

5  6 

44 

112 

22 

56 

35 

8 

715 

1  716 

Number  of  codes  within  each 
supplementary  block  for  n=26 


F igure  7.  2 
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CHAPTER  VIII 


A  PARTIAL  SEARCH  FOR  KERNELS  OF  LENGTH  34 
The  next  possible  kernels  for  investigation  after  length  26  were 
those  of  length  34.  As  will  be  shown  in  this  chapter  an  exhaustive 
search  for  kernels  of  length  34  was  not  feasible.  However,  an  impor¬ 
tant  subset  of  possible  kernels  was  exhaustively  searched  and  a  scheme 
which  can  be  adapted  to  the  general  search  was  programmed.  No  ker¬ 
nels  were  found  in  this  partial  search.  To  use  the  same  method  of 
attack  as  that  used  on  n=26  did  not  seem  too  feasible,  both  because  of 
the  length  of  time  involved  in  the  search  and  because  there  would  be  no 
guide  to  searching,  beyond  random  picking  of  possibilities  until  all 
possibilities  were  exhausted.  The  length  of  time  involved  for  an  ex¬ 
haustive  search  of  n=34  was  estimated  fairly  reasonably  by  checking 
the  number  of  possible  codes  of  length  26  against  the  number  of  pos¬ 
sible  codes  of  length  34,  since  the  time  for  the  26  search  was  known 

with  reasonable  accuracy.  For  n=26  there  were  71  5x7 1  5x1  7 1  6x1  7 1  6  7^ 

12 

1.5x10  code  possibilities  at  one  point  in  the  screening  process.  At 
the  same  point  in  the  screening  process  for  n=34  codes  there  were 
1  23  76x1 23  76x1  9448x1  9448  ^  50  8x1  0^  possibilities,  which  is  a  ratio  of 

4 

4x10  to  one  as  compared  to  n=26.  A  time  estimate  from  the  26 

4  6 

length  code  to  the  34  length  code  was  therefore  75x4x10  =3x10  hours. 
This  length  of  time  was  of  course  not  practical  for  a  search  of  an  ex¬ 
haustive  nature. 

The  use  of  the  cyclic  complementary  property  as  a  screen  was 

pointed  out  by  Dr.  Golay  as  highly  practical  because  some  reasonable 

19 

guess  could  be  made  as  to  which  cyclic  sets  might  contain  a  kernel. 

This  would  lead  to  a  hybrid  "by-hand"  and  computer  search.  It  would 

also  be  easier  to  document  the  areas  searched  by  the  cyclic  method  as 

compared  to  the  supplementary  method  and  would  therefore  avoid  the 
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duplication  of  effort  by  future  n=34  length  searchers*  As  will  be  dis¬ 
cussed  at  some  length,  there  is  in  addition  an  even  more  important 
reason  leading  to  the  decision  to  substitute  the  cyclic  complementary 
property  for  the  supplementary  property  as  a  screen. 

The  generation  of  cyclic  sets  which  are  to  be  utilized  in  the  search 

for  new  kernels,  is  accomplished  in  a  somewhat  unusual  fashion.  This 

method  of  classifying,  puts  within  the  same  set  all  possible  numbers 

which  are  cyclic  permutations  of  one  another,  as  is  to  be  expected; 

however,  there  are  also  contained  within  the  same  set  all  numbers 

which  are  formed  by  the  removal  of  every  other  bit  in  a  cyclic  fashion 

until  all  bits  are  used.  As  an  example  the  five  bit  number  11010  will 

be  used  to  show  all  the  numbers  contained  within  its  set.  First, 

selecting  every  other  bit,  starting  with  the  first  bit,  and  repeating  on 

the  result  gives:  1.  11010 

2.  10011 

3.  10101 

4.  11100 

1.  11010 

and  second,  by  permuting  the  bits  of  each  of  these  in  the  normal  cyclic 

permutation  manner  yields 

1.  11010,  01101,  10110,  01011,  10101 

2.  loon,  nooi,  moo,  oino,  ooiu. 

It  is  not  necessary  to  generate  the  cyclic  permutations  from  3  and  4 
since  they  are  contained  within  1  and  2. 

Counting  cyclically  the  number  of  likes  for  each  spacing  gives  the 
code  number  133  1  for  1  and  the  code  number  3113  for  20  These 
code  numbers  start  with  a  spacing  of  1  and  include  up  through  a  spacing 
of  n- 1  or  four  in  this  particular  example.,  The  code  number  3113 
means  that  at  a  spacing  of  one  there  were  3  likes  in  10011  counted 
cyclically,  at  a  spacing  of  two  there  was  1  like,  at  a  spacing  of  three 
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there  was  1  like,  and  at  a  spacing  of  four  there  were  3  likes.  The 
example  just  given  is  somewhat  trivial  since  all  binary  numbers  of 
length  5  with  3  ones  fall  into  the  same  set.  However,  this  is  not  the 
case  for  n=26  as  there  are  7  sets  for  codes  of  length  13  with  9 
ones  and  76  sets  for  codes  of  length  13  with  7  ones. 

An  examination  of  the  kernel  of  length  26,  shows  that  its  sequence 
quadruple  of  codes  came  not  from  four  different  cyclic  sets  as  one 
would  expect  but  came  instead  from  just  two  cyclic  sets. 

Using  the  kernel  form 

A=01 00 1101 11 10 101 11 100 11 1010 

B  =  10 1 1  001  00001  1 1 1 1  1 1  00  1 1 1  0 10  one  can  decompose  this  into 

1=00101 11110111 
11=00101  10011  101 
111=  1  101001110111 
IV=  0  0101  11  100010 

Using  I  to  form  its  cyclic  set  gives: 

1.  0010111110111 
2.  0111111001101 

3.  0111011111010 

4.  0101100111111 

1.  0010111110111 

Now  take  3  and  write  it  twic  e  ,  this  gives  01  1101  11(11010011  101  1  1) 

1 1010.  Picking  off  from  the  ninth  to  the  twentieth  second  bit  gives  a 

code  which  is  the  same  as  III  above  Therefore  I  and  III  are 

members  of  the  same  cyclic  set.  Similarly  take  II  and  form  its 

cyclic  set;  this  gives' 

1.  0010110011101 
2.  0110111001010 

3.  0111000101011 

4.  0100001110111 

5.  0001111100101 
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6.  001 1011011100 

7.  010111001 1010 

8.  0010100111011 
9.  0110101000111 

1  0.  0111011100001 
1  1.  0101001111100 
1  2.  00011  101  101  10 
1.  00101 10011101 

To  make  a  comparison  between  II  and  IV  the  complement  of  IV  will 
have  to  be  used  in  order  that  it  will  have  7  ones  as  does  II. 

IV=  1101000011101 

Writing  4  from  this  list  of  12  twice  gives  0100001  1101(110100001  1  101)11. 
Picking  off  from  the  twelfth  to  the  twenty  fifth  bit  gives  TV,  which  in¬ 
dicates  that  II  and  IV  are  from  the  same  cyclic  seta 

A  check  of  kernels  of  length  10  to  see  if  their  quadruples  came 
from  cyclic  sets  in  pairs,  as  the  26  kernel  did,  proved  that  both  ker¬ 
nels  of  length  10  did  in  fact  come  in  pairs.  This  was  meaningless 
however,  since  as  was  seen  in  the  example  of  length  five  with  three 
ones,  there  was  only  one  cyclic  set  possible  for  each  pair. 

In  the  hope  that  if  the  kernel  existed  in  the  n=34  case,  it  would 
also  be  formed  with  I,  III  from  one  cyclic  set  and  II,  IV  from  another 
cyclic  set,  the  program  was  modified.  This  change  along  with  giving  an 
exhaustive  search  for  code  possibilities  taken  two  at  a  time  from  cyclic 
sets  would  also  give  an  estimate  of  the  computer  run  time  for  the  ex¬ 
haustive  search  in  general  using  the  cyclic  complementary  property. 

This  was  the  primary  reason  for  the  change  in  the  program. 

Eight  octal  digits  were  used  to  designate  the  cyclic  complementary 
number  for  each  of  the  12,376  possible  half  codes  of  length  17  with 
11  ones  and  for  each  of  the  19*448  possible  half  codes  of  length  17 
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Flow  Chart  for 

Cyclic  Complementary  Designation. 
FIGURE  8,  1 
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with  10  ones.  Only  eight  digits  were  necessary  as  each  cyclic  count 
of  spacing  i  actually  counts  both  the  i  spacing  and  the  n+l-i  spacing 
simultaneously.  This  was  demonstrated  by  step  8  in  the  proof  of 
Theorem  6.  2.  The  example  of  the  code  of  length  five  given  earlier  in 
this  chapter  showed  this  symmetric  property,  as  the  code  numbers  were 
1331  and  3113.  The  first  two  digits  of  these  cyclic  complementary 
designators  contain  all  the  information  available. 

The  revision  made  to  the  computer  program  was  not  too  extensive 
as  only  blocks  2,  3,4,5  in  Figure  7.1  were  modified  to  handle  the 
cyclic  complementary  property  rather  than  the  supplementary  property. 
Of  course  all  blocks  were  modified  to  handle  codes  of  length  34  rather 
than  codes  of  length  26. 

One  difficulty  encountered  was  the  lack  of  computer  memory  for  the 
large  blocks  of  data,  and  it  was  necessary  to  store  the  cyclic  sets  on 
magnetic  tape  rather  than  in  the  main  memory  core  as  the  26  case  was 
handled.  The  cyclic  sets  were  then  called  in  two  at  a  time  for  checking. 
This  slowed  the  search  down  a  little,  but  by  no  more  than  1%  of  the 
total  time.  Figure  8,  1  is  a  block  diagram  of  the  process  of  obtaining 
the  cyclic  complementary  numbers  and  it  will  be  worthwhile  to  remember 
the  memory  size  restriction  when  reviewing  the  procedure  used. 

A  general  description  of  Figure  8.  1  will  show  the  differences  in 
the  method  used  in  the  search  for  n=34  kernels  as  compared  to  that 
used  for  n=26  kernels.  The  procedure  in  the  n=34  data  formation  was 
identical  with  that  of  n=26  as  far  as  the  initial  possible  number  genera¬ 
tion  was  concerned.  (Block  l  in  Figures  7.  1  and  8.  1)0  There  was 
an  immediate  departure  from  the  old  method  in  blocks  2,3,  and  4 
where  one  of  the  17  bit  numbers  with  11  ones  was  shifted  left  17  bits 
and  then  added  to  itself.  This  formed  a  34  bit  number  which  was 
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actually  the  17  bit  number  written  twice  in  time  sequence.  A  mask 
was  then  put  over  the  last  17  bits  and  this  masked  number  was  used, 
with  an  equality  search  of  the  12,3  76  generated  17  bit  numbers,  to 
identify  this  particular  number  in  the  list*  After  the  number  equal 
to  the  reference  number  was  located  in  the  generated  list  it  was  tagged 
a  one,  indicating  cyclic  set  one.  The  whole  34  bit  number  was  then 
shifted  right  one  bit  and  a  mask  used  to  again  pick  off  the  last  17  bits. 
The  search  through  the  list  and  tagging  with  the  group  number  was  re¬ 
peated.  After  cycling  through  all  17  possible  cyclic  codes  and  tagging 
each,  the  original  34  bit  number  was  now  repeatedly  long  right  shifted 
one  bit  and  then  right  shifted  one  bit.  A  long  right  shift  saves  the  bit 
which  is  pushed  off  the  A  register  and  stores  it  in  the  Q  register, 
while  the  right  shift  just  pushes  the  extra  bit  off  the  end  of  the  A  reg¬ 
ister.  Therefore  every  other  bit  of  the  34  code  was  saved  in  Q  and 
this  formed  a  new  1  7  bit  code  within  the  same  complementary  cyclic 
set.  This  code  was  then  formed  into  a  double  length  code  and  the  list 
of  all  possible  17  bit  numbers  with  11  ones  was  searched  and  tagged 
with  a  one  and  the  process  which  was  first  described  to  generate  cyclic 
permutations  was  repeated.  This  process  of  first  forming  17  cyclic 
permutations  of  a  code  and  then  taking  every  other  bit,  was  continued 
until  the  possible  code  which  was  originally  operated  upon  reappeared. 
Then  cyclic  set  2  was  started  and  tagged  accordingly,  and  so  on 
through  98  sets.  The  19,448  codes  of  length  17  with  7  ones  were 
classified  in  exactly  the  same  manner  and  formed  150  cyclic  comple¬ 
mentary  possibility  sets. 

After  each  series  of  sets  were  formed  a  mask  was  used  to  search 

the  set  tag  numbers  and  pick  off  each  set  for  transfer  to  magnetic  tape. 

To  insure  that  the  pick-off  was  exhaustive,  a  zero  was  inserted  on  the 
list  in  place  of  each  removed  tagged  number. 
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The  sets  were  now  ready  to  be  called  back  into  the  computer  for  use 
in  the  search.  Normally  only  one  set  was  called  in  at  a.  time  because 
the  pair  from  the  I,  III  set  was  checked  against  all  the  II,  IV  sets 
before  a  new  I,  III  set  was  called  into  the  computer^  The  I,  III  sets 
contained  the  11  weight  codes  while  the  II,  IV  set  contained  the  7 
weight  codes.  When  a  set  was  called  into  the  computer  each  cyclic 
permutation  subset  of  17  possible  codes  was  given  a  cyclic  comple¬ 
mentary  number  to  categorize  the  entire  subsetc  Each  group  contained 
either  4  or  8  of  these  cyclic  subsets.  These  subset  numbers  were 
formed  by  counting  the  likes  for  the  I,  III  codes  and  the  unlikes  for  the 
II,  IV  codes  for  the  eight  unique  spacings,  The  codes  within  the  subsets 
were  checked  against  each  other  only  if  the  likes  equaled  the  unlikes  of 
the  cyclic  complementary  subset  numbers  for  all  spacing,  or  Lj  + 


Lm  =  un 


+  u 


IV' 


Figure  8.  1  is  joined  to  Figure  70  1  just  beyond  block  5  and  takes 
the  place  of  blocks  1,2, 3, 4, 5  in  Figure  7.1.  The  same  feedback  paths 
shown  in  Figure  7.  1  to  insure  the  checking  of  all  codes  within  a  set 
and  to  feed  in  a  new  set  after  a  check  has  been  completed,  are  still  in 
operation  for  the  same  purpose.  Appendix  VI  contains  the  computer 
program  for  Figure  8.  I* 

The  total  run  time  for  the  n=34  search,  taking  I,  III  from  the 

same  cyclic  set  and  II,  IV  from  another  cyclic  set  was  approximately 

15  hours.  No  codes  were  found.  An  extrapolation  from  this  run  time 

5 

to  the  exhaustive  run  time  for  n=34  gives  15x150x98  ^2.2x10  hours, 

which  although  one  order  of  magnitude  less  than  the  supplementary 

method,  is  still  not  reasonable. 

A  second  partial  search  for  n=34  kernels  was  attempted  after 

observing  that  the  kernel  of  length  1  0 
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A=0 10 1000011 
B=0000 1001 10 


with  its  center  quad  removed  yields 

Ax=  01010011 
B1=000001 10 

which  is  a  complementary  pair  of  length  8.  Breaking  A1,  B'  into  its 
(I  II  III  IV)  configuration  gives 

1=0001 

11=1011 

111=0001 

IV=0100 

where  (I,  II)  and  (III.,  IV)  are  both  complementary  pairs 

An  attempt  was  made  to  draw  an  exact  parallel  from  this  observa- 

tion  to  codes  of  length  34.  One  possible  solution  to  equation  2.5  which 

determines  the  number  of  ones  necessary  in  each  of  the  kernels  of 

length  34  is  (21,18).  If  this  same  equation  is  applied  to  codes  of 

length  17,  two  of  the  possible  solutions  are  (11,10)  and  (11,7).  The 

(11,  10)  pair  when  interlaced  will  total  21  ones,  and  similarly  the 

(11,  7)  pair  when  interlaced  will  total  18  ones.  Whereas  the  kernel  of 

length  10  with  the  quad  removed  had  a  sequence  quadruple  form  with 

codes  of  length  4,  the  34  code  with  the  quad  removed  would  have  a 

sequence  quadruple  form  with  codes  of  length  16  Appendix  II  gives 

a  list  of  all  possible  codes  of  length  16  with  weight  10  in  the  A  code 

and  weight  6  in  the  B  code0  This  list  does  not  include  the  results  of 

the  time  reversal  transformations  which  have  to  be  included  since  weight 

is  invariant  under  T ±  ,  T2  and  T,  It  was  observed  that  if  these  codes 

of  length  16  were  interlaced  to  form  codes  of  length  32  and  then  sep- 

era.ted  in  the  middle  and  the  quad  A=-ll,  B  =  10  inserted,  all  conditions 

as  far  as  the  number  of  ones  required  in  both  interlace  codes  of  17  to 
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form  codes  of  length  34  would  be  satisfied.  An  example  might  serve 

to  clarify  the  preceding  statements, 

C=0101 1111110001 10 
D=10010011000010i0 

is  a  complementary  pair  of  length  16„  Using  the  interlace  method  to 

form  a  pair  of  length  32  gives 

A=01 100011 101011 11 1010000001 101100 
B= 001 101 101 11 1 10101 11 10101001 1 1001  . 

Splitting  this  A,B  pair  in  the  center  and  adding  the  quad  to  this  center 

give  s 

A1“01  100011  10101  111  11  1010000001101100 

B1=00 1101 1011111010  10  1111010100111001  . 

Breaking  this  pair  into  standard  (I  II  III  IV)  form  gives 

1=0101 11 1111 1000110 
11=01010000111001001 

111=0101  11  1  11  1  1000110 
IV=10101111000110110  . 

I  has  11  ones.  Interlacing  I  with  II*  which  has  7  ones,  yields 
a  total  of  18  ones.  Similarly  III  has  11  ones,  and  IV  has  10  ones. 
When  these  are  interlaced  the  code  of  length  34  has  21  ones.  In  the 
actual  computer  run  the  codes  of  length  16  with  6  ones  were  converted, 
by  complementing,  into  codes  with  10  ones.  Therefore  both  the  A  and 
B  codes  contained  10  ones.  An  exhaustive  list  of  these  and  their  time 
reversals  was  generated.  Each  code  of  length  16  was  then  split  in  the 
center  and  a  zero  inserted  for  the  ninth  bit,  forming  an  over  all  code  of 
length  17  with  10  ones,  A  second  list  was  then  formed,  using  a  one 
rather  than  a  zero  for  the  inserted  ninth  bit,  forming  codes  of  length  17 
with  11  ones. 

These  two  lists  of  numbers  were  then  fed  into  the  program  of 
Figure  7.  1  at  block  two,  and  the  regular  supplementary  check  program 


was  used.  The  codes  with  11  ones  were  used  for  (I,  III)  and  the 
codes  with  10  ones  were  used  for  (11*  IV)Q  When  these  codes  were 
combined  in  the  interlace  scheme,  I  interlaced  with  II  gave  21  ones 
and  III  interlaced  with  IV  gave  18  ones.  All  four  interlace  combina¬ 
tions  (I  II  III  TV),(II  I  III  "JV),( I  II  IV  III)  and  (II  I  TV  III)  were 
run  exhaustively.  The  run  time  was  approximately  three  hours  and  no 
codes  were  found  These  partial  searches  for  n=34  are  documented 
here  because  they  represent  two  fairly  obvious  approaches  to  the  prob¬ 
lem,  and  it  would  be  very  wasteful  of  time  for  someone  to  duplicate  this 
effort.  The  program  as  given  in  the  appendix  can  with  a  slight  modifica¬ 
tion  be  used  for  an  exhaustive  search  for  kernels  of  length  34  or  for  a 
partial  search  based  upon  other  "judicious  guesses" 


no 


CHAPTER  DC 
CONC  LUSIONS 


The  study  which  has  been  presented  here  had  two  principal  objec¬ 
tives.  The  first  objective  was  to  make  an  exhaustive  search  for  kernels 
of  length  26  to  determine  whether  they  existed,  and  if  they  did  exist, 
to  determine  the  number  of  such  kernels*  The  second  objective  was  to 
develope  a  methodical  scheme  of  code  decomposition  from  composite 
codes  back  to  their  generating  kernels.  We  shall  summarize  in  this 
chapter  the  ways  in  which  each  of  the  major  portions  of  this  dissertation 
is  related  to  the  accomplishment  of  these  objectives.  Also  a  number  of 
suggestions  for  further  research  will  be  presented. 

The  operations  group  and  its  family  of  theorems  were  originated  to 

formalize  the  elimination  of  possible  redundancies  in  the  search  for  new 

* 

kernels.  The  predetermination  of  the  number  of  ones  in  each  of  the 
sequence  quadruple  vectors  was  a  vital  screen  that  eliminated  many 
such  redundancies.  Also,  the  starting  of  both  the  I  and  III  vectors 
from  the  same  possible  supplementary  code  groups  eliminated  an  ad¬ 
ditional  half  of  the  possible  codes.  The  theorems  in  the  operations 
group  chapter  also  were  useful  in  the  code  decomposition  part  of  the 
problem,  since  they  show  the  necessary  and  sufficient  conditions  for 
the  (I,  II)  and  (III,  IV)  pairs  to  be  complementary*  Theorem  3.4, 
which  states  that  if  the  sequence  quadruple  form  is  (I  II  I  TTJ  then 
(I,  II)  must  be  a  complementary  pair,  is  most  helpful  since  the  process 
of  checking  for  complementarity  is  very  tedious  by  hand  methods. 

Hamming  distances  were  the  next  major  topic  considered.  They 
were  found  to  be  useful  in  a  number  of  the  proofs  in  the  Hamming  vector 
chapter.  They  were  also  a  great  aid  in  checking  hand  decomposition, 
since  it  was  easy  to  determine  whether  D(I*  II)=D(III ,  IV)=n /4  for  compo¬ 
site  codes  or  whether  they  had  the  characteristic  distance  in  the  case  of 
kernels.  Ill 


Hamming  vectors  were  the  most  powerful  tool  developed  in  this 
paper  for  code  decomposition,,  The  theorems  cover  all  known  possible 
decompositions  and  also  give  strong  support  to  the  conjecture  that  ker¬ 
nels  of  length  50  do  not  exist*  The  Hamming  vector  concept  should 
allow  a  fuller  understanding  of  complementary  sequences  in  general, 
due  to  its  characteristic  qualities  such  as  anti- symmetry,  and  also  be¬ 
cause  there  are  only  four  possible  configurations  of  Hamming  vectors 
for  any  code  pair.  Similarly  in  sequence  quadruple  form  it  is  easy  to 
check  whether  a  decomposed  code  may  be  complementary  by  using  its 
Hamming  vector. 

Several  exhaustive  searches  fo  r  codes  of  various  lengths  were  con¬ 
ducted  during  this  research0  The  most  important  of  these  was  the 
search  for  kernels  of  length  26.  This  search  revealed  that  only  one 
kernel  existed*  An  exhaustive  search  for  codes  of  length  16  and  20 
revealed  that  all  codes  of  these  lengths  were  generated  from  shorter 
codes  by  standard  methods.  This  gives  considerable  strength  to  the 
conjecture  that  these  generating  methods  are  the  only  ones* 

The  possibilities  for  future  research  in  this  field  of  complementary 
sequences  are  many.  Some  of  these  follow  directly  from  this  paper 
while  others  lie  quite  far  afield  from  the  ideas  exposed  here.  An  inter¬ 
esting  extension  of  the  present  work  would  be  an  exhaustive  search  for 
kernels  of  length  34a  It  is  felt  that  the  cyclic  complementary  screen 
which  was  used  in  the  partial  search  is  not  powerful  enough,  and  a 
better  screen  is  needed  before  this  search  can  be  conducted  on  an  ex¬ 
haustive  basis  using  todayss  compute  rs0  The  writer  had  considered  the 
combination  of  both  the  cyclic  and  supplementary  properties  as  screens 
in  the  same  search.  However,  a  full  study  was  not  made  of  the  gain 


which  might  be  obtained  by  this  screening  procedure*.  The  programming 
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of  such  a  screen  would  not  prove  too  difficult  and,  since  the  cyclic 
property  would  logically  be  the  first  screen  of  the  two  to  be  used,  the 
areas  of  search  could  still  be  catalogued  quite  easily. 

If  it  would  be  possible  to  determine,  before  the  search,  the  charac¬ 
teristic  Hamming  distances  D(I,  II)  and  D(II,  IV)  of  kernels,  these 
would  act  as  a  most  potent  screen  in  the  search  for  new  kernels.  It  is 
therefore  felt  that  further  research  on  characteristic  Hamming  dis¬ 
tances  of  kernels  could  prove  quite  fruitful  for  future  applications. 

It  was  conjectured  that  kernels  of  length  50,  130  and  338  do  not 
exist,  based  on  several  pieces  of  evidence,  which  although  strong  are 
not  conclusive.  Further  investigation  along  these  lines  might  offer  a 

conclusive  proof  without  the  necessity  for  an  exhaustive  search,  since 

* 

an  exhaustive  search  for  these  length  codes  approaches  the  impossible. 

During  the  investigation  of  kernels  of  length  50  the  writer  noticed  an 

oddity  which  might  be  worthy  of  further  inve stigation. 

One  kernel  of  length  10  is 
A= 1001010001 
B= 10000001 10. 

Divide  both  the  A  and  B  codes  into  segments  of  length  2,  then 
A=  1 0  01  01  00  01 

B=10  00  00  01  10. 

Now  let  the  symbols  10=A,01=A,  11-B,  and  00-B.  Forming  a  code 

of  length  50  by  this  method  of  using  symbols  of  portions  of  the  quad 
to  represent  kernels  of  length  10  gives 
Kx= AAABA 

k2=abbaa. 

This  code  pair  of  length  50  satisfies  the  parity  check,  and  also  equa¬ 
tion  2.  5  for  the  number  of  ones  in  both  full  codes  and  in  the  four  half 
codes.  Several  examples  were  checked  in  this  fashion  for  both  kernels 
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of  length  10,  and  all  satisfied  the  same  conditions  as  did  this  example, 
but  of  course  none  of  them  satisfied  the  necessary  condition  for  com¬ 
plementarity.  Along  with  helping  to  search  for  codes  of  length  50,  this 
offers  a  possibility  for  research  in  that  all  kernels  seem  to  have 
various  characteristics  in  common^  However,  the  quad  seems  more 
flexible.  Of  the  four  known  kernels,  the  quad  is  the  only  one  thus  far 
discovered  that  can  be  used  to  build  up  composite  codes  by  one  method 
of  generation  and  then  these  codes  can  be  decomposed  by  a  different 
method.  Therefore  a  study  of  kernel  characteristics  might  lead  to 
easier  ways  of  generating  new  kernels,  or  at  least  of  determining 
whether  they  exist  for  the  various  possible  lengths. 

Another  possible  field  of  research  is  cyclic  complementary  codes. 

* 

A  communications  system  is  more  likely  to  use  a  cyclic  complementary 
carrier  modulation  instead  of  straight  complementary  modulation  due  to 
the  continuous  nature  of  the  carrier.  Therefore  an  investigation  of 
these  codes  seems  quite  in  order. 

Two  other  topics  worthy  of  research  in  the  complementary  sequence 
field  are  the  various  correlation  functions  of  complementary  sequences 
and  the  frequency  spectra  of  these  sequences,  For  secure  communica¬ 
tions  purposes  it  is  important  to  have  a  code  that  is  noise-like,  A 
uniform  distribution  of  ones  and  zeros  gives  a  noise- like  appearance  to 
the  code.  For  example,  the  complementary  pair 

A=1 1010001 
B= 11011110 

is  a  typical  complementary  code  pair  of  length  8C  If  the  bits  of  the  A 
code  are  considered  in  a  cyclic  fashion  one,  two,  a.nd  three  bits  at  a 
time,  the  results  are  remarkably  uniform  as  is  seen  in  Table  9-  1. 

This  uniformity  in  output  would  make  it  extremely  difficult  for  someone 


monitoring  the  code  to  decide  if  he  had  noise  or  an  actual  signal  on  his 
receiver.  However,  it  is  seen  that  the  B  code  is  not  nearly  so  uniform 


and  would  be  easier  to  detect  as  a  signal. 
A  code 


Bits 


Number 


Two  Bits 


Number  Three  Bits 


1  4 

0  4 


1 1 

2 

1  1  1 

10 

2 

1  10 

01 

2 

101 

00 

2 

Oil 

100 

010 

001 

000 

B  code 


1  1 

4 

111 

10 

2 

1  10 

0  1 

2 

101 

00 

0 

Oil 

100 

010 

001 

000 

Noise  unifoi mity  of  a  complementary 
pair  of  length  8, 

TABLE  9  1 


Numbe  r 

1 

1 

1 

1 

1 

1 

1 

1 


2 

2 

2 

2 

0 

0 

0 

0 


An  investigation  of  a  large  number  of  codes  by  the  above  technique  would 
be  useful  in  determining  the  most  noise-like  code  pairs,  where  noise¬ 
like  is  defined  as  a  uniform  output  with  bits  taken  one,  two,  three,  etc. 
at  a  time.  This  might  in  turn  lead  to  some  general  characteristic  of 
most- noise- like  complementary  sequence  pairs.  *■’ 


Another  important  study  would  be  the  autocorrelation  function  of 
each  code  pair.  This  is  synonymous  with  the  output  of  a  code  detected 
by  its  matched  filter.  For  example,  the  autocorrelation  of  the  code  A 
given  above  is 

1  0  -3  0  -1  0  -1  8  -1  0  -1  0  -3  0  1 

where  the  outputs  in  time  are  read  from  left  to  right  The  autocorrela¬ 
tion  of  the  B  code  is 

-  1  0  3  0  1  0  1  8  1  0  1  0  3  0  -1. 

The  autocorrelation  functions  are  important  from  two  standpoints.  First, 
if  a  portion  of  the  system  should  need  to  be  shut  down  for  repairs  or 
routine  maintenance,  a  good  autocorrelation  function  would  allow  opera¬ 
tion  on  just  one  carrier*  Second,  the  frequency  spectrum  of  the  trans¬ 
mitted  signal  is  the  Fourier  transform  of  this  autocorrelation  function. 
Therefore,  one  way  to  study  the  spectrum  is  first  to  autocorrelate  and 
then  to  take  the  Fourier  transform  of  the  transmitted  signals0  The  more 
widespread  the  spectrum,  the  more  difficult  is  a  jamming  procedure. 

A  study  of  the  crosscorrelation  functions  of  complementary  sequence 
pairs  is  also  very  important  because  of  fhe  difficulty  of  isolating  the  A 
code  RF  carrier  from  the  B  code  receive  rc  It  was  pointed  out  in 
chapter  4  that  because  the  A  code  and  the  B  code  were  orthogonal,  at 
exact  match  cross  talk  would  be  no  problem.  However,  when  the  signals 
are  not  at  exact  match  cross  talk  might  be  a  large  source  of  noise.  An 
investigation  of  the  c  ros  sc  o  rre  lation  func  tions  of  the  A  and  B  code 
pairs  might  lead  to  a  class  of  complementary  sequences  where  the 
problem  of  crosstalk  is  minimized. 

Another  class  of  correlation  functions  which  are  important  are 
those  where  cyclic  errors  are  created  in  the  received  signal  detected  by 
the  matched  filter.  These  errors  could  be  caused  by  a  linear  shift  in 


the  phase  of  the  received  signal,,  There  may  be  certain  classes  of 
complementary  sequences  where  the  loss  in  coding  signal- to -noi se  ratio 
is  minimized  when  noise  is  caused  by  an  error  of  this  type.  This  par¬ 
ticular  class  of  signals  would  be  extremely  important  in  an  air  search 
radar  application. 

A  study  of  the  maximum  number  of  different  composite  codes  for 
any  length  would  be  an  important  contribution  to  the  field,  because  for 
anti- jam  reasons  it  is  very  important  to  be  able  to  change  codes,  and 
the  more  codes  with  good  characteristics  that  are  available  the  better. 

The  last  study  to  be  recommended  is  an  analysis  of  the  autocorrela¬ 
tion  functions,  frequency  spectra,  cross  correlation  functions,  and 
linear  phase  shift  correlation  functions  for  codes  of  lengths  from  100 
to  ZOO.  The  purpose  of  this  investigation  would  be  to  determine  if 
certain  kernels  are  capable  of  forming  composite  codes  of  better  char¬ 
acteristics  than  other  kernels  for  radar  applications. 

In  conclusion,  the  writer  feels  that  complementary  sequence  pairs 
will  prove  very  important  in  future  applications  to  both  radar  and  com¬ 
munications  schemes.  However,  much  research  must  be  done  to  pick 
optimal  classes  of  these  codes  for  such  applications.  This  paper  has 
offered  many  tools  to  be  used  by  future  investigators  of  complementary 
sequences.  These  tools  have  found  important  applications  in  their 
present  state,  but  additional  investigations  along  the  lines  of  operations 
groups,  Hamming  weights  and  Hamming  vectors  should  further  in¬ 
crease  their  usefulness. 
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GLOSSARY  OF  SYMBOLS  AND  TERMS 


(A,  B) 

A  complementa  ry  pair  of  binary  sequences, 
see  equation  Za  2. 

(I  II  III  IV) 

A  complementary  pair  of  binary  sequences  in 
sequence  quadruple  form,  see  equation  3.  1. 

D(U,V) 

The  Hamming  distance  of  two  binary  vectors, 
see  equation  4e  10 

H(U,  V) 

The  Hamming  vector  of  two  binary  vectors, 
see  the  first  page  of  Chapter  5e 

* 

Means  interlace  two  Hamming  vectors  in  a  pre¬ 
scribed  manner,  see  Theorem  5.11,  step  2. 

Is  transformed  into. 

A_ 

The  time  inverse  of  a  binary  sequence  or 
vector  A. 

Alte  r 

Means  to  complement  every  other  bit  of  both 
codes  of  a  complementary  code  pair. 

Complementary 

code 

A  pair  of  binary  sequences  which  satisfy 
equatipn  2„  2,  also  called  a  complementary 
sequence  pair,  a  complementary  sequence 
and  sometime  just  a  pair  or  code. 

Complement  a 
code 

Means  change  the  sign  or  take  the  complement 
of  a  binary  sequence 

Composite  code 

A  complementary  code  which  is  reducible  to  a 
shorter  code  Length  by  standard  methods. 

Cyclic  comple¬ 
mentary  code 

A  pair  of  binary  sequence  which  satisfy 
equation  6.  3. 

Hamming  Distance 

The  number  of  bits  in  which  two  binary  sequence 
differ,  notation  is  DfA,B^?  see  equation  4.  1. 

Hamming  vector 

A  binary  vector  or  sequence  formed  by  the 
modulo  2  sum  of  binary  sequences. 

Hamming  weight 

Hamming  distance  of  a  binary  sequence  or  vector 
from  the  null  vector.  It  is  the  number  of  ones 
in  a  binary  vector. 

Ke  rnel 

A  set  of  complementary  codes  which  is  irreducible. 
The  term  is  sometimes  used  for  a  member  of  the 
s  et. 

Quad 

A  kernel  of  length  2. 

Sequence 

Quadruple 

A  complementary  code  pair  expressed  as  in 
equation  30  1 0 

Supple  mentary 
code 

A  quadruple  of  binary  sequences  which  satisfy 
equation  60  1  * 

Time  inverse 

Means  to  make  the  first  bit  of  a  binary  sequence  the 
last  bit,  the  second  bit  next  to  last,  and  so  on;  also 
called  time  reverse. 

Weight 

See  Hamming  weight0 
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EXPLANATION  OF  FORMAT  OF  THE  APPENDICIES 


Appendix  I  on  matched  filte  rs  is  s elf-explanatory.  Appendix  II 
is  the  result  of  an  exhaustive  search  for  codes  of  lengh  16.  There  were 
96  code  pairs  found  by  the  computer  after  the  program  screens  had 
eliminated  most  of  the  redundancies.  The  format  of  Appendix  II  gives 
the  A  code,  the  B  code,  H(A,  B),  and  H[Ti  (A,B)]  which  signifies 
time  reversing  the  A  code.  It  is  worthwhile  to  note  that  each  code  pair 
and  the  pair  formed  by  the  time  reverse  of  A  were  decomposable. 

From  these  192  possible  decompositions  48  were  interlace,  48  were 
time  sequence,  and  96  were  2  special  method,  with  48  each  of  two 
different  kinds. 

Appendix  III  is  the  result  of  an  exhaustive  search  for  codes  of 
length  20.  The  format  of  Appendix  III  gives  the  A  code,  the  B  code, 
if  it  is  required,  the  time  reversal  operation  used  to  put  H(A,  B)  in 
standard  form,  H(A,B)  and  from  which  of  the  2  kernels  of  length  10 
the  generating  codes  came.  All  of  the  24  code  pairs  or  their  time 
inverse  were  decomposable.  There  are  12  from  each  kernel,  8  of 
each  of  these  interlace  and  the  other  4  are  time  sequence* 

Appendix  IV  is  a  diagrm  of  all  the  subgroups  of  the  operations 
group.  It  starts  with  the  subgroup  of  order  1  and  works  up  to 
those  of  order  16.  An  x  in  the  row  indicates  that  the  operation  at 
the  top  of  the  column  is  a  member  of  the  subgroup* 

Appendix  V  is  the  CDC  1604  computer  program  in  AR  format 
and  also  in  machine  language  for  the  26  length  search*  There  are 
slight  modifications  in  the  various  sections  of  the  program  required 
for  other  length  codes.  Appendix  VI  is  also  in  both  AR  and  machine 
language,  and  is  for  codes  of  length  34.  This  program  parallels 
Figure  8.1,  whereas  the  program  in  Appendix  V  parallels  Fig.  7.  1. 
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APPENDIX  I 


MATCHED  FILTERS 

A  Matched  Filter  is  by  definition  a  filter  which  maximizes  the  peak 
signal- to- noise  ratio.  Deterioration  of  the  signal  wave  form  is  accepted 
in  order  to  obtain  the  desired  maximum  ratio. 

One  tool  needed  for  the  derivation  of  the  matched  filter  character¬ 
istics  is  the  Schwarz  inequality  which  is  a  special  case  of  the  Holder 
ZZ 

inequality.  One  representation  of  the  Schwarz  inequality  is 


j" 


x(w)  y  (w)  dw  |  <  J  *  x(w)  |^dw  j  |  y(w)  udw 


(A  l.l) 


It  is  to  be  noted  that  the  equality  holds  when  y(w)  is  the  complex 
conjugate  of  x(w)c  Rewriting  equation  A  l.  1  gives 


Too  Z 

I  J  x(w)  y(w)  dw  1 

1  mm  00  ' 

03  2  °°  ? 

J  x(w)  I  dw  J  |y(w)  dw 


<  l 


(A  1.2) 


The  derivation  for  the  matched  filter  given  here  closely  parallels 
that  of  Turin. 

Let  f(t)  be  a  signal  impressed  across  a  filter  whose  frequency 
characteristic  is  H(jw)c  F(jw),  the  frequency  spectrum  of  f(t),  is 
given  by  the  Fourier  transform  of  f(t)c 


F(jw)  =  J 


f(t)  e  -jwtdt. 


G(jw)  is  the  frequency  spectrum  of  the  output  after  passing  through 
the  filter 

G(jw)  -  H(jw)  F(jw), 


or  transforming  back  to  the  time  domain  the  output  signal  voltage  is 


g(t) 


0  oo 

=  J  H(jw)  F(jw)  e  ^NVtdf. 


(Al.  3) 


Since  at  some  time,  t  =  A,  g ( t )  must  be  a  maximum. 


max  g(t)  =  g ( A  )  = 


J  H(jw)  F(jw)  e  JwAdf. 


(A  1.4) 


To  complete  the  derivation  it  will  be  necessary  to  obtain  the  total 
noise  power,  and  the  total  power  in  the  signal  If  the  noise  is  assumed 
to  be  white  noise  of  N  watts  per  cycle,  the  output  noise  power  density  is 
Nq  =  N  |  H(jw)  |2, 

and  the  output  noise  power  is  therefore 


CO  n  CO 

J  NdfsJ  N[H(jw)rdf. 

—  oo  o  —  oo  [  * 

The  total  energy  contained  in  the  signal  is 


E 


J"  '2 


f  (t)  dt. 


(A  1.5) 


which  is  also 


r  oo 

E  =  Jm  |  F(jw)  |  <If ,  (A  1.6) 

by  Plancherel’s  Theorem.  It  is  to  be  noted  that  E  is  a  constant  not 
dependent  upon  the  filter  used 

The  output  power  signal- to- noi se  ratio  for  the  maximum  output  is 
the  refo  re 
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MAlil 

N 


C  oo 

|  J  F(jw)  H(jw)  eJwAdf 


N 


Dividing  this  equation  by  a  constant  will  not  effect  the  time  of  maximum 
signal  output;  therefore  dividing  by  E,  equation  A  1.  6  ,  gives 

00 


N  |  g(A)  |  2  _  j  Jm  F(jw)  H(jw)  e  JwA  df  |  ? 


F  I  H(jw)  !2df  |  |  F(jw)  I2  df  (A  1.  7) 


If  the  right  hand  side  of  equation  (A  1.7)  is  compared -to  the 
Schwarz  inequality  (A  1*2)  with  x(w)  =  F(jw)  e  y(w)  =  H(jw), 

the  two  expressions  are  the  same.  But  for  the  expression  to  be  max¬ 
imum  y(w)  must  equal  the  complex  conjugate  of  x(w),  therefore 
H(jw)  =  F(-jw)  e  'JwAs 

A  matched  filter  is  therefore  a  filter  whose  frequency  character¬ 
istic  is  the  complex  conjugate  of  the  signal  spectrum  to  which  is  matched. 
The  output  signal  for  a  given  signal  applied  to  a  matched  filter  is 
the  refo re 

oc 

gw-L  iF<jw>t2« jw<t-A)df . 

at  t  =  A  the  output  is 


CO 


This  is  the  same  as  equation  A  1.  6  which  is  all  the  energy  con¬ 
tained  in  the  signal.  Therefore  at  exact  match  all  the  energy  in  the 
pulse  is  in  the  output  signal. 
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APPENDIX  II 


Exhaustive  Search  for  Codes 
of  Length  1  6. 
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APPENDIX  V 

REMARKS  ON  COMPUTER  PROGRAM 

In  Appendix  V  the  addresses  from  60000  to  60035  are  block  1 
of  Figure  7.  1,  All  numbers  used  in  the  program  are  octal.  The  program 
as  written  is  for  codes  of  length  26,  so  most  of  the  program  deals  with 
half  length  codes  of  length  13,  Words  to  be  changed  for  different  length 
codes  are  60000  (mask  of  ones  for  half  length  code  as  determined  by 
weight),  60001  (length  of  half  code  minus  one  [upper  half],  and  length 
of  memory  word  minus  length  of  code  [  lower  half]) ,  60004  (length  of 
half  code),  60034  (last  sequence  in  code  list).  For  exactly  the  same 
reasons  respectively  change  60015,  60016,  60021  and  60032. 

Addresses  1  7000  to  17045  separate  the  sequences  into  blocks  by 
supplementary  number.  Ii\dex  register  2  must  be  entered  with  100 
before  starting.  Addresses  to  be  changed  here  are  1  7000  (mask  of  first 
3  and  last  3  bits),  1  7005  (contains  the  total  number  of  sequences  for 
both  half  length  codes),  17007  (contains  the  number  of  sequences  with 
7  ones),  17015  (contains  a  mask  for  half  length  codes),  17026  (contains 
the  number  of  sequences  with  9  ones),  1  7032,  17033  (have  for  their 
instruction  address  the  base  word  50000  plus  the  number  of  half  codes 
with  7  ones),  17034  (contains  a  mask  for  code  length). 

Addresses  10000  to  10106  are  the  repetitive  part  of  the  program. 
The  only  memory  address  needing  modification  here  for  different  length 
codes  is  10106  which  contains  SVN  a  mask  equal  to  the  one  half  code 
length. 

Addresses  7200  to  7252  contain  the  unpacking  subroutine  and  the 
like-unlike  subroutine.  7201,7207,721  5,  and  7223  all  contain  the 
half  code  length.  Similarly  7204,7207,7220,7223  and  7230  all  contain 
the  full  code  length  or  the  length  modified  by  one  or  two  to  fit  the  program 


needs . 
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The  15000  to  15100  series  is  all  possible  bit  arrangements  for 
the  first  3  and  the  last  3  bits  of  codes  of  length  13,  while  16000  to 
16100  is  a  list  of  supplementary  block  numbers  in  one  to  one  correspon¬ 
dence  with  this  list. 

The  10300  and  10400  series  respectively  contain  the  number  of 
codes  in  each  supplementary  block  for  the  half  length  codes  of  weight 
9  and  weight  7. 

16000  to  16020  is  the  final  list  of  the  18  supplementary  block 
numbers  and  is  in  one  to  one  correspondence  with  the  10300  and  10400 
series . 
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INF4 

12 

2 

10000 

INF4 

LDA 

2 

A 

75 

4 

0001  6 

SLJ 

4 

THIRD 

12 

0 

50101 

LDA 

0 

TE5 

05 

0 

00006 

ALS 

0 

6 

14 

0 

50071 

ADD 

0 

DUM5 

20 

0 

50041 

STA 

0 

INFS 

12 

3 

1  1  000 

INFS 

LDA 

3 

B 

75 

4 

0001  5 

SLJ 

4 

SECOND 

1  2 

0 

50073 

LDA 

0 

TE3 

05 

0 

00006 

ALS 

0 

6 

14 

0 

50072 

ADD 

0 

DUM6 

20 

0 

50044 

STA 

0 

INF  6 

12 

1 

10000 

INF6 

LDA 

1 

A 

75 

4 

00014 

SLJ 

4 

FIRST 

57 

1 

50074 

SIL 

1 

FAG 

57 

2 

50075 

SIL 

2 

SAG 

57 

3 

50076 

SIL 

3 

LAG 

57 

4 

50077 

SIL 

4 

GAG 

75 

0 

00013 

SLJ 

0 

BEGINEQ 

50 

0 

00000 

ENI 

0 

0 

53 

1 

50074 

REL 

LTL 

1 

FAG 

53 

2 

50075 

L1L 

2 

SAG 

53 

3 

50076 

LIL 

3 

LAG 

53 

4 

50077 

LIL 

4 

GAG 

54 

00020 

INSKIP 

ISK 

3 

20 

75 

0 

50030 

SLJ 

0 

INF  1 

155 


50053 

50054 

50055 

50056 

50057 

50060 

50061 

50062 

50063 

50064 

50065 

50066 

50067 

50070 

50071 

50072 

50073 

50074 

50075 

50076 

50077 

50100 

50101 

50102 

50103 

50104 

50105 


54 

2 

00020 

ISK 

75 

0 

50030 

SLJ 

54 

1 

00020 

ISK 

75 

0 

50030 

SLJ 

52 

1 

50073 

LIU 

53 

2 

50100 

LIL 

52 

3 

50101 

LIU 

52 

4 

50102 

LIU 

54 

3 

00000 

INFA 

ISK 

75 

0 

50007 

SLJ 

54 

2 

00000 

INFB 

ISK 

75 

0 

50007 

SLJ 

54 

1 

00000 

INFC 

ISK 

75 

0 

50007 

SLJ 

76 

0 

01111 

SLS 

50 

0 

00000 

ENI 

53 

4 

63101 

SWHE 

LIL 

75 

0 

50057 

SLJ 

12 

1 

10000 

DUM1 

LDA 

42 

2 

10000 

SCM 

42 

3 

11000 

DUM2 

SCM 

50 

4 

00020 

ENI 

64 

4 

1  1000, 

DUM3 

EQS 

75 

0 

50052 

SLJ 

00 

0 

00000 

TRENSE 

OCT 

00 

3 

77777 

1  2 

2 

100000 

DUM4 

LDA 

75 

4 

00012 

SLJ 

1  2 

3 

11000 

DUM5 

LDA 

75 

4 

0001  1 

SLJ 

12 

1 

10000 

DUM6 

LDA 

75 

4 

00010 

SLJ 

00 

0 

00000 

TE  3 

OCT 

00 

0 

00000 

00 

0 

00000 

FAG 

OCT 

00 

0 

00000 

00 

0 

00000 

SAG 

OCT 

00 

0 

00000 

00 

0 

00000 

LAG 

OCT 

00 

0 

00000 

00 

0 

00000 

GAG 

OCT 

00 

0 

00000 

00 

0 

00000 

TE4 

OCT 

00 

0 

00000 

00 

0 

00000 

TE  5 

OCT 

00 

0 

00000 

00 

0 

00000 

TE6 

OCT 

00 

0 

00000 

00 

0 

00000 

DEL 

OCT 

00 

0 

00000 

B 

EQU 

00 

0 

00000 

DOL 

OCT 

00 

0 

00000 

42 

0 

07457 

DUM7 

SCM 

22 

0 

07454 

AJP 

1  56 


2  20 
0  INF  1 
1  20 
0  INF  1 

1  TE3 

2  TE4 

3  TE5 

4  TE6 

3  0 

0  JMP2 
2  0 

0  JMP2 
1  0 

0  JMP2 
0  1111 
0  0 

4  63101 
0  INFA 

1  A 

2  A 

3  B 

4  20 
4  B 

0  INSKIP 
377777 

2  A 

4  THIRD 

3  B 

4  SECOND 
1  A 

4  FIRST 
0 

0 

0 

0 

0 

0 

0 

0 

0 

1  1000 
0 

0  PE  RE 
0  FRERE 


50106 

50107 

50110 

51000 

51001 

51002 

51003 

51004 

51005 

51006 

51007 

51010 

51011 

51012 

51013 

51014 

51015 

51016 

51017 

51020 

51021 

51022 

51023 

51024 

51025 

51026 

51027 


42 

0 

07457 

DUM10 

SCM 

0 

PERE 

22 

1 

07454 

A  JP 

1 

FRERE 

00 

0 

00000 

RETNA 

OCT 

0 

00 

0 

00000 

00 

0 

00000 

ANTER 

OCT 

0 

00 

0 

00000 

ORG 

51000 

75 

0 

00000 

INPU 

SLJ 

0 

0 

50 

0 

00000 

ENI 

0 

0 

57 

6 

51020 

S1L 

6 

CONT 

50 

6 

00004 

ENJ 

6 

4 

50 

0 

00000 

ENI 

0 

0 

10 

0 

00000 

ENA 

0 

0 

74 

0 

52031 

EXF 

0 

52031 

74 

7 

52000 

EXF 

7 

52000 

20 

0 

57777 

STA 

0 

57777 

10 

0 

60001 

ENA 

0 

60001 

61 

0 

00005 

SAL 

0 

5 

74 

5 

57777 

EXF 

5 

57777 

12 

0 

57777 

LDA 

0 

57777 

22 

0 

51006 

A  JP 

0 

/ 

14 

0 

51014 

ADD 

0 

CONSI 

20 

0 

00005 , 

STA 

0 

5 

50 

0 

00000 

ENI 

0 

0 

74 

7 

52000 

EXF 

7 

52000 

74 

7 

52003 

EXF 

7 

52003 

75 

0 

5101  5 

SLJ 

0 

ERR 

74 

7 

52005 

EXF 

7 

52005 

75 

0 

51015 

SLJ 

0 

ERR 

53 

6 

51020 

LIL 

6 

CONT 

75 

0 

51000 

SLJ 

0 

INPU 

00 

0 

00001 

CONSI 

ZRO 

0 

1 

00 

0 

00000 

ZRO 

0 

0 

55 

6 

5101  6 

ERR 

IJP 

6 

/+ 1 

75 

0 

51036 

SLJ 

0 

EREX 

74 

0 

52006 

EXF 

0 

52006 

74 

7 

52000 

EXF 

7 

52000 

75 

0 

51002 

SLJ 

0 

INPU+2 

50 

0 

00000 

ENI 

0 

0 

00 

0 

00000 

CONT 

OCT 

0 

00 

0 

00000 

50 

0 

00000 

ENJ. 

0 

0 

50 

0 

00000 

50 

0 

00000 

OUTPUT 

ENI 

0 

0 

74 

0 

52000 

EXF 

0 

52000 

57 

1 

51030 

SIL 

1 

/  +5 

55 

1 

51025 

IJP 

1 

/ +2 

75 

0 

51031 

SLJ 

0 

RTN2 

50 

0 

00000 

ENI 

0 

0 

10 

0 

60000 

ENA 

0 

60000 

61 

0 

00005 

SAL 

0 

5 

74 

0 

52031 

EXF 

0 

52031 

74 

7 

52000 

EXF 

7 

52000 

74 

5 

57777 

EXF 

5 

57777 

74 

7 

52000 

EXF 

7 

52000 

157 


51030 

51031 

51032 

51033 

51034 

51035 

51036 

51037 

51040 

51041 


55 

1 

51027 

IJP 

1 

/-I 

50 

0 

00000 

ENI 

0 

0 

50 

0 

00000 

RTN2 

ENI 

0 

0 

76 

0 

02134 

SLS 

0 

2134 

75 

0 

00000 

VERT 

SLJ 

0 

0 

50 

0 

00000 

ENI 

0 

0 

12 

0 

50106 

LDA 

0 

DUM10 

75 

0 

51032 

SLJ 

0 

VERT 

50 

5 

00010 

TERN 

ENI 

5 

10 

12 

0 

07734 

LDA 

0 

TEM2 

20 

0 

63101 

STA 

0 

63101 

75 

0 

07700 

SLJ 

0 

YEA 

12 

0 

50107 

EREX 

LDA 

0 

RETNA 

53 

5 

50110 

LIL 

5 

ANTER 

20 

5 

65000 

STA 

5 

HOPET 

51 

5 

00001 

INI 

5 

1 

57 

5 

50110 

SIL 

5 

ANTER 

50 

5 

00000 

ENI 

5 

0 

75 

0 

07400 

SLJ 

0 

ERAS 

50 

0 

00000 

ENI 

0 

0- 

HO  PET 

EQU 

END 

65000 

158 
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